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Kurzfassung 


Nanoverbundmaterialien bestehen aus einem Träger-Material, auch als Mat- 
rix-Material bezeichnet, und einem nanoskaligen Füllstoff und spielen eine 
zunehmend wichtige Rolle in vielen Bereichen des industriellen, medizini- 
schen oder des alltäglichen Umfelds. Die Attraktivität von Nanomaterialien 
beruht auf deren einzigartigen physikalischen Eigenschaften, die über jene von 
klassischen Materialien hinausgehen. So ergeben sich, beispielsweise durch 
Oberflächeneffekte dieser Stoffe, Materialeigenschaften, die im Wesentlichen 
von deren Oberflächenbeschaffenheit sowie von Form und Größe der Nano- 
strukturen abhängen. 


Eine wesentliche Herausforderung in der Entwicklung von Nanopartikeln und 
Nanoverbundmaterialien, als Nanokomposit-Materialien bezeichnet, stellt die 
Sicherstellung von Partikelgrößen und -formen, sowie deren Verteilung im 
umgebenden Matrix-Material dar. Diese Eigenschaften werden während des 
Produktionsprozesses aktiv beeinflusst, beispielsweise durch Scherkräfte in 
einem Dispergierextruder. Hierbei werden sowohl Partikelagglomerationen 
aufgebrochen, als auch einzelne Partikel zerteilt. Je nach Anwendung können 
diese Vorgänge beabsichtigt oder unerwünscht sein. So ist für ein elektrisch 
leitfähiges Polymer mit Kohlenstoffnanoröhren ein Netzwerk vereinzelter Na- 
noröhren notwendig, die dabei selbst mit möglichst großem Aspektverhältnis 
innerhalb des Materials vorliegen sollen, um die Bildung eines leitfähigen 
Netzwerks zu fördern. Die Prozesskontrolle ist konventionellerweise nur nach- 
träglich möglich, indem einzelne Ausschnitte von wenigen hundert Mikrome- 
tern Größe mit elektronenmikroskopischen Verfahren untersucht werden. Das 
führt zu einer Vielzahl von Problemen, welche die insbesondere mittelständi- 
schen Hersteller von Spezialkunststoffen schwer bewältigen können: Einer- 
seits ergibt sich durch aufwändige Laboruntersuchung, die oft außer Haus 
stattfindet, eine große Zeitverzögerung bei der Weiterentwicklung der Produk- 
tionsverfahren, andererseits liefern bildgebende Verfahren mit mikroskopi- 
scher Auflösung nur Informationen über die Materialbeschaffenheit eben jenes 
kleinen untersuchten Ausschnitts. Darüber hinaus unterscheidet sich die Größe 
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von Nanopartikeln und ihren Agglomeraten teilweise um mehrere Größenord- 
nungen, was die Charakterisierung vor Herausforderungen stellt. Eine weitere 
große Hürde stellen die Investitions- und Unterhaltungskosten für ein speziali- 
siertes Materiallabor dar. Viele Hersteller verzichten daher auf mikroskopi- 
sche Untersuchung und optimieren ihre Prozesse unter Beobachtung von mak- 
roskopischen Materialparametern. Eine tiefere Kenntnis der Nano- und 
Mikrostruktur wird dabei nicht erlangt, und weiteres Entwicklungspotential 
bleibt begrenzt. 


Im Rahmen dieser Arbeit wird die Anwendung von optischer Kohä- 
renztomographie (OCT) auf die Nanokomposit- und Nanopartikelcharakteri- 
sierung vorgestellt. OCT ist ein Verfahren zur dreidimensionalen Bildgebung 
mit mikroskopischer Auflösung. Mit einem typischen Messbereich von eini- 
gen Kubikmillimetern stellt OCT eine Erweiterung von Sonographie, Compu- 
tertomographie und Magnetresonanztomographie in Richtung kleinerer Auflö- 
sungen und Messbereiche dar. OCT basiert auf einer Abtastung der Probe mit 
einem Laserstrahl und dem Überlagerungsempfang des von der Probe rückge- 
streuten Lichts. Übliche Anwendungsbereiche sind weitestgehend in der Me- 
dizintechnik zu finden, insbesondere in der Augenheilkunde. 


Zur Charakterisierung von Nanokomposit-Materialien wird in dieser Arbeit 
ein multiskaliger Ansatz verfolgt. Während Partikel und Agglomerationen im 
Mikro- und Millimeterbereich bildgebend analysiert werden, ergänzt eine 
streumodell-basierte Auswertemethode die Analyse hin zu kleinsten Partikel- 
größen im Nanometerbereich. In diesem Zusammenhang werden der Durch- 
messer und das Aspektverhältnis von sphärischen und zylinderförmigen Parti- 
keln in homogen verteilten, monodispersen Partikelproben mit einer 
Genauigkeit von wenigen Nanometern analysiert und der Dispergierungszu- 
stand von Nanokompositen wird sowohl offline als auch inline, d.h. während 
des Produktionsvorgangs, bestimmt. Hierbei zeigt sich die Robustheit dieses 
Messverfahrens, welches in der Lage ist, hinreichend große Materialproben zu 
untersuchen, um zuverlässige Rückschlüsse auf die globalen Materialparame- 
ter zu erlauben. Weiterhin ist dieses zerstörungsfreie Verfahren für flüssige 
und feste Materialproben geeignet und kommt ohne Probenvorbereitung aus. 
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Fiir Anwendungen von OCT abseits der Medizintechnik ergeben sich andere 
Schwerpunkte: Auflösung, Empfindlichkeit und Messgeschwindigkeit werden 
ergänzt durch Anforderungen hinsichtlich Robustheit, geringen Systemkosten, 
Parallelisierbarkeit und kleiner Baugröße. Konventionelle freistrahl- oder fa- 
seroptische Systeme sind bauartbedingt limitiert bezüglich minimaler Größe 
und Kosten und benötigen eine exakte geometrische Justage und Einstellung 
der Polarisation. Temperaturschwankungen und Vibrationen verändern opti- 
sche Kopplungseigenschaften und Polarisationszustände, weshalb OCT- 
Systeme derzeit oft in geschützten Umgebungen und von geschultem Personal 
zu bedienen sind. 


Die photonische Integration von Kommunikations- und Messsystemen bietet 
robuste, günstige und kompakte Alternativen zu herkömmlichen Aufbauten. 
Durch photolithographische Herstellungsverfahren kann eine Vielzahl von 
Systemen gemeinsam auf einem Chip untergebracht werden. Im Rahmen die- 
ser Arbeit werden siliziumphotonisch integrierte OCT-Systeme untersucht. 
Dort sind sowohl Interferometer als auch differentielle Photodetektoren auf 
einem gemeinsamen Chip untergebracht. In Systemcharakterisierungen wird 
gezeigt, dass bereits kleine Probenreflektivitäten von —64 dB detektierbar sind. 
Die integriert-optischen Systeme ermöglichen damit die Untersuchung von 
biologischen und technischen Proben. 


Kapitel 1 führt in die Problemstellung ein, indem zunächst auf Nanomateria- 
lien, deren Anwendungsgebiete und auf die Notwendigkeit von geeigneten 
Charakterisierungsverfahren eingegangen wird. Darauf folgend wird das 
OCT-Messprinzip eingeführt und schließlich werden die Möglichkeiten der 
Siliziumphotonik zur Integration von Messsystemen beschrieben. 


Kapitel 2 geht auf die Grundlagen der OCT und der polarisationsauflösenden 
OCT ein. Dabei wird das Messsignal eines OCT-Systems mathematisch be- 
schrieben und Systemeinflüsse auf Auflösung, Messbereich und Empfindlich- 
keit werden untersucht. Grundlagen der optischen Polarisation und von polari- 
sationsempfindlichen OCT-Systemen werden vorgestellt und beschrieben. 
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Kapitel 3 stellt Ergebnisse von Experimenten in denen OCT zur multiskaligen 
Untersuchung von Nanokompositen verwendet wird vor. Dabei werden Nano- 
partikel und Partikelagglomerate in Größenordnungen von Nano- bis Millime- 
tern untersucht. Große Strukturen werden klassisch bildgebend analysiert, 
während Nanopartikel mittels eines optischen Streumodells untersucht wer- 
den. Die entwickelte Messmethodik wird dabei sowohl in einer Laborumge- 
bung als auch im Produktionsumfeld angewandt. Eine eigens entwickelte 
OCT-Messsonde erlaubt dabei Messungen während der Nanokomposit- 
Produktion. 


Kapitel 4 beschreibt eine Erweiterung des in Kapitel 3 vorgestellten Streumo- 
dells auf Polarisationseffekte. Im Vergleich mit elektromagnetischen Streusi- 
mulationen wird mit polarisationsauflösender OCT gezeigt, dass neben der 
Partikelgröße auch die Form von Nanopartikeln bestimmt werden kann. 


Kapitel 5 geht auf die siliziumphotonische Integration von OCT-Systemen ein. 
In diesem Zusammenhang werden zwei Systeme mit integriertem Interferome- 
ter und integrierten Photodioden vorgestellt. Die Systeme nutzen entweder 
einen internen oder einen externen Referenzpfad, wodurch unterschiedliche 
Abstände zwischen Chip und untersuchter Probe ermöglicht werden. Eine um- 
fassende Analyse spezifischer Störsignalquellen zeigt, welche Aspekte beim 
Design integriert-photonischer Systeme besonders berücksichtigt werden müs- 
sen. Dreidimensionale Aufnahmen von biologischen und technischen Objek- 
ten belegen die Funktionalität der Systeme. 


Kapitel 6 fasst die Ergebnisse dieser Arbeit zusammen und gibt einen Aus- 
blick auf zukünftige Entwicklungen. 


Die im Rahmen dieser Arbeit entstandenen Ergebnisse wurden in mehreren 
wissenschaftlichen Fachzeitschriften publiziert sowie auf internationalen Kon- 
ferenzen vorgetragen. Eine Auflistung aller Veröffentlichungen ist auf Sei- 
te 181 ff. zu finden. 
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Preface 


Nanocomposite materials consist of a host material and a nanoscale filler ma- 
terial and play an increasingly important role in various fields of applications, 
reaching from medical diagnostics and therapy to industrial applications and 
everyday necessities. The attractiveness of nanomaterials arises from their 
unique physical properties, which go beyond those of classical materials. Sur- 
face effects, for instance, lead to material properties, which depend essentially 
on the nanostructure of the surface. 


A major challenge in the development of nanoparticles and nanocomposite 
materials is to achieve the desired particle sizes and shapes, and the dispersion 
in the host material. These properties are actively influenced during the pro- 
duction process, e.g., by application of shear forces in a compounding extrud- 
er. Thereby, particle agglomerates are split, and single particles are fractured 
or separated. Depending on the application, these processes are either intended 
or unwanted. In the case of, e.g., an electrically conductive polymer with car- 
bon nanotubes, a network of dispersed nanotubes is necessary, while the single 
tubes preferably keep their large aspect ratio and thus enable the formation of 
a conductive mesh. Process control is only possible afterwards, using an elec- 
tron microscope for the inspection of small sample details with few hundred 
micrometres in size. This leads to a couple of problems, which represent chal- 
lenges, especially for small and medium-sized enterprises: On the one hand, 
elaborate laboratory characterization, which often takes place off-site, increas- 
es the duration of process development. On the other hand, microscopic imag- 
ing methods only provide information on the small analysed spot. In addition, 
the size of nanoparticles and their agglomerates differs by orders of magni- 
tude, which is a problem for many characterization methods. A further obsta- 
cle is the high investment and maintenance cost of a specialized material la- 
boratory. For this reason, many manufacturers go without microscopic 
characterization and optimize their processes instead with respect to macro- 
scopic material properties. A deeper insight into the nano- and microstructure 
is not obtained, and future potential for material development is thus limited. 
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In this work, the application of optical coherence tomography (OCT) for 
nanocomposite and nanoparticle characterization is investigated. OCT is a 
three-dimensional imaging method with microscopic resolution. Having a typ- 
ical imaging range of a few cubic millimetres, OCT complements ultrasonic, 
X-ray computed tomographic, and magnetic resonance tomographic tech- 
niques with respect to higher resolution. OCT relies on scanning the sample 
with a laser beam and on interferometrically detecting the light backscattered 
from the sample. Typical applications are to be found in medical diagnostics, 
especially in ophthalmology. 


For the characterization of nanocomposites, a multi-scale approach is pursued. 
While particles and agglomerates in the micrometre and millimetre regime are 
analysed based on images, an analysis method relying on a dedicated light 
scattering model extends the range of measurable particle sizes towards nano- 
particles. With this technique, diameter and aspect ratio of spherical and cy- 
lindrical particles in homogeneous monodisperse samples are analysed with 
nanometre accuracy. Further, the dispersion state of nanocomposites is charac- 
terized both off-line and in-line during the production process. In this context, 
the robustness of the measurement principle becomes apparent, as large sam- 
ple volumes can be analysed, and conclusions on global material parameters 
are possible. This non-destructive measurement principle proves to be suitable 
for liquid and solid materials and does not need any sample preparation. 


Industrial applications pose new challenges to OCT systems that differ from 
those of medical applications: Resolution, sensitivity and imaging speed are 
complemented by requirements with respect to robustness, small system cost, 
parallelization and small system size. Conventional free-space or fibre-optic 
systems are bulky and expensive and need exact optical alignment. However, 
temperature changes and vibrations change optical coupling properties and 
polarization states, and this is the reason why current OCT systems often have 
to be operated in protected environments and by skilled personnel. 


Photonic integrated systems in communication and measurement technologies 
offer robust, cheap and compact alternatives to traditional setups. Photolitho- 
graphic production methods allow a dense integration of a multitude of 
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systems together on a single chip. In this work, silicon photonic integrated 
OCT systems are investigated. These systems comprise interferometer and 
balanced photodetectors integrated on a common chip. System characteriza- 
tions show a sensitivity to even weak sample reflections of only —64 dB. As an 
application example, we examine biological and technical samples. 


Chapter I begins with a description of the problem by introducing nano- 
materials, their applications and the need for characterization techniques. This 
is followed by an introduction into the measurement principle of OCT. Fur- 
ther, the silicon photonics platform for system integration and miniaturization 
is presented. 


Chapter 2 expands on the fundamentals of OCT and polarization-sensitive 
OCT. An analytic expression of an OCT signal is deduced from scalar electric 
field equations. System influences on resolution, measurement range and sen- 
sitivity are investigated. Basics of optical field polarization are presented, and 
concepts for polarization-sensitive OCT systems are shown. 


Chapter 3 presents the application of multi-scale OCT for characterizing 
nanocomposite materials. Nanoparticles and particle agglomerates are ana- 
lysed, with sizes ranging from nanometres to millimetres. Large structures are 
analysed classically by OCT imaging, while nanoparticles are characterized by 
application of a dedicated optical scattering model to the measured backscatter 
data. The measurement technique is applied both in laboratory and production 
environments. A specifically developed OCT probe enables in-line measure- 
ments during the nanocomposite production process. 


Chapter 4 extends the optical scattering model from Chapter 3 towards polari- 
zation effects. Comparison of electromagnetic scattering simulations and po- 
larization-sensitive OCT measurements show that along with the size of the 
nanoparticles also their shape can be inferred. 


Chapter 5 describes the silicon photonic integration of OCT systems. In this 
context, two systems are presented, both having interferometers and balanced 
photodetectors integrated on a single chip. One system uses an internal inte- 
grated reference path, while the reference path of the other system is external 
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to the chip. This enables measurements at different distances of the sample 
from the photonic chip. A comprehensive analysis of specific signal artefacts 
reveals challenges that require special attention when designing integrated 
OCT systems. Three-dimensional images of biological and technical objects 
confirm the functionality of the systems. 


Chapter 6 summarizes the results of this work and outlines possible future de- 
velopments. 


The results obtained in this work have been published in peer-reviewed scien- 
tific journals and at international conferences. A full list of publications is 
found from page 181 onwards. 
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Achievements of the present work 


In this thesis, optical coherence tomography (OCT) systems are developed and 
applied to nanocomposite and nanoparticle analysis. This comprises fibre- 
optic swept-source OCT, both of conventional and polarization-sensitive type, 
as well as integrated silicon photonic systems. The material analysis aims at 
nanoparticles and nanocomposite materials and pursues both traditional im- 
age-based methods and scattering-model-based analysis. A summary of the 
main achievements within this work is presented in the following. 


Polarization-sensitive optical coherence tomography system 


Setup of a fibre-optic polarization-sensitive swept-source optical coherence 
tomography (OCT) system and development of a mathematical model for po- 
larization calibration. 


Light scattering model 


Development of a dedicated light scattering model, which allows the analysis 
of nanoparticles with OCT. This model enables fast and reliable characteriza- 
tion of scatterers beyond the diffraction limit [P2, P3, P5]. 


Nanocomposite dispersion state analysis 


Using OCT, nanocomposites are characterized with respect to their dispersion 
state. This comprises image-based size and shape analysis of particle agglom- 
erates inside the volume of a sample, and scattering-model-based analysis on 
the nanometre scale [P2, P6]. 


In-line process monitoring 
Application of the OCT system with a dedicated probe at a nanocomposite 
compounding extruder pilot line. The in-line measurement allows for the first 


time real-time inspection of the nanocomposite production process using OCT 
[P2]. 
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Achievements of the present work 


Size and shape determination of nanoparticles 


From measurements with the OCT system, and with the help of Mie’s scatter- 
ing theory, diameters of spherical nanoparticles in monodisperse samples are 
determined with an accuracy in the nanometre regime [P2, P3, P5]. Polariza- 
tion-sensitive OCT measurements of gold nanorods are compared to a number 
of numerical scattering simulations, and both diameter and aspect ratio of the 
nanorods are determined in good agreement with the nominal values. 


Silicon photonic optical coherence tomography 


Demonstration of the first fully integrated silicon photonic optical coherence 
tomography system. This comprises an interferometer, an optional reference 
path and a balanced photodetector. The swept-source laser is external to the 
chip. Proof-of-principle measurements of biological tissue and non-biological 
samples are shown. For robust broadband fibre-chip coupling of multiple par- 
allel fibres, polymer lenses are fabricated at the edge of the silicon chip, using 
3D laser printing [P1, P4, P7-10]. 
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1 Introduction 


Optical coherence tomography (OCT) is a measurement technique for three- 
dimensional volumetric imaging of samples. Its high measurement speed 
makes imaging possible, and resolution in the micrometre regime can replace 
biopsies in histological tissue examinations. OCT relies on low-coherence in- 
terferometry, which shares the fundamental measurement principle: Obtaining 
depth information of a scattering or reflecting object by means of a broadband 
phase-sensitive interferometric evaluation of the backscattered light. Time- 
domain OCT (TD-OCT) systems make use of the short coherence length of 
spectrally broadband light. By varying the length of the reference path, the 
measured layer inside the sample is selected by looking at the cross- 
correlation peak of the reference and the sample field. Fourier-domain sys- 
tems, as compared to TD-OCT, observe the interference signal as a function of 
optical frequency. For this purpose, either the emission wavelength of the laser 
is tuned over time, or the detection of a scattered broadband source takes place 
with a spectrometer. A general advantage of all OCT systems is the inherent 
signal-to-noise outperformance compared to direct detection techniques, 
where, in contrast, the backscattered light power is measured directly on a 
photodiode: In OCT, the backscattered light interferes with a relatively strong 
reference field, which leads to a strong beat signal. The detection principle 
makes the measurement of samples with reflectivities below -100 dB with 
respect to the incident light power possible. This allows the analysis of weakly 
scattering samples like low-index contrast biologic tissues, or strongly absorb- 
ing or scattering samples, where the signal strength decays significantly with 
increasing depth. 


First reported by Huang ef al. [1], OCT gained importance in medical applica- 
tions, like ophthalmology [2,3], dermatology [4], endoscopy [5], or in cardio- 
vascular imaging [6]. In these fields, OCT complements the conventional im- 
aging techniques, like ultrasonic imaging, magnetic resonance tomography, or 
X-ray computed tomography. OCT increasingly replaces invasive histological 
examinations, and assists minimally invasive surgery. In ophthalmology, 
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which is still the most important application field of OCT, the examination of 
the eye is accompanied by OCT images of the anterior eye [7]. Due to light 
absorption of water for wavelengths larger than 1 um inside the vitreous hu- 
mour of the eye, wavelengths around 850 nm are applied for posterior eye im- 
aging [8]. 


1.1 Nanomaterials 


The emerging multi-billion dollar market of nanomaterials is driven by a huge 
variety of applications ranging from medical diagnostics and therapy [9] to 
mechanical and civil engineering to energy storage and life sciences. This in- 
cludes the investigation of nanocomposite polymers with enhanced mechani- 
cal or electronic properties [10,11], functional coatings [12], and anode mate- 
rials for lithium-ion batteries [13]. Beyond that, nanoparticles are attractive for 
scientific examinations. Their size is in or below the range of the optical 
wavelength, and their surface-to-volume ratio is large. In contrast to bulk ma- 
terials, size and shape of the particles govern their physical properties. Thus, 
the nano-size dimensions and the distinct particle shapes become the essential 
parameters in achieving specific material characteristics. For example, gold 
nanoparticles as chemical catalysts show a significant reaction-rate- 
dependence on their size and shape [14], and also the cellular uptake of gold 
nanorods into biological tissue, when used as imaging markers, turns out to be 
shape- and size-dependent [9]. Figure 1.1 shows electron-microscopic images 
of various nanoparticles. Multi-wall carbon nanotubes (MWCNT) are depicted 
in Figure 1.1 (a) and (b) in dispersed and agglomerated form. Nanospheres 
with defined size are shown in Figure 1.1 (c). Due to their spherical shape, 
their light scattering properties are well described by Mie’s scattering theory. 
Gold nanorods in Figure 1.1 (d) are an interesting class of nanoparticles, be- 
cause their shape-dependent plasmonic resonances allow wavelength-selective 
excitation, which can be exploited for instance in tumour therapy [15]. 
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Figure 1.1: Electron-microscopic views of nanoparticles: (a, b) Multi-wall carbon nanotubes 
(MWCNT) dispersed (a) and agglomerated (b). The tubes form agglomerates because of their 
relatively strong adhesive forces and their curly shape. The spatial separation of the individual 
nanoparticles is a key challenge in nanocomposite engineering. (c) Polystyrene nanospheres as 
used as particle size standard and Mie scatterer. (d) Gold nanorods with well-defined 
aspect ratios. 


Widespread industrial application of nanomaterials is still limited by the lack 
of adequate analysis tools. In particular, achieving and maintaining a well- 
dispersed particle distribution in nanocomposites is difficult without in-line or 
at least on-line analysis methods for process control. In-line methods analyse 
the main stream of the product material directly during production, and on-line 
methods analyse a continuously taken sample stream. Nevertheless, size and 
shape characterization of various types of nanoparticles mainly relies on mi- 
croscopic imaging [16,17]. Optical or electron microscopy [18-20] need labo- 
rious, costly sample preparation and do not permit fast extraction of nanoscale 
structural information from statistically relevant sample volumes. Light scat- 
tering techniques, in contrast to imaging methods, are viable tools for deter- 
mining particle size distributions from large sample volumes. Static light 
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scattering (SLS) uses angle-resolved and/or spectrally resolved detection of 
scattered light [21,22] and needs a rather complex optical setup, when large 
scattering angles are to be taken into account. Furthermore, surface reflections 
can distort the measurement. Dynamic light scattering (DLS) methods make 
use of temporal fluctuations of interference patterns of scattered light to calcu- 
late the hydrodynamic diameter of particles within the respective solvent [23]. 
Relying on the Brownian motion of the particles in liquid media, this tech- 
nique is not applicable to solid samples. Both SLS and DLS preferentially 
provide information on particle size, while particle shape cannot be deter- 
mined [24]. Particle size distributions can also be measured using centrifugal 
sedimentation techniques, where particles with different dimension are sepa- 
rated by differences in their sedimentation speed in presence of centrifugal 
forces [25]. While this technique is applicable to a broad size range, high ex- 
perimental effort is necessary, only liquid samples can be analyzed, and no 
information on particle shape is obtained. 


None of the aforementioned techniques can meet the requirements associated 
with industrial process development and quality control, where robustness of 
the measurement method, fast analysis and the capability to characterize rep- 
resentative sample volumes without laborious sample preparation are key. In 
addition, for in-line process control, good mechanical and thermal robustness 
of the measurement system are important, as well as the possibility to integrate 
the measurement system into the processing line. 


In contrast to the established measurement techniques, optical coherence to- 
mography (OCT) provides three-dimensional imaging data from a bulk sample 
without the need for expensive sample preparation. In the following, we show 
that OCT represents an attractive tool for fast and robust nanomaterial charac- 
terization and that the technique can be applied to a wide range of particle and 
agglomerate sizes both on the micrometre and on the nanometre scale. Using a 
theoretical model of light scattering in the sample, this work demonstrates ac- 
curate determination of diameter and aspect ratio of particles in homogeneous 
monodisperse samples. Particle agglomerates with sizes of up to hundreds of 
micrometres are detected by applying image processing techniques to the OCT 
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data. Both, the scattering-model-based method and the image-based method 
can be performed in situ, without prior sample preparation, in both liquid and 
solid materials, and are applicable to laboratory investigations as well as to in- 
line process control in industrial production. 


1.2 | Measurement principle of OCT 


The three-dimensional measurement principle of OCT relies on scanning a 
thin laser beam transversely over the sample with a diffraction-limited lateral 
resolution (x, y dimension). Simultaneously, the depth-dependent amplitude 
backscattering profile r (z) of the sample is recorded. Resolution in depth 
(axial resolution, z-dimension) is achieved by measurement and evaluation of 
the broadband interference of backscattered light from the sample with light 
reflected from a reference target. Figure 1.2 shows the measurement concept 
as an example for a time-domain OCT system: Light from a light source, cov- 
ering a broad wavelength range, is split by a beam splitter into two parts, 
which travel along the sample and the reference path. Inside the sample, light 
is scattered back at different depths, according to the amplitude backscattering 
profile r,(z)=./R,(z) of the sample. The part travelling along the reference 
path is reflected by a mirror. At the beam splitter, the backscattered light from 
the sample and the backreflected light from the reference mirror are superim- 
posed, and the photodetector measures an interference pattern, depending on 
the length difference z,—z, of the reference mirror position and the position 
of the respective sample reflection. In a time-domain OCT system, the coher- 
ence length, see Appendix A.6, of the light is small, such that only sample re- 
flections from the layer with same optical distance to the beam splitter as the 
reference mirror are correlated with each other and generate an interference 
signal. The photocurrent is recorded by an analogue-to-digital converter for 
different positions z, of the reference mirror. 
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Figure 1.2: Interferometric measurement principle of time-domain optical coherence tomo- 
graphy (TD-OCT): A light source emits light (electric field Eo) over a wide optical spectrum. 
The optical power is split at the beam splitter and a first part is sent to a reference path whereas 
a second part is sent to a sample path. Depending on the sample structure, light (electric 
field E,) is scattered back from different depths of the sample, here, Zs,1, Zs,2, Zs,3, measured rela- 
tive to the beam splitter centre O with zo = 0, and interferes with the beam (electric field E,) 
from the reference mirror at position z,. The photodetector receives the electric field Eq and 
measures the interferogram obtained by scanning the position of the reference mirror, which 
contains information on the backscattering profile R, (z) of the sample. 


Generally, time-domain OCT (TD-OCT) is distinguished from Fourier- 
domain OCT (FD-OCT), which differs in the way it measures the depth- 
dependent backscattering profile: The Fourier relation of the autocorrelation 
function of a signal and its optical power spectral density is known as the 
Wiener-Khintchine theorem [26]. An analogue relation holds for the TD-OCT 
interference signal, which is the cross correlation of the fields returning from 
sample and reference path and which is linked by a Fourier relation to the op- 
tical cross power spectral density. Thus, it is possible to investigate the com- 
plex backscattering factor of the sample in frequency domain, using a spectral- 
ly resolved detection of the interference signal, rather than in time-domain 
with a moving mirror. 
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The spectrally resolved detection of the interference signal can be performed 
in two ways: A first type is referred to as spectral-domain OCT (SD-OCT) 
[27], where a broadband light source is used for illumination of sample and 
reference path, and the frequency-resolved detection of the interference signal 
is accomplished by a spectrometer, see Figure 1.3 (a). A second type is called 
swept-source OCT (SS-OCT) [28], where a narrow-band laser source is tuned 
in frequency over a broad frequency range and thereby probes the sample with 
various frequencies, see Figure 1.3 (b). 
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Figure 1.3: Spectral-domain and swept-source OCT, (a), (b), respectively. Both systems evalu- 
ate the interference signal of the light backscattered from the sample and the light of the refer- 
ence path resolved in optical frequency and do not need a scanning reference mirror. (a) SD- 
OCT operates with a broadband light source and uses a spectrometer for spectrally resolved 
detection of the interference signal, here illustrated by a diffraction grating and a detector array. 
(b) SS-OCT relies on a swept-laser source, which scans its optical frequency rapidly over time, 
as well as a photodetector for detection. 


As time-domain OCT systems can be built up with conventional opto- 
mechanical components and a broadband light source, it is regarded to be the 
simplest OCT configuration. For this purpose, most practical systems use spe- 
cially developed and fast moveable optical delay lines for scanning the refer- 
ence path length. 


However, with rising scanning speed, the mechanical effort to build robust fast 
moveable components increases and the averaging time in the detector 
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decreases, which leads to a decrease in signal-to-noise ratio (SNR). Typically, 
TD-OCT systems run at depth-scan repetition rates in the low kHz range and 
are able to detect power reflection levels of around —110 dB with respect to 
the optical power incident on the sample [29]. The low scanning speed is, e.g., 
a major problem in biomedical OCT applications, since fast and involuntary 
movements of the imaged object, e.g. the human eye, disturb 
the measurement. 


Both FD-OCT types, SD-OCT and SS-OCT, avoid mechanically movable 
parts, but require more advanced optical instrumentation as compared to TD- 
OCT, since spectrometers or swept-source lasers are necessary. FD-OCT fur- 
ther needs more data processing effort for inverse Fourier transform of the 
measurement signal, but offers higher operation speed, and higher sensitivity, 
which means a lower minimum detectable sample reflectivity level. The sensi- 
tivity advantage is inherent [30,31], since the Fourier transform reduces the 
noise level by spectral filtering. These systems are able to reach shot-noise 
limited performance and can detect around —120 dB sample reflectivity. Fur- 
ther, without need for mechanically scanned delay lines, much higher scan- 
ning speeds become possible: Choi et al. showed a 10 MHz depth-scan repeti- 
tion rate with an SD-OCT system and Wieser, Klein et al. demonstrated SS- 
OCT systems operating at repetition rates between 3.2 MHz and 
20 MHz [32-34]. 


1.3 Silicon photonics 


While conventional OCT systems are predominantly built from discrete com- 
ponents, photonic integration facilitates compact systems that stand out due to 
robustness and that lend themselves to cost-efficient mass production. Espe- 
cially in endoscopic applications as well as in industrial surface and distance 
metrology, integrated systems can comply with the demands in size, mechani- 
cal robustness, temperature stability and cost. 


Silicon as the basis for photonic integrated circuits (PIC) gained attention as it 
proves to be a material platform, which allows the reuse of CMOS fabrication 
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infrastructure. A plurality of foundry services has emerged within the last dec- 
ade [35-37], primarily driven by applications in optical communications. The 
foundries offer process design kits (PDK) with libraries of standardized silicon 
photonic devices such as passive components, electro-optic modulators, and 
germanium-based photodetectors. 


Other material platforms, e.g., based on III-V material systems, inherently 
contain active materials for laser-source co-integration. Nevertheless, the low- 
cost and high-yield process flow in silicon photonic fabrication is an important 
advantage when aiming at large-scale photonic integrated systems. 


Besides compatibility to CMOS fabrication lines, large-scale and inexpensive 
system-on-chip co-integration of photonics and electronics becomes possi- 
ble [38]. Figure 1.4 shows a schematic cross-sectional drawing of a silicon-on- 
insulator (SOI) chip with passive photonic and active electro- 
optic components. 
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Figure 1.4: Schematic cross section of a silicon photonic chip, including passive optical com- 
ponents and active electro-optical components. The silicon-on-insulator (SOI) wafer with func- 
tional silicon (Si), germanium (Ge), metal (M1, M2) layers and p-, n-doped silicon allows, 
among others, the co-integration of waveguides, grating couplers, photodetectors 
and modulators. 
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This chapter provides the theoretical background of optical coherence tomo- 
graphy (OCT) including a theoretical treatment on OCT signal calculation. 
Building on this, relations for achievable resolution and sensitivity are formu- 
lated, which arise from technical limitations in wavelength range, lens aper- 
ture, scanning speed and noise of available components. A section on sensitiv- 
ity analysis compiles influences of different noise sources and suggests 
operating parameters for optimum OCT performance. 


In a section comprising fundamentals on optical polarization, polarization- 
sensitive OCT (PS-OCT) is introduced, together with PS-OCT variations and 
applications thereof. Further, the individual PS-OCT setup built in the course 
of this work is presented, which is used for the polarization-resolved meas- 
urements shown in Chapter 4. 


2.1 Optical coherence tomography 


In this chapter, the signal obtained from an OCT measurement is derived ana- 
lytically. The emphasis is on the calculation of the depth-dependent backscat- 
tering signal, which is the central element of OCT theory. Information on lat- 
eral resolution can be found in Section 2.1.2. 


The signal calculation can either be performed in the time-domain [39,40] or 
in the frequency-domain [41]. In the following, the signals and fields are de- 
rived in time-domain, starting from the electric-field impulse responses of the 
sample and of the reference path, and convolving them with the time- 
dependent source field. Superposition of the fields from sample and reference 
path leads to interference. The resulting intensity is detected by a photodetec- 
tor and is described with the help of coherence functions. These allow the 
analysis of the depth resolution and the measurement range. The results of this 
calculation are applied to TD-OCT and, using Fourier-transform, to the FD- 
OCT types SS-OCT and SD-OCT. 
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Complementing these considerations related to depth resolution, Section 2.1.2 
is dedicated to lateral resolution and focal measurement range, which both re- 
late to the optical beam properties. 


In addition, the effect of sampling is addressed in Section 2.1.3, together with 
properties of the discrete Fourier transform, since the OCT signal is processed 
in the digital domain. 


Finally, in Section 2.1.4 the sensitivity of an OCT measurement, describing 
the minimum measurable sample reflectivity, is derived from the noise contri- 
butions of the individual system components. The concept of balanced detec- 
tion for sensitivity improvement is introduced, and the dependence of the sen- 
sitivity on reference and sample arm power is derived. Optimum operating 
points for different measurement tasks are suggested. 


2.1.1 OCT signal calculation 


The light source in OCT, see Figure 1.2 and Figure 1.3, can be described by 
the time-dependent electric field of the emitted light, which is considered to be 
a plane wave and which can be represented by the complex scalar electric field 


E(t) = Ap (t)i. (2.1) 


In this relation, A, (t) is the complex slowly-varying random amplitude, 7 the 
time and @&,=2r.f, the optical angular carrier frequency of the source with 
carrier frequency f,. In the case of a SS-OCT system, the optical angular fre- 
quency is additionally varied over time. The light is split up by the beam split- 
ter at position O with z,=0 on the distance axis z, which describes the dis- 
tance of a reflection with respect to the beam splitter location. In the 
following, this axis is also referred to as the depth axis. In both sample and 
reference path, light is scattered and reflected back, which can be described by 
the corresponding field reflectivity profiles r (z)=4/R, (z); 7.(z)=JR,(<)- 
The reference-path reflectivity profile is given by a ö-function, see Appen- 
dix A.4 for a mathematical definition of the ö-function, with field reflectivity 
mg at the depth of the mirror position z,, 
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r.(z) =m (z-z). (2.2) 


The field impulse responses of sample and reference path describe their tem- 
poral back reflection behaviour and are derived from the field reflectivity pro- 
files using the respective speed of light in sample and reference path, 
V; = Cong» Vy =Con,, With the vacuum speed of light cy and the effective refrac- 
tive indices n,,n,. For simplicity, we assume the effective refractive indices 
n, and n, to be constant along z in the following. The depth positions 
z = (vr eae correspond to round trip times 7, „, where the factor of 1/2 
accounts for forward and backward propagation in sample and reference path, 
respectively. The field impulse responses of sample and reference path thus 
are given by 


(2.3) 


The factor of dz/dt = v, , / 2 accounts for the conversion of a reflectivity per 
depth to a reflectivity per time. With Eq. (2.2), the return time 7, correspond- 
ing to the reference mirror in the reference path, and z, = (v, / 2)r, , the impulse 
response of the reference path results in 


h(t.) =r ol“ -))e 


(2.4) 


SR olt -7.). 


The fields returning from sample and reference path can be calculated by a 
convolution of the incident fields with the impulse responses. The electric 
fields returning from sample and reference path are denoted as Æ, and E,, 
respectively, and considered in front of the detector in the detection path, see 
Figure 1.2. Each of the fields is once reflected at the beam splitter with a field 
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reflection factor of j/ ay and once transmitted through the beam splitter with 
a field transmission factor of 1/ 2 . This leads to a total factor of j/2 for both 
of the fields, 


E, (t) = 5 Eo(t) * hg (t) 
(2.5) 
E 


j 
E(t) = 5 Bolt) *h (t) 
where * denotes a convolution in the time domain, see Appendix A.3. 
Eqs. (2.5) can be expressed as 


+ +00 


E,(t)=2 h(t JA (t-t,) "ar, 


(2.6) 


E,(t)=2 | h,(t,) Ag(t—t,) 1a, 


= İn A(t- r, Jal), 


In these equations, the effect of sample and reference on the electric field be- 
comes visible: The amplitude is scattered back with factor of m or 
r ((v,/2)t,) along with a phase shift —@pr, or —@pt,. The variables t,t, de- 
note the respective return times from sample path and reference path. 


Both fields from sample and reference path superimpose on the photodetector, 
which generates a photocurrent proportional to the averaged squared absolute 
value of the complex electric field 


|El =|E, + Ef 
S E E 27) 
=|E,| +|E,| +2Re| EE}. 


This average is a time average, and we assume that the fields are described by 
ergodic processes, see Appendix A.6. The square leads to mixing terms, of 
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which the sample-reference mixing term contains the depth-dependent scatter- 
ing information of the sample. In the following, all terms are evaluated sepa- 
rately, starting with the sample-reference mixing term, which is the actual 
OCT signal. 


With Eq. (2.6) the sample-reference mixing term is given by 


2Re{ EE; | 


= si Re Í h, (1; A(t -r Jel) A (t z r, Je Holt ay (2.8) 


= T rel [n (1; )A (t —t,) Ag (t = “ene 


—00 


The integrand contains a correlation of two fields with temporal distance 
T, —t, and can thus by expressed in terms of a complex coherence function, 


see Appendix A.6, 
Tl ~.)=4o(t—%) 45 (t-7,) FO) , 29) 
resulting in 
1 +0 
2Ref EE |} = rie f WC )Re{EL (7 =4)} di, (2.10) 


This term can be rephrased using the real-valued coherence function 


Ki (t, -1,)=ŻRefE, (r, -=t )} (2.11) 


leading to 


+00 
2Re{ EE, | = Fig J I(t.) Ku (T -1)di. (2.12) 
The real-valued coherence function can be expressed by the normalized real- 
valued coherence function k,(r,-1)=Kı(r.-t,)/K,(0), with 
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Kı (0)=|4 (£)? . In the case of quasi-monochromatic sources (fy > 6f ), 
where the instantaneous source bandwidth Sf is significantly smaller than the 
carrier frequency, which is typically valid for OCT light sources, the coher- 
ence function can be split up in a baseband coherence function Kk, (c) and a 
harmonic term, 


Ki (t, -t,) = Kz, (0), (2, -6) 


(2.13) 
=K, OLS (T -t, )cos(@, (7, -t )). 
Hence, the sample-reference interference term is expressed as 
2Ref EE, | 
+00 
=K, (0) f h(t, )KL (7, -t,) cos(@ (t-t )) dt, (2.14) 


X 


— +0 
= EMO f n( Sa hala —t,) cos (a(t; ~t,)) dt,. 
This equation shows that the signal is proportional to the square of the source 
amplitude Ay (t) and to the reference mirror reflectivity rg. Further, the signal 
is given by an integration of all reflectivities r, (t; v, / 2) in the optical sample 
path. The contribution strength of each reflection is governed by the baseband 
coherence function K (t, -t,) and thus depends on the difference in path- 
length distance with respect to the reference path. The cosine-term leads to a 
harmonic oscillation of the signal when tuning the carrier frequency or chang- 
ing the reference mirror position. 


The OCT signal, Eq. (2.7), also contains a part that is proportional to the opti- 
cal power returning from the reference arm. With Eq. (2.6) this part is 
given by 


i 


. 5 2 
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In contrast to the sample-reference interference, the mixing of reference with 
reference fields results in a constant term, since only reflections from a single 
depth occur. 


Moreover, the OCT signal in Eq. (2.7), contains a part that is proportional to 
the optical power returning from the sample arm. Using Eq. (2.6), the mixing 
of fields from all sample reflections is calculated to 


. +00 


E*|=|2 Í h,(t,) Ap lt -i)a 


Z f n (e)a (e-e) ag (2.16) 
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using the integration variables t ,t; for discrimination of the two integrals. 
The harmonic sample-sample interference term depends on the path length, 
and thus on the temporal differences t,—f{ between the individual sample 
backscatterers. In order to calculate the interference signal, the integral is first 
written down twice, then in the second part, the integration variables are ex- 
changed, which results in an expression of the complex conjugate (CC) of the 
first part. Hence, an expression of the real value results, 
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Using the definition of the real-valued coherence function in the baseband, 
Eqs. (2.9), (2.11) and (2.13), the sample-sample interference term is 
expressed as 


—_— +00 00 
IEPES J a (a)i(K)Re {Eol -janar 
BEN (2.18) 
= Slash Í Í h(t, hg (£) my (th -t )cos(@ (f -1,))dr,dr.. 
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In comparison to the sample-reference interference, Eq. (2.14), the sample- 
sample interference depends on the mutual return time differences from re- 
flexions inside the sample. Due to sample reflectivities that are typically 
weaker than the reference reflectivity, the sample-sample mixing term is also 
weak, compared to the sample-reference interference. 


The calculation leads to an entire OCT signal, which is dependent on the opti- 
cal return time differences in both reference and sample path. By inserting 


2.1 Optical coherence tomography 


Eqs. (2.14), (2.15) and (2.18) into Eq. (2.7), the squared electric field on the 
photodetector amounts to 


esa) 2 
Efa Sr 
const. term 
ee 
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sample-reference interference: OCT signal 


In the following, the sample-reference interference signal i, = 2Re}E,E, + is in 
the focus, which represents the actual OCT measurement signal, compare also 
Eq. (2.14). Constant factors like the area and the responsivity of the photode- 
tector are omitted. We describe this signal as an effective photocurrent i., 
since it is proportional to a detected photocurrent i,, where the subscript de- 
notes correlation of sample and reference path. The time integral in Eq. (2.14) 
is now evaluated for the specific cases of first a TD-OCT system and second 
an FD-OCT system. 


In case of a TD-OCT system, the light source emits broadband light, which 
means the complex slowly-varying random amplitude A(t) is broadband 
noise, and thus the coherence time T, of the light is small. The coherence time 
is the time span, in which the normalized baseband coherence function is ap- 
proximately one, k(t, -1,) xl, if Ir, -t <r,. For larger time differences 
Ir, A >T,, the coherence function is approximately zero. The measurement 
signal in Eq. (2.14) contains the once (a sample impulse response with 
the coherence function, /,( )*[ x. (¢, )cos(@pr, ) )]- In an idealized case, 
A(t ) represents noise of ae 1/c,. With the approximation of t, > 0, 
the convolution approaches the term h, (ers Accordingly, the measurement 
signal in Eq. (2.14) becomes 
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v, (2.20) 
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This shows that the signal current is in direct proportion to the impulse re- 
sponse and to the reflectivity profile of the sample. In the realistic case of a 
limited bandwidth of the source, the convolution of the sample impulse re- 
sponse with the coherence function leads to a limited depth measurement reso- 
lution of z =v;r.. 


In the case of an FD-OCT system, both in the sense of an SD-OCT or in the 
sense of an SS-OCT, narrow-band light is applied, either by using a tunable 
narrowband light source (SS-OCT) or by spectrally filtering the detected light 
(SD-OCT). Thus, with a narrow bandwidth, a long coherence length allows 
the interference of backscattered sample light with light travelling along the 
reference path for long path-length differences. As a consequence, the refer- 
ence mirror does not have to be moved any more, and interference from vari- 
ous depths in the sample can be measured simultaneously. In an ideal case, the 
coherence time T, of the detected light is infinite and the baseband coherence 
function Ko(t) 1 equals one for all time differences T. The measurement 
signal in Eq. (2.14) then results in 


+00 


i! pp (©) = Ki (0) re Í h,(t,) cos(@(z, -1,)) dt. (2.21) 


oo 


In FD-OCT, the frequency-dependent photocurrent i pp (a) is in direct pro- 
portion to the Fourier transform of the temporal field impulse response of the 
sample. The following calculation shows this relation. 


The cosine in Eq. (2.21) can be expressed by exponential functions 
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which represent Fourier transforms (FT, see Appendix A.1) of the sample im- 
pulse response, 
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+00 ; . +00 l (2.23) 
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The complex field transfer function h, (f) of the sample is the Fourier trans- 
form of the temporal field impulse response 
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h, (f)= J h (i)e Par, . (2.24) 
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The Fourier integrals in Eq. (2.23) can then be written as 


' 1 7 -j27 for, 7 j27 for 
kro (fo) => Ki (0) re [As(-fo)e ehe +h, (fo)e? fo | (2.25) 


Equation (2.25) expresses the measurement signal of a FD-OCT system and 
shows that the measurement signal is proportional to the complex field trans- 
fer function of the sample. 


In contrast to the general case of a sample with an arbitrary complex impulse 
response, the following section illustrates the case of a plane mirror as a test 
sample in an FD-OCT measurement. The test sample leads to a single reflec- 
tion with a return time t, ı and a real-valued field amplitude reflectivity of 7, ,, 
leading to an impulse response of he s(t.) =" 5(t, -1.1)- Accordingly, the 
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-j2r ft, ı 


complex field transfer function becomes h sl (f)= re and the meas- 


urement signal of such a sample is 
„i 1 ~i2afo(t ts, J27 fo(tr-ts, 
isro (Jo) = 5 Ki (O)re [re le) + He le) 


=K (0) Rh1 cos(27 fo (z, - tea) 


This equation describes a real-valued signal, which is dependent on the optical 


(2.26) 


center frequency fo of the detected spectral slice. The FD-OCT measurement is 
performed for a range of optical frequencies, ( fn — Af /2) < fo < (fm + Af/2), 
with centre frequency f, and measurement frequency span Af. For the de- 


termination of the time- or equivalent depth-dependent signal, an inverse Fou- 
rier transform can be applied. Because the input spectrum is real-valued, this 
transformation results in a complex symmetric time-domain signal with com- 
ponents at positive and complex-conjugate components at negative times, with 
respect to the reference arm return time 7,. A sample reflection with a positive 
return time t,,;—7z,>0, thus appears also at the negative time difference 
he 
negative times or depths is often referred to as depth degeneracy [42]. 


and vice versa. The inability to discriminate between positive and 


For SS-OCT, the carrier frequency fo (t) is subject to a temporal variation, 
and the signal associated with a single reflection within the sample, Eq. (2.26), 
is given by an oscillation along the carrier frequency axis, 
cos(27 fo(t)(r, -t,.)). The return time difference 7,— t, between a sample 
reflection and the reference path return time causes an intermediate frequency, 
whose value allows the determination of the time difference. If the optical car- 
rier frequency fọ is increased over time, see Figure 2.1, with a linear tuning 
slope Sw =Af/T,,, and frequency ramp duration T,,,, starting at t= ftg, the 
time-dependent carrier frequency is given by fo (t) =( fmn -Af /2) + Sswlt to). 
This dependence allows to calculate the electrical signal frequencies fig) OC- 
curring in the SS-OCT measurement signal 
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The measured signal has a phase offset @. Note that in practical tunable la- 
sers, the different frequency points might be subject to different offset phases. 
This signal model, however, assumes a chirped laser source with a continuous 
phase. Any initial phase offset would be present in both sample and reference 
arm and would thus cancel out in the measurement. The frequency of the 
measured signal is 


Soig,1 = (z, ae = (z tea Af [Tose : (2.28) 


The temporal dependence of the carrier frequencies fo, returning from refer- 
ence and sample path are depicted schematically in Figure 2.1. 


The frequency ramps are assumed to be linear in time with a frequency span 
of Af and a ramp time Tsw, giving a slope Ssw = Af/ Tsw. Depending on the re- 
turn time difference 7. — ts, of reference and sample arm, the intermediate fre- 
quency of the OCT signal by heterodyne down-conversion amounts to 
Jigi =(= ts.1) Ssw- 


In order to estimate the signal frequencies in a real system, an exemplary case 
is investigated with a measurement frequency span of 
Af, = (AA, / A )co =17.5 THz, which corresponds to a wavelength span of 
AA, =100 nm at a centre wavelength of A, =1310 nm. Here, the frequency 


ramp time is assumed to be 7,,,;=>500us and the return time difference of 


sw, 1 7 
reference and sample path is given by 7, -1,, =2z /c,. Sample and reference 
path velocities are assumed to match the vacuum speed of light, and the depth 
difference of both paths is z,=1mm. In this case, the measured signal fre- 


quency is expected to be frig = (221/0 (Af Tw) = 233 kHz. Thus, 
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photodiode bandwidths and analogue-to-digital converter sampling speed 
should be adapted to the expected signal frequency. 


teeth fall = te) 


Jm Tei. © E oe figi = (T g ts,1) Ssw 
-° 


Figure 2.1: Frequency ramp used for an SS-OCT scan. The carrier frequency fo is increased 
linearly over time f. One scan of frequency range Af centred around a frequency fm takes the 
time Tsw, thus the frequency slope is ssw = Af/ Tsw. The carrier frequency fo(t — 7,) (solid) of the 
light returning from the reference path with round-trip time 7 and the carrier frequency 
Jolt- tsı) (dashed) of the light returning from the sample path with round-trip time tsı have a 
frequency difference of fsig.1, which gives the intermediate frequency of the mixing product. 


The limited frequency span has an influence on the depth resolution of the 
measurement. In practice, the frequency span Af is limited due to either the 
source tuning range (SS-OCT) or the emission bandwidth of the superlumi- 
nescent diode (SD-OCT). In order to account for this, the measurement signal 
in Eq. (2.25) is considered as a spectrum, multiplied by a windowing function 
in frequency, w( fy) = tect Af ( fo- TAP here given by a frequency-shifted 
rectangular function, see Appendix A.1 for the definition and Figure 2.2 for a 
visualization. Accordingly, the measurement signal is denoted by the follow- 
ing equation 


7 l i, j27 for, D -j27 for, 
itro (fo) =5 Ki (0)RW (Fo) Bs (foe? + hs (fo)e P | 2.29) 
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Im Jo 


Figure 2.2: Window function W(fo) describing the bandwidth-limited scanning range of the 
FD-OCT signal, centred around fm. In this case, the window is given by a rectangular function 
of width Af. 


In order to obtain the temporal response of the sample from the measurement 
signal, an inverse Fourier transform with respect to fọ is applied to Eq. (2.29), 


i, l T p J27 for, „J2r 
ito (t,)= 5K. (Om | [WeH)E (A) reg, 
x (2.30) 
+ | W( fo) bg (fo) reg, |. 


The windowing function and the measurement signal can be transformed sepa- 
rately. The transformation of the windowing function 


W(t,)= | W(fo)e?* at, (2.31) 


is in this case the transformation of a frequency-shifted rectangular function 
and is expressed by the following modulated sinc-function 


+0 


W(t,)= Í rect ar (fo - fm er sd, 


oo 


= Af sinc(zt,Af Joints 


(2.32) 


The time-domain width of the sinc is reciprocal to the measurement bandwidth 
Af. The exponential phase modulation can be omitted, as the representation of 
the measurement as an A-scan usually contains only the absolute value. 
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The former multiplication of measurement signal and windowing function, see 
Eq. (2.30), is now expressed by a convolution of their individual transfor- 
mation products 


7 1 17 P j27 fot, n j2 7 
ice (t)= > Ki (0) re W(t,)* AOG Sot: ei? Sols af. 
u (2.33) 
+W (1) * Í h,(—fy)e P eP Adf ; 


The inverse transform of h,( A) is identical to a time-shift of the 
sample impulse response h,(t,) by the reference path return time 7, see 
Eq. (A.4). With the auxiliary variables ģ =t, +7, and ¢, =t, —t, the integrals 
in Eq. (2.33) are 


a (2.34) 
Í h, (- fo)e df = h, ( Ey) 7 h, (T, t) 


and the entire Eq. (2.33) results in 


& 1 = en 
opp (t,) = 5K (0) rig W(t.) *A, (T, +1,)+W(t,) *A,(z,-1,) }- (2.35) 


Expressed in sample reflectivities, this reads 
(2.36) 


The inverse transform of the measurement signal corresponds to the sample 
reflectivity profile mirrored at z,, but is mainly influenced by three effects: 
Depth degeneracy, window artefacts and coherence decay. These are ex- 
plained in detail in the following section. 
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The measured profile r, (z + Z, ), r (z u z) occurs twice, which is referred to 
as “depth-degeneracy”. This effect is related to the real-valued input spectrum, 
Eq. (2.21), leading to a symmetric result of the inverse Fourier transform. Pos- 
itive and negative path-length differences cannot be discriminated and super- 
impose. Further, since each singular reflection is broadened by the convolu- 
tion with the window response W (t), this causes a limited resolution 
ôt, =1/Af in sample response time and accordingly a limited depth resolution 


öz=v,/(2Af). 


Up to now, the coherence time +, of the source was assumed to be infinite. 
Now the case is regarded, where the baseband coherence function x; ( C= t) ; 
see Eq. (2.19), is limited to a finite width z,.. This time is typically large 
enough to allow coherence lengths in the millimetre range, but significantly 
limits the imaging depth in the sample. With this influence, Eq. (2.35) can be 
written as 


iro (t)=5K (0) LOLA (eh le +) 
+W(1,)* x, (7, ~t, Jh (T, -t,)|. 


The impulse response of the sample h,(t,) is multiplied by K (t, -t,) and 


(2.37) 


the signal intensity decreases, once the path-length difference of sample and 
reference path increases. The measurable range of sample response times is 
limited by the coherence length Ar, =7,, leading to a maximum measurement 
depth of Az =7,v,/2. Figure 2.3 visualizes the different effects on a FD-OCT 
measurement signal. The amplitudes of all signal components are plotted ver- 
sus response time f,. 


The original sample response is assumed to consist of three infinitely sharp 
responses at times f1» 1,2, t3 (black). The reference path return time is T, 
(dotted vertical line). The measured sample response, see Eq. (2.37) (absolute 
value in blue), is subject to mirroring at the reference path return time, to a 
decay with the coherence function K (7, =t) (red), and to broadening by 
convolution with the time response of the spectral window W(t) (absolute 
value in green), here assumed to be of a sinc-shape, i.e. by a rectangular 
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spectrum, Figure 2.2. Note that an actual OCT signal would contain parts from 
sample-sample and reference-reference interference as well, which are both 
not included in this investigation. For a mitigation of the undesired signal 
components emerging from interference of fields from the same path, most 
practical measurement setups include balanced detection, where these signal 
components are suppressed, see Section 2.1.4. Here, only the signal originat- 
itr (ts) 


ing from sample-reference interference is shown. 


Amplitude 


i Original sample response 
aai, r 


ancol | Wren 


Measured sample response 


irn (f) 


0 ts 0 A fsa ts 153 


Figure 2.3: (a) Visualization of the absolute value of the time response of the spectral window 
Wi) , here assumed of sinc-shape. (b) Visualization of an original (black) and measured by 
FD-OCT (blue) sample response, h, (i) I FD (1) 
ized as three sharp peaks. The measured profile is subject to three effects: mirroring of the orig- 


, respectively. The original profile is ideal- 


inal profile at the reference return time 7, , decay of the coherence function KL ( = ), (red), 
and broadening of each peak by the time response of the spectral window W (t) i 


2.1.2 Lateral resolution and focal measurement range 


In OCT, three-dimensional images are obtained either by laterally scanning a 
laser beam over the sample or by moving the sample position laterally, both in 
addition to the depth-resolved backscatter measurement. Typically, the beams 
are considered to be Gaussian beams. In a conventional fibre-based system, 
standard single-mode fibres (SMF) guide light towards and away from the 
OCT scan head. The scan head includes scanning mirrors and a scan lens. The 
fundamental mode in these step-index fibres is the linearly polarized LPo; fi- 
bre-mode [43]. Its radial field distribution can be approximated by a Gaussian 
distribution. This Gaussian shape is conserved after transmission through the 
collimator and the scan lens objective [44] and results in a focussed Gaussian 
free-space beam in the sample, see Figure 2.4. 
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In the following, the lateral resolution and the axial measurement range of an 
OCT system are derived from the properties of a Gaussian beam, see Appen- 
dix A.5. The complex field amplitude of the dominant transversal field com- 
ponent of a Gaussian beam, as depicted with its radial amplitude and its axial 
profile in Figure 2.5, is given by 


2 
Tee Tol 2.38 
Ag (r.z) Ago wa) E ( ) 
Accordingly, the intensity distribution of a Gaussian beam is 
2 
Ig(r,2) ~ [Ac (2) 
2 
ADNE 2.39 
A co re 
> 2 (z) 


with the radial distance r form the beam axis and the axial position z. 


Divergent beam; Collimated beam | Focussed beam 


Single-mode fibre 


Sample 


N 
Collimator Sie Focal 
lens plane 


Figure 2.4: Schematical drawing of the Gaussian beam, which is emitted by the single-mode 
fibre. The beam is then collimated by the collimator lens and finally focussed into the sample 
by the scan lens. The focal plane of the scan lens defines the plane of best lateral resolution. In 
the drawing, the sample refractive index is assumed to match the surrounding refractive index. 
A scan mirror (not shown) can be placed in the collimated beam section. 


The maximum intensity in the focal spot is given by Jgo and the radius of the 
beam waist is wo. The beam profile w(z) widens with axial distance from the 
position z of the beam waist according to 
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(2.40) 


see Eq. (A.23) in Appendix A.5. For an analysis of the lateral resolution of the 
OCT system, we first calculate the focal response of the system, which is de- 
fined by both illumination and reception of the scattered power by the same 
optical system. This can be explained by a point scatterer as measurement tar- 
get, which is exposed to the incident light, with an intensity distribution given 
by Eq. (2.39). Additionally, the efficiency for coupling scattered light back 
into the reception optics is governed by the same relation and reduces with 
radial and axial distance from the focal spot. Thus, the focal power response 
function 77 (r, z) is given by the squared normalized intensity 


Ialr2)_ wô Aa l (2.41) 
Io w*(z) 


From the focal response function, two important system parameters can be 


n (17,2) = 


derived: the lateral resolution and the axial measurement range. We define the 
lateral resolution in the focus at z = zę as the span between the points, where 
the power response is 1/e?. The lateral resolution thus is given by 


dr =V2 wo. (2.42) 


For the exemplary fibre-based OCT system described in this work, the beam 
waist is wọ =18.8 um, and the lateral resolution at the beam waist is 
or =26.6um. 


The focal measurement range Az; is determined from the decay of the focal 
power response function 77; (r, z) along the beam axis (r = 0) to 1/4. The span 
between both points with 7; (0,z) =1/4 is given by the confocal parameter b 
of the underlying Gaussian beam 
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The focal measurement range in the system used in this work is Az; =1.7 mm. 
The wider the beam waist and thus the worse the lateral resolution is, the 
longer is the focal measurement range. For an optimum OCT operation, it is 
useful to design the system such that the decay caused by coherence effects 
and the decay due to the axial focal response function are in a comparable 
range and hence do not unnecessarily limit each other. By choosing the coher- 
ence point in depth apart from the focus at z;, both functions can even counter- 
act and form a rather flat depth-response of the system. 


exp[-(r/w(z))2] 


Figure 2.5: Gaussian beam illuminating the sample. Left: Normalized radial amplitude profile. 


Right: Axial profile of the beam with the beam waist of radius wo at the position z and the 
depth-evolution of the beam radius w(z). 


Some OCT systems use Bessel beams instead of Gaussian beams in order to 
improve the focal measurement range in scattering media. These beams are 
characterized by an energy flux from the outside lobes to the inner main lobe 
and are thus said to have self-healing properties [45] in scattering media. Bes- 
sel beams can be generated by sending a Gaussian beam through an axicon. 


2.1.3 Sampling and DFT for acquisition of FD-OCT signals 


For processing, display and storage, the OCT voltage signal is digitized, either 
by sampling and quantization of the photodetector output signal using an ana- 
logue-to-digital converter (ADC) in SS-OCT, or by digitizing the spectrometer 


31 


2 Fundamentals 


output, e.g., with a charge-coupled device (CCD) followed by an ADC, in SD- 
OCT. Sampling the signal in optical frequency steps fo with a limited number 
of samples M and a step size of Of, covers a limited range of the optical 
spectrum Afy=Mof,, see Figure 2.6 (a). Since positive and negative fre- 
quencies cannot be discriminated, only the positive part of the spectrum is 
sampled. The signal is processed by an (inverse) Fourier transform, see 
Eq. (2.30), which is performed digitally after sampling, using a discrete (in- 
verse) Fourier transform (IDFT) algorithm, see Appendix A.2. The transfor- 
mation with an (DDFT generally links discretization in one domain with peri- 
odicity in the other domain. Specifically, the limited spectral observation 
range Af, leads to discrete data points separated by the resolution in return 
time ofy = 1/ (2Af,), and the step size ôf} in frequency-space leads to a lim- 
ited unambiguity range Af, =1/(26 fa) see Figure 2.6 (b). The number of 
points M stays the same before and after the (DDFT, however, after an 
(DDFT, the data points are distributed over positive and negative sample path 
return times f,. The number of relevant data points is reduced to M/2, since 
the complex amplitudes at positive and negative times are complex conjugate 
with respect to each other, thus containing redundant information. This fact 
corresponds to the depth-degeneracy as artefact in OCT signal reconstruction 
introduced in Chapter 2.1.1. 


Figure 2.6: Illustration of discrete-Fourier-transform-(DFT-)related quantities in both frequen- 
cy domain and time domain. (a) Sampling points 0, 1, ... , M - 1 in the frequency-domain, with 
total sampling range Af}, step size 6 f,, and total number of sampling points M. (b) Data 
points 0,1,..., M-1 after inverse DFT (IDFT) with one-sided unambiguity range of 
Ar, = 1/(26 fa ), discretization Sty = 1/( 2Af,), and the number of points M. 
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The system design is chosen reasonably if the sampling-related temporal reso- 
lution ör, and the coherence-related temporal resolution dt, are in the same 
range. Furthermore, the sampling related unambiguity range Aty should not 
be significantly smaller than the coherence-related temporal range Af, . Oth- 
erwise, the measurement signal is superimposed by periodic backscatter from 
sample return times exceeding Af, . 


2.1.4 Sensitivity considerations and noise analysis 


The coherent mixing process of reflected sample light with reference-arm light 
is key to the high sensitivity of OCT, since the measured signal power in the 
electrical domain is proportional to both, the sample arm power and the typi- 
cally higher reference arm power, see Eq. (2.19). In order to achieve systems 
with high sensitivity, one would apply strong reference arm powers to lift the 
signal power above the noise floor of the electronics. However, with rising 
power, light source related noise contributions also increase, like the relative 
intensity noise from the laser, while the signal power might run into a detec- 
tion limit, given by, e.g., the maximum input voltage of the analogue-to-digital 
converter. Hence, a comprehensive noise analysis becomes important to 
achieve the optimum operating point. This section is focussed on the situation 
in fibre-based SS-OCT systems. 


The following considerations on noise and sensitivity in an SS-OCT system 
rely on the configuration as shown in Figure 2.7, which is exemplary for the 
SS-OCT systems used in this work. The wavelength-swept source (SS, Micron 
Optics s3) with a wavelength range of 1260 — 1370 nm and relative intensity 
noise RIN4p =—100 dB Hz! at 2 MHz is connected to a fibre-based interfer- 
ometer, consisting of two 3-dB couplers (CPL1,2). The first coupler CPL1 
splits the light power from the source among the sample and the reference 
path. The sample path (SP) includes a free-space propagation section between 
fibre collimator (FC) and sample, and the reference path (RP) is completely 
fibre-based. Optical path lengths of RP and SP are considered to be matched. 
The second 3-dB coupler CPL2 acts as combiner and feeds the interfering 
light to both inputs of the balanced detector (BD, Thorlabs PDB430C), which 
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features a responsivity of o=1 A/W, an optical noise-equivalent power 
NEP =7 pw/ Hz , a common-mode rejection ratio (CMRR) of ¢ =—25 dB, 
and a transimpedance gain of Z, =5000 V/A for a 50 Q load. Note that, for 
this device, the transimpedance gain is Z, =1x10* V/A for a load with a 
high input impedance. This number reduces by a factor of 2 for a 50 Q -load 
since the output resistance of the transimpedance amplifier also amounts to 
500. The NEP is defined by an optical input power 


Fopt,NEP = NEP\ Afa > (2.44) 


which leads to a signal-to-noise ratio at the detector output of one, assuming 
an electrical measurement bandwidth Af,,. 


Figure 2.7: Swept-source OCT setup considered for sensitivity calculations. SS: swept-source 
laser (Micron Optics s3, 1260-1370 nm, RIN =1x10"'° 1/Hz), CPL1,2: 3-dB-coupler, SP: 
sample path, RP: reference path, FC: fibre collimator, BD: balanced detector (Thorlabs 
PDB430C, responsivity of o=1 A/W, noise of NEP = 7pw/VHz , CMRR of ¢ =-25 dB 
and transimpedance gain of Z, =5000 V/A for a 50 Q load), ADC: analogue-to-digital con- 
verter (Alazar ATS660, 16 bit, voltage step size Öl, =392 uV). 


The difference signal of the BD is digitized by an analogue-to-digital convert- 
er (ADC, Alazar ATS660) with 16 bit resolution and a voltage step size of 
Oug = 392 uV in the low-noise operating mode, corresponding to an effective 
number of 12 bits. The polarizations of the light from the RP and SP are as- 
sumed to match and are hence not regarded explicitly. Deflection and 
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focussing optics in the sample path are neglected. The sample power reflectiv- 
ity R, (z) = r (z) is the squared sample field reflectivity and is defined by the 
backscattering into the fibre-optic system, thus including coupling losses into 
the optical system. 


Using this setup configuration, the following section introduces the concept of 
balanced detection and evaluates the impact of the different noise contribu- 
tions on the overall SNR. As will be shown in this section, highest sensitivity, 
allowing measurements of lowest sample reflections, is only achieved with 
this system using rather low reference path powers, when RIN contributions 
from the light of the reference path are small. This situation can be different 
for sources with low RIN or for perfect balancing over all wavelengths. Real 
systems, however, make this comprehensive noise consideration necessary. In 
contrast, reaching the maximum SNR and thus exploiting the dynamic range 
of the system is shown to be possible with higher sample path powers and 
with rather relaxed constraints to the reference path power. 


A technical solution to reduce noise and to remove constant parts in a signal is 
balanced detection. The following section briefly introduces this concept and 
points out its benefits. A balanced detector (BD) comprises two photodiodes 
PD1 and PD2 connected electrically antiparallel, as well as a transimpedance 
amplifier (TIA), see Figure2.8. The fields E,,= J2E,, E,o = j2E, 
returning from reference and sample path, respectively, in front of the combin- 
ing coupler, are denoted with pre-factors of 2, = W2 ‚since E, and E, are 
defined as the respective fields after the coupler (CPL2) and in front of PD1. 
The photodiodes PD1,2 generate the respective photocurrents ig) from the 
respective electric fields E,, and Ep entering the photodiodes. Due to the 
antiparallel connection of both photodiodes (PD1,2), the difference current 
Lig BD © laı — 4a2_ Can be tapped and amplified by a transimpedance amplifier 
(TIA) with real-valued transimpedance gain Z,, leading to the output voltage 


Usig BD = ~Za sig BD- (2.45) 
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Balanced detector 


Figure 2.8: Balanced detection concept with coupler (CPL2) and balanced detector. The input 
fields Eo = V2E, and E o=- WE, come from reference and sample path, respectively. A 
balanced detector comprises two photodiodes (PD1,2) connected electrically antiparallel, and a 
transimpedance amplifier (TIA). The photodiodes PD1,2 generate the respective photocurrents 
igi. In comparison to detection of the signal at one port of the interferometer, both ports with 
the respective electric input fields Æq; and Ey. are converted into currents. The antiparallel 
connection of both photodiodes (PD1,2) leads to a difference current iio BD iq) ig at the 
TIA input and to a voltage Usie BD = —Zylsig BD at the output of the TIA. 


The optical fields on both ports 


Ea = E, +E, (2 46) 
Ew =j(E, E,) 


are a superposition of the fields from sample and reference path. The photodi- 

ode currents are proportional to the squared absolute values of the respective 

2 2 
Eaz 


electric fields |En , which are calculated as 


’ 


ia ~ (Eal F Ef % JE T 2Ref EE) (2.47) 


in ~ [Evol IE + [ED -2Ref Ee, 


With the respective optical powers P., P., arriving from reference and sample 
path, and measured in front of each of the photodetectors, for optimum coher- 
ence GE -1,)=1), and for a single reflection in the sample path, which 
leads to a phase difference of Ap=@y (7, = t) at the optical angular frequency 
@, see Eq. (2.14), the photodiode currents in Eq. (2.47) can be expressed as 
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iu = p(B. + P, +2 JPP, cos(Ap)), a 
in = p(P. +R -2PR cos(Ag)). 


The subtraction in the balanced detector leads to a rejection of signal parts, 
which arrive with equal power on both detector inputs. Depending on the cou- 
pling ratio of the coupler and on the similarity of the photodiodes the rejection 
is not perfect but can be described by a common-mode rejection ratio (CMRR) 
¢ of the electrical signal powers. Using the CMRR, the signal current after 
balanced detection can be expressed as 


Iig BD = lat "lan + ei BE ig) 


lg T la 


(2.49) 


and can thus be described by the difference of the individual photodiode cur- 
rents ig; and ig), plus a small not-suppressed additive part {f (ig) +ia2). This 
analysis concentrates on the calculation of the OCT signal current isig gp after 
balanced detection, while noise contributions will be treated separately. Insert- 
ing Eq. (2.47) or (2.48) into Eq. (2.49) leads to 


lie, BD ~4Re| EE} + elle + EP) 


iig BD = P(4V PR cos(Ag) + 246 (P, + P,)), 


(2.50) 


where p denotes the responsivity of the photodetector. This shows the ad- 
vantage of balanced detection in OCT systems: The sample-reference interfer- 
ence term Re E E,| remains strong, while the parts that do not contain in- 
formation on the depth-dependent reflectivity of the sample are attenuated. 


For quantification of the dynamic properties of OCT measurements, three 
terms need to be defined: First, sensitivity is defined as a system property, ex- 
pressed in sample power reflectivity. The sensitivity quantifies the lowest 
measurable reflectivity of a sample and is linked to the noise floor of the 
measurement system in case the measurement signal is weak enough to not 
raise the noise floor itself, e.g., as a result of RIN. Second, the signal-to-noise 
ratio (SNR) relates the signal power level to the level of the noise floor. Third, 
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the dynamic range (DR) is a system property, describing the maximum 
achievable SNR. Usually, the DR is smaller than the total range of measurable 
signal strengths. The range of measurable signal strengths is limited on the 
lower side by the sensitivity and on the upper side by the smallest individual 
saturation value of a system component. Typically, ADC input range or pho- 
todetector saturation define the upper limit. By variation of the reference arm 
power, the DR can be shifted inside the total range, see Figure 2.9 (a). 


The lower limit of the measurement range is defined by the noise level. For a 
quantitative description, the following section is dedicated to an analysis of 
the individual noise terms in swept-source OCT systems. Figure 2.9 (b) shows 
the origin of the most prominent noise types in OCT: Relative intensity noise 
(RIN) from the swept source (SS), shot noise arising from the corpuscular na- 
ture of light, thermal noise from the amplification circuits in detector and 
ADC, and quantization noise from the digitization process in the ADC. 


OCT inter- 
detector 
ferometer 


RIN thermal noise 


sensitivity quantization 


shot noise no 
Figure 2.9: Sensitivity, dynamic range and noise sources in an OCT system. (a) Definition: 
sensitivity, dynamic range (DR), and total range. The total possible range of a measurement 
signal (expressed in sample power reflectivity R, light blue) is limited on the lower side by the 
sensitivity of the system and on the upper side by the saturation values of the individual system 
components like ADC or PD. The dynamic range (DR, dark blue) is within this range and typi- 
cally smaller. By variation of the reference arm power, the DR can be shifted within the total 
range. (b) Origin of the different noise sources in an OCT system: Relative intensity noise 
(RIN) from the frequency-swept source (SS), shot noise arising from the corpuscular nature of 
light, thermal noise from the amplification circuits in detector, and ADC, as well as quantiza- 
tion noise in the ADC. 


Relative intensity noise (RIN) describes noise originating from intensity fluc- 
tuations of a laser diode [46]. In swept-source lasers, these fluctuations do not 
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only come from random phase noise as in conventional continuous-wave la- 
sers. Instead, longitudinal mode hopping and gain competition between differ- 
ent longitudinal modes happen during the fast wavelength sweep and contrib- 
ute to the RIN [47]. RIN relates the intensity fluctuations to the average 
intensity and can be expressed as [44] 


RIN, = [RIN(f)df = (2.51) 

0 h 
with the one-sided power spectrum RIN( f ), which can be determined by a 
spectral analysis of the photocurrent in the baseband. To this end, the temporal 
fluctuations of the laser power are detected by a photodiode and spectrally an- 
alysed by an electrical spectrum analyzer. The spectral RIN( f ) is specified in 
units of Hz. The optical output power is given by A, and 


öR,(t)= Py(t)—- A(t) defines the mean-free power fluctuations over time t. 


In the following, we assume that with small electrical frequency ranges, the 
relative intensity noise spectrum is independent of the frequency, thus 
RIN( f ) =RIN. If RIN is known by measurement, the resulting second mo- 
ment of the noise current 


N 2 
lig! = pP? RIN P? Afa (2.52) 


is given by multiplication of the RIN with the optical power P, on the photo- 
detector, with the measurement bandwidth Af, and with the detector respon- 
sivity p. Expressed in signal current i, = pP, the relation follows as 


lien = RIN if Afa. (2.53) 


Shot noise is inherent in photon generation and detection systems and arises 
from the discrete random arrival times of the individual photons at the detec- 
tor. As a result, the generated photocurrent shows random fluctuations. Shot 
noise is a Poisson process and follows the Poisson distribution [44]. The re- 
sulting noise current can be expressed with its time averaged second central 
moment as 
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eee 
lin! = 2goisAfa (2.54) 


with the elementary charge qo, the average photocurrent is and the measure- 
ment bandwidth Af. The average noise power is proportional to the signal 
current, which is in accordance with the fact that the expectation value of a 
Poisson distribution equals its variance. 


Detector noise comprises thermal noise originating from the detector electron- 
ics. Often detector noise is specified in noise-equivalent power (NEP), see 
Eq. (2.44), expressing an assumed optical input power, which leads to a detec- 
tor output current with an SNR of one in the electrical bandwidth of 1 Hz. 
With the general expression for the photocurrent ig = PF 
the averaged second cen- 


dependent on de- 
tector responsivity p and optical input power F,,,, 
tral moment of the detector noise current inside the measurement bandwidth 


Af., can be expressed as 


il” = o?NEP?Af... (2.55) 
Quantization noise has its origin in the assignment of discrete quantization 
levels to the analogue input signal with continuous values. Although, strictly 
speaking, quantization is not a random process, the quantization error can be 
described by stochastic means and is part of this investigation. Dependent on 
the step size, the quantization error can be larger or smaller. This error is 
equally distributed within one voltage step OU, [48]. Because of the uniformly 
distributed quantization error n,, this results in a central second moment of 


the quantization error of' 


! The quantized signal Ug = uo +n, differs from the input signal u, by the equally distrib- 
uted quantization error n, € [- Ou, / 2, Ou, / 2]. Accordingly, the probability density func- 


tion inside the interval is Pa = 1 ug. The central second moment (variance) of the mean- 


free quantization error is then calculated to be 


—~ ou, /2 Su, /2 du Sul Su 
= Ir, = J my Pad =f in| == ae = 
£ ou, L3 -6u4/2 3u | 8 8 12 
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2 
=—, (2.56) 


For ease of comparison to other noise sources, the quantization error in volt- 
ages, Mq, is expressed as it was an effective photocurrent 1, =, / Z,» using the 
real-valued transimpedance Z, of the transimpedance amplifier. The resulting 
central second moment of the equivalent noise current is then 
given by 


= 2 
; Öug 


il EA 


(2.57) 


In the following, all noise contributions are determined for the case of bal- 
anced detection. 


In order to calculate the relative intensity noise contributions of the different 
signal currents arising from interference of sample path with sample path re- 
flections (subscript SS), reference path with reference path reflections (sub- 
script RR), and sample path with reference path reflections (subscript SR), the 
signal current in Eq. (2.50) is rephrased to 


Iig BD = 2p JEP, +2pJeP, +4p| P,P, cos(Ag). (2.58) 
— r enn end 
iss RR ise 


The corresponding second moments of the RIN-related noise currents are de- 
rived from Eq. (2.53) 


: 2 2 2 2 

liein.ss| =RIN iss Afa = 4° RIN SR Afo 

: 2 2 2 2 

lirIN.RR | =RIN ikr Afa =4p° RIN F Afa (2.59) 
; 2 ; 

liein,sr| =RIN isk Afa =16p° RIN P.P.Af, cos” (Ag) 


and sum up to the overall second moment of the electrical RIN 
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Ba = ce] k likINRR i + irın.sr j (2.60) 
= 4p°RIN| ¢ (P? + P?) +4P,P, cos? (Ag) | Afa. 


The second moment of the total RIN current is thus dependent on the optical 
powers returning from both reference and sample arm, however, for perfectly 
balanced detectors, only the mixing product of PP, plays a significant role. 
For typical detectors with ¢ > 0, the total RIN current also depends on the 
individual return powers and on ¢ . In most cases, in OCT, the reference arm 
returns a higher optical power level than the sample arm, thus the strong refer- 
ence-reference RIN contribution in the measurement signal can be suppressed 
by choosing a well-balanced combiner (CPL2 in Figure 2.8) and balanced de- 
tector such that they provide a good common-mode rejection (small ¢). 


The shot noise, see Eq. (2.54), of both photodiodes can be expressed as 


2 — 
isn. = 2qotaiAfe = 2qop(P. +R +2) RR cos (A) Afa 


(2.61) 
. [2 a 
in.| = 2g9ig2Afeı = 2qop(P. +R-2/RR cos(Ap))Af.ı 
and adds up in noise power, as both noise sources are independent, 
| 2 | Pal Ê 
lial” = isn. t in.| =4 qp(P, +P.) Afa- (2.62) 


In the OCT measurement, the variances of all noise currents add up to a total 
noise current variance 


Ta 4. 2 
la + lan + ln] + enf 
ou; 3 5 
Spt + p “NEP Afa +4qop (P, + P.) Afar (2.63) 


a 


+ 4p°RIN| c (PÈ +P?)+4PP cos? (Ap) | Af. 


In analogy to the noise power, the photocurrent signal power can also be ex- 
pressed dependent on the incident optical powers returning from sample arm 
and reference arm. Using Eq.(2.50) for the case of full coherence 
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KL (7, =t) =] between reference path light and light returning from one sin- 
gle sample reflection, as well as assuming perfect common-mode rejection 
(¢ =0) leads to a balanced signal current of 


iso pp =4P\V BP, cos(Ap). (2.64) 


For comparison to noise powers, the mean of the squared signal photocurrent 
is calculated to 


eb = Bann =(40/BP cos(Ap i) )) =16p° PR P. cos? (Ag). (2.65) 


This signal current is amplified by the transimpedance amplifier, see Fig- 
ure 2.8 and Eq. (2.45). For the calculation of the signal level, the upper signal 
limit of the transimpedance amplifier inside the photodetector and the limit of 
the ADC need to be taken into account. In the specific setup discussed in this 
p =+1.8 V output 
voltage, and the ADC is limited at max apc =+0.5 V when operated in the 


work, the amplifier of the photodetector saturates at Umax p 
low-noise amplifier-bypass mode. The lower of these limits (Umax apc) causes 
a corresponding maximum measurable photocurrent imax = Umax, ADC /Za> 
which limits the photocurrent to 


lie Jim = min {i lig, BD> ! hoary (2.66) 


This limitation is usually not important for moderate sample and reference arm 
powers, but plays a role when strongly reflecting elements are placed into the 
beam at perpendicular incidence and thus lead to strong sample arm return 
power. To keep the signal level below the saturation limit, either the transim- 
pedance gain Z, is reduced or an additional electrical attenuation is used. 


In the following, signal and noise contributions are evaluated basing on the 
aforementioned models and for the fibre-based SS-OCT system discussed in 
this work. Polarization effects are neglected and optimum coherence is as- 
sumed. Figure 2.10 shows the OCT signal and the various noise contributions 
as a function of reference arm power, when the sample reflectivity 
R,=-50dB is low. The vertical axis represents the mean squared photocur- 
rents on a logarithmic scale. The sample power of P, = 1.4 nW incident on the 
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photodetectors is derived from the output power of the swept source of 10 mW 
(+10 dBm), less 3 dB for the power split at each of the three passes through a 
coupler (CPL1 forward, CPL1 reverse and CPL2), less 50 dB at the sample, 
see Figure 2.7 for a sketch of the setup. Finally, the sample return power is 
subject to additional 9.5 dB measured coupling loss in the free-space optics 
part. The responsivity of the photodetector is o=1 A/W. The electrical signal 
power (black solid line) increases linearly with reference arm power, see 
Eq. (2.65). The ADC is not yet saturated at a sample arm reflectivity of 
-50 dB and for the reference arm powers investigated here. Quantization 
noise (brown curve) and detector noise (green curve) stay constant versus ref- 
erence arm power and are calculated using Eqs. (2.57) and (2.55), assuming an 
ADC voltage step size of du, =392 uV . The photodetector is specified with 
noise-equivalent power NEP = Tpw/ Hz, CMRR of ¢ =-25 dB and tran- 
simpedance gain of Z, =5000 V/A at a 50 Q load. The detection bandwidth 
is limited by an anti-aliasing filter to Af, =2 MHz. Shot noise power, 
Eq. (2.62), increases linearly with reference arm power (red curve), and thus 
stays in a certain fixed ratio to the signal power. The laser RIN amounts to 
1x107'° 1/Hz at f =2 MHz, and because of the limited CMRR and the weak 
sample arm return power, the P? term in Eq. (2.60) dominates the shape of 
the RIN curve (blue), which is thus quadratically dependent on the reference 
arm power. In total, all noise contributions sum up, Eq. (2.63), and are dis- 
played as the black dashed curve, limited by quantization noise for low refer- 
ence arm powers and by RIN for medium and high powers. 


The total noise is lowest for weak reference arm powers, however, at this point 
the measurement signal is weak as well. In order to find an optimum operating 
point of the system, the signal-to-noise ratio is calculated according to 


2 
l; . 
SNR = a n, (2.67) 
.2 2 
lial 
given by the ratio of averaged signal power uae lim Šima; and averaged noise 


power lin . Since the OCT signal is Fourier transformed by a M-point DFT, 
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the actual noise bandwidth is Af / (M / 2), thus decreasing the noise power by a 
factor of M /2, compare [30]. 


10" 

De Signal of -50 dB reflector 
ze} 

p : 

3 10 '*¥ Total noise 

3 Quantization noise 

p= 

£10" 
2 Detector noise 
a noise 

LQ 

2107 

ze} 

Cc 

oO 

Se 

510“ 

a 0 20 40 60 80 100 120 140 


Reference arm power P, at one diode [uW] 


Figure 2.10: Comparison of the OCT signal power and various noise contributions versus ref- 
erence arm power. All numbers are calculated for the specific fibre-based OCT system in Fig- 
ure 2.7. Micron Optics SS with RIN =1x 191° 1/Hz at f =2 MHz; balanced photodetector 
with responsivity of o=1 A/W, noise of NEP=7pW/VHz, CMRR of ¢=-25 dB and 
transimpedance gain of Z, =5000 V/A at 50 Q load; the ADC voltage step size amounts to 
Jug = 392 uV . The detection bandwidth is limited by an anti-aliasing filter to Af, = 2 MHz. 
The sample offers a reflectivity of R, =-50 dB, leading to a sample arm power on one photo- 
diode of P, =1.4 nW . The signal power (solid black) increases linearly with reference arm 
power and is plotted for a -50 dB reflector in the sample arm. Quantization noise (brown 
curve) and detector noise (green) are not dependent on the reference arm power. Shot noise 
(red) increases linearly with reference arm power (axis in log scale) and the RIN (blue) increas- 
es quadratically. The total noise (dashed black) is limited by quantization noise and RIN. 


Besides noise currents and SNR, the minimum measurable sample reflectivity 
R, is of interest. This quantity is equivalent to the noise floor at given sample 
and reference arm return powers and is expressed in terms of sample reflectiv- 
ities. At given sample reflectivity Rs of a first reflector, a certain SNR results, 
due to signal power dependent noise contributions. Thus, the minimum meas- 
urable sample reflectivity of a second reflector is given by 
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a's (2.68) 
” SNR’ ' 
This specification would be unnecessary, if the noise floor would be inde- 
pendent of the sample reflectivity, however, shot noise and RIN depend on 


received optical power and make this definition essential. 


The sensitivity S of the OCT system is defined as the minimum measurable 
sample reflectivity R,, at optimum sample and reference arm return powers, 


S=min{R,}. (2.69) 


Both noise floor R, and SNR are analysed in their dependence on sample and 
reference arm powers, see Figure 2.11. In Figure 2.11 (a), the colour-coded 
minimum measurable sample reflectivity R, (noise floor) in dB-units is plotted 
versus reference arm power P, (horizontal axis) and versus sample reflectivity 
Rs (vertical axis). The respective sample arm power Ps on one photodiode is 
plotted as second vertical axis. The sensitivity of the system of S =—100 dB 
is achieved for 16 uW reference arm power and sample arm reflectivities be- 
low —45 dB. At higher reference arm powers and especially at higher sample 
arm powers the noise level increases. The SNR of the OCT signal is given in 
Figure 2.11 (b) in colour-coded dB-units for varying sample and reference arm 
powers, denoted on vertical and horizontal axes, respectively. The SNR in- 
creases with sample arm power, when RIN from the reference arm is negligi- 
ble and the ADC is not saturated. The SNR reaches a plateau at a maximum of 
66 dB, defining the dynamic range (DR) of the system. 


Figure 2.11 (c-e) show the impact of all different noise contributions on the 
total noise floor R, in sample reflectivity units: The noise floor of the system 
is calculated in total and plotted together with the individual contributions of 
the various noise sources. Three specific sample reflectivities R, of 
-50 dB, —30 dB and 0 dB are considered, see Figure 2.11 (c), (d) and (e), 
respectively. 


At low and medium sample reflectivities, Figure 2.11 (c) and (d), quantization 
noise (brown curve) and detector noise (green) are constant. However, the 
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signal power increases with reference arm power, so the contribution of quan- 
tization noise and detector noise to the noise floor, expressed in sample reflec- 
tivity, decreases. At high sample reflectivities, Figure 2.11 (e), and above a 
reference path return power of P. =9 uW, the TIA gain needs to be reduced 
with increasing reference arm power, to avoid ADC saturation. This reduces 
thermal noise, whereas quantization noise stays constant. 


For a given sample reflectivity, the shot noise increases linearly with reference 
arm power, and the same is true for the signal power. The influence of shot 
noise on the noise floor in sample reflectivity hence stays constant as long as 
the ADC is not saturated. The situation is different for strong sample arm 
powers, Figure 2.11 (e), where the shot noise is not dominated by the refer- 
ence arm power anymore. In the regime of low reference arm powers, the 
squared sample arm power P? dominates the shot noise, leading to a shot 
noise contribution which decays with increasing reference arm power, as the 
signal power ~ P,P. increases. 


In the investigated setup, in most cases, RIN dominates the noise. For low and 
medium sample reflectivities, Figure 2.11 (c, d), RIN increases quadratically 
with reference arm power P,, while the signal power increases linearly with P,. 
As a consequence, the noise floor expressed in sample reflectivity units in- 
creases with the reference arm power. For high sample reflectivities, see Fig- 
ure 2.11 (e) and Eq. (2.60), RIN increases linearly with P,, and the same is 
true for the signal power, compare Eq. (2.65). The noise floor expressed in 
sample reflectivity units is hence independent of the reference arm power. 


The noise floor expressed in sample reflectivity units is lowest for reference 
arm powers of around P.=16uW, leading to the sensitivity of 
S=R,=-100dB for a sample reflectivity of R,=—50dB, see Fig- 
ure 2.11 (c). For a sample reflectivity of R =—30 dB, Figure 2.11 (d), a noise 
floor of R,=—95 dB can be achieved. At a strong sample reflectivity of 
Rs =0 dB, Figure 2.11 (e), the minimum noise floor is R, =—66 dB. 


A challenge in an OCT measurement is to reach both good sensitivity values 
and high SNR. Both goals can be achieved at a moderate reference arm power 
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around P.=16 uW on each photodiode. In this case, weakly reflecting sam- 
ples make OCT operation at the sensitivity limit possible, and strongly reflect- 
ing samples with R, >—25 dB still exploit the entire dynamic range of the 
system. 
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Figure 2.11: Minimum measurable sample reflectivity (noise floor) R, and signal-to-noise ratio 
(SNR) of the exemplary OCT system in Figure 2.7. (a) Colour-coded noise floor R, of the OCT 
system. The horizontal axis denotes the reference arm power P, incident on each photodiode, 
and the vertical axes refer to sample reflectivity Rs in dB units and the corresponding sample 
arm powers Ps on each photodiode. The sensitivity of S = —100 dB is achieved for P, = 16 uW 
reference arm power and sample arm reflectivities Rs below —45 dB (shaded area). (b) The 
SNR increases with sample reflectivity and reaches a plateau at 66 dB for sample reflectivities 
above approximately —20 dB. This maximum SNR defines the dynamic range (DR) of the sys- 
tem (shaded area gives full DR at ADC saturation). (c) Noise floor R, and individual contribu- 
tions to the noise floor R, in sample reflectivity units at a low sample reflectivity of 
Rs =-50 dB, plotted versus reference arm power. Quantization noise (brown line) and detector 
noise (green) have less influence with raising reference arm power, shot noise (red) stays 
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constant and the RIN influence increases linearly (note the log-scale on vertical axis). The total 
noise floor (black dashed) is limited by quantization noise at low reference arm powers and by 
RIN at high powers. A sensitivity of S = R, =—100 dB is achieved for reference arm powers 
of P.=16 uW. (d) Noise floor R, for sample arm reflectivity of Rs = —30 dB. The minimum 
measurable sample reflection is reduced to R, = —95 dB. (e) Noise floor R, for sample reflectiv- 
ity of Rs = 0 dB. RIN dominates the noise floor at Ra = -66 dB. From the ADC saturation limit 
around P. =9 uW onwards, the signal is clamped to the maximum range, e.g. by dynamic re- 
duction of the TIA gain Z, or by electrical attenuation. Here, the shot noise is dominated by 
strong sample reflected light and with reduction of the TIA gain the shot noise contribution is 
reduced, like the contribution of the detector noise. 


2.2  Polarization-sensitive optical 
coherence tomography 


This section first introduces the optical polarization and the Jones and Stokes 
formalisms, which serve for polarization calculations. Then polarization- 
sensitive OCT (PS-OCT) is introduced, followed by its typical application 
fields. Finally, two system variants, one with free-space bulk optical elements 
and one with fibre-coupled components are presented. In the following sec- 
tion, the underlining of complex variables is omitted for better readability. 


2.2.1 Polarization of light 


Light propagates as electromagnetic wave and its physical properties are de- 
scribed by Maxwell’s equations [49], linking the electric and magnetic field 
vectors E, H to the electric displacement D , the electric material polariza- 
tion P and the magnetic induction B. The current density and the space 
charge density are assumed to be zero. The vacuum permittivity is given by £ 
and the magnetic permittivity by 4. The relative magnetic permittivity is 
here assumed to be yw, =1 and the material is regarded as isotropic, linear and 
time-invariant. In the following, we assume complex sinusoidal time-domain 
quantities with the dependency ¥(7,t)=(7,@)e!”, where Y can be re- 
placed by E, H, D, B, P and Y can be replaced by E, H, D, B, P. The breve 
denotes the individual complex time-domain amplitudes. Accordingly, Max- 
well’s equations for the complex time-domain amplitudes Y (7,0) are 
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Vx H(F,@)=joD(F,o) Vx E(7,@) =-joB(7,a) 
V-D(7,a@) =0 V-B(r,o)=0 (2.70) 
D(7,@) = &E(7,@)+P(7,a) B(7,@)=H (7,0) . 


In linear media that are isotropic and homogeneous, the electric displacement 
D(F,@) can be expressed by the relative permittivity ¢,(@) =n" (œ), leading 
to D(F, o)= ££, (@)E (7,@), with the real-valued refractive index n(a). 
From Maxwell’s equations, Eq. (2.70), the wave equation for electromagnetic 
waves in homogeneous isotropic media can be derived, 


V7E(F,@) =-n(@) ky(@) E(F,0), (2.71) 


where co =1/ RENT is the vacuum speed of light”. An analogue expression 
exists for the magnetic field H (7,0). The wave equation can be solved by an 
equation describing the electric field of a wave, compare the scalar equation 
in Eq. (2.1), 


_ Z ; E ERT 

E(r.1)=E(r,0)ei® = Ae), (2.72) 
with the complex E-field amplitude vector A= A,é, + A,e, + A,e,, the propa- 
gation direction F=Xxe,+ye,+ze., and the wave vector 


k=k,e,+k,e,+k,e,. This wave propagates as a transverse wave, and the 


? Starting from V x E(F,@) =-j@B(7,@), and applying a (Vx)-operator to both sides of 
the equation, results in Vx(VxE(F, o)) =Vx(-j oB(F,o)), which can be rephrased to 
v(V-E(7,0))-V’E(r,0)=-jow(VxH). With no space charge, and in isotropic and 
homogeneous media V:D(F,0)=&8.(@)V-E(F,o)=0 holds. Further, 
Vx H(F,@) = joD(r,®) and D(7,@) = &é,() E(7,) lead to 
-V*E(7,@) =@ mE, (@)E(F,@) or, with ky(w) = = and co =1/,Jé os . is rephrased 


0 
to the wave equation V7E(7,@)= -n(o) ko (of E(F,o). 
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direction of its oscillation can be defined within a plane perpendicular to the 
direction of propagation [50] and is referred to as polarization’. 


The Jones and the Stokes formalism are conventional approaches for describ- 
ing the state of polarization of a light wave. This description includes changes 
of the state of polarization, which are induced by optical elements. In theory, 
both calculi could be applied for describing an optical non-depolarizing sys- 
tem, however in practice each of them is advantageous in different situa- 
tions [51]. The Stokes formalism is able to handle depolarized light, while the 
Jones formalism cannot. Stokes vectors are based on a series of intensity 
measurements after certain polarization-optical elements, Jones vectors calcu- 
late from electromagnetic theory. A graphical representation of the Stokes pa- 
rameters is the Poincaré sphere. However, only the Jones calculus allows di- 
rect access to the optical phase. 


The Jones vector J is normalized, J - J "= 1, and describes any fully polarized 
electric field that consists of two complex field amplitudes oriented along the 
x and y-direction, having a common phase 6, a phase difference 2¢, and ampli- 
tudes cos (a), sin(@), 


-jọ 
-hye Daa (2.73) 


etl? sin(a) 


With the help of the Jones vector and its components Jy, Jy, and the scalar elec- 
tric field amplitude A, any fully polarized plane wave, here propagating in z- 
direction, can be expressed as 


* with = V-D(7,@)=0, and thus V-E(7,@)=0, it follows that 


2% el (orka) +Ê A de) +Ê Aei (0ta) =0. Applying the derivation leads to 
ax * dy ° oz ` RR 
JKR,E,-jk,E,-jk,E, =0, which is equivalent to —j k : E =0. Hence, the electric 


field E of a plane wave is transverse to the direction of propagation k. This holds analo- 


gously for the magnetic field. 
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i _ : ; L ; J 
E(F,t)=E(F,@)e!” = Ae 2E) F = geiler) i ; (2.74) 


Polarization describes the relationship of the x and y-component of the com- 
plex wave vector: Linearly polarized light is defined by equal phases in the 
two orthogonal components, which have a phase difference of 29=0. Circu- 
lar polarization is observed, if the phases are different by zj 2 (p=+r/4) and 
the amplitudes are same. If the relative vertical amplitude J, is zero and the 
relative horizontal amplitude J, is non-zero, the wave is horizontally polar- 
ized — vertical polarization is defined analogously with J,=0, J, #0. Any 
intermediate form with non-equal amplitudes and phases is called elliptical 
polarization. 


Stokes vectors, in contrast to Jones vectors, contain four real elements, repre- 
senting the total intensity /,,, and the intensities measured with a horizontal 
polarizer por» a vertical polarizer /,.,, as well as intensities measured with a 
+45°- and a +135° -linear polarizer 145, 1,35, and the intensity measured with 
a right-circularly and a left-circularly polarizing filter [p7, I,z. According to 
these measurements the Stokes vector S, ‚ results in 


Sy Liot 

= Ky IL -I 

S, > 2: x hor ver (2.75) 
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Both the Stokes and the Jones formalisms allow to describe changes in the 
state of polarization by using matrices that transform the input polarization 
vector to the output vector. In case of multiple optical elements passed, the 
matrices have to be concatenated by matrix multiplication. In the Stokes cal- 
culus, the transformation matrices are named Mueller matrices, in the Jones 
calculus, Jones matrices. Accordingly, light of polarization J passing optical 
elements with Jones matrices J; „1s transformed to light with Jones vector 


Jn =I SJ. (2.76) 
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Expressed in Stokes vectors, light of polarization Sı passes elements with 
Mueller matrices M, „ and is transformed to 


S,=M, MS. (2.77) 


Both calculi are linked and Mueller matrices can be derived from Jones matri- 
ces [52]. The reverse calculation of Jones matrices from Mueller matrices only 
succeeds in the case without depolarization. 


The purpose of polarization considerations within this work is the description 
of polarization effects in optical media including the system itself and the 
sample under test. As OCT measurements record only the coherent part of the 
backscattered light, depolarized light will not be detected. Hence, the Jones 
calculus is sufficient for the description of PS-OCT and it allows direct access 
to the horizontal and vertical field components. 


2.2.2 PS-OCT concepts 


The very general aim of a PS-OCT measurement is the depth-dependent de- 
termination of polarization properties of a sample [53]. Depending on the ap- 
plication, it suffices to only determine birefringence [54], dichroism, the opti- 
cal axis orientation [55] or a set of Stokes parameters [56] for a single incident 
polarization. This leads to simplifications in the optical setup and reduced 
measurement time, as compared to more comprehensive PS-OCT systems 
with capability of extracting the full Jones [57] or Mueller [58] matrix of the 
material. Early systems were based on conventional TD-OCT approaches [53], 
whereas later polarization-sensitive systems also adopt the SD-OCT or SS- 
OCT concept [59]. 


2.2.3 Applications of PS-OCT 


PS-OCT systems can serve for various applications. In dental OCT, buried 
caries lesions can be detected by PS-OCT, as demineralization accompanying 
carious defects leads to anisotropically ordered structures, which cause cross- 
polarized backscatter signals [60]. Dermatologic applications comprise analy- 
sis of burn depths that can be determined by characterization of the 
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birefringence difference between unaffected skin layers and layers with former 
heat impact [61]. This is possible, since collagen proteins in the skin are origi- 
nally birefringent, but denaturate from temperatures of around 60°C onwards, 
while losing their polarization properties. PS-OCT is further applied to the de- 
tection of basal cell carcinoma, which lead to a local alteration of skin struc- 
ture and thus to local changes in dermal birefringence [62]. Recent studies 
combined this analysis with a machine-learning algorithm and achieved a de- 
tection sensitivity and specificity of 93 % and 87 %, respectively [63]. An im- 
portant application field of PS-OCT is ophthalmology, aiming at, e.g., early- 
stage detection of glaucoma. This degeneration of the retinal nerve fibre layer 
is in some cases preceded by a necrosis of ganglion cells inducing a replace- 
ment of nerve axons by tissue of a lower degree of organization. Hence, the 
birefringence of the retinal nerve fibre layer changes. A PS-OCT measurement 
can detect this birefringence modification and thus allows for early diagnosis 
of glaucoma [64]. As a use case in material sciences, PS-OCT is applied to 
determine mechanical stress in glass/epoxy composites, where the birefrin- 
gence of the material is stress dependent in both orientation and magnitude 
[65]. The PS-OCT analysis offers the prospect of early failure detection. 


2.2.4 Measurement principle 


A free-space optical setup of a PS-OCT system is depicted in Figure 2.12. 
With a single measurement of both received polarizations, this setup is capa- 
ble of determining the birefringent properties of the sample [53]. 


In contrast to the conventional free-space TD-OCT, see Figure 1.2, the inci- 
dent light is horizontally polarized by a polarizer (Pol) before being split 
among sample and reference arm. The part towards the sample passes a quar- 
ter-wave plate (QWP) at 45°, which leads to circularly polarized light incident 
on the sample. This can be seen from the following analysis: 


A QWP with its fast axis in the horizontal direction is described by the Jones 
matrix [66] 
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im(1 0 
Jow =e 4 | ) (2.78) 
Using the rotation matrix [67] 


SOE M a i (2.79) 


—sinY cosg 


the optical element, here the rotated QWP with fast axis at an angle 9 to the 
horizontal direction, can be expressed by the Jones matrix 


Jowr (9) =R(-3)JowrR (3) 
i; { cos” 9+ jsin’ 9 (1—j)sin cos 9 (2.80) 


=e i 
(1- j)sin 9cos9 sin? 9+ jcos? 9 
Mirror f 
L——QWP (22.5°) 
Broadband 


QWP (45°) Sample 


Photo detector 


Figure 2.12: Free-space polarization-sensitive OCT setup. Pol: Polarizer, BS: (non-polarizing) 
beam splitter, QWP: quarter-wave plate at given orientations of the fast axis with respect to the 
horizontal direction, PBS: polarization beam splitter. The setup is a modification of the conven- 
tional TD-OCT setup shown in Figure 1.2. The sample comprises various reflections at depths 
Zs,15 Zs2 and zs3. The polarizer lets horizontally polarized light pass, and the sample is illuminat- 
ed by circular polarization. The light passing the 22.5°-QWP in the reference arm forth and 
back returns to the BS linearly polarized at —45°. The interference signal is split at the PBS in 
horizontal and vertical polarization and is detected. 
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With = 45° this results in 


1 Zfl+j 1-j 
45°)=—e 4 : 2.81 
Applied to horizontally polarized input light, as is the case here, the output 
polarization travelling towards the sample is 


Telen 1+j 1-j\/1 -l7 1+j _v2 j (2.82) 
73. (-j ae Ol 2° Se 2U 


which represents a circular polarization state, characterized by equal ampli- 
tudes in each of the orthogonal axes and 90° optical phase delay. Circular po- 
larization is of advantage over linear polarization, as an arbitrary linear polari- 
zation could coincide with the optical axis of the sample. In this unwanted 
case, no birefringence information could be extracted, as the difference in 
propagation constants of slow and fast axes would not be measurable. 


The light returning from the reference arm has passed a QWP at 22.5° in for- 
ward and backward direction, which results in —45° linearly polarized light, as 
is explained in the following calculation: 


The reference arm represents a concatenation of QWP, plane mirror and QWP 
in reverse direction of propagation, leading to the total Jones matrix of the ref- 
erence arm 


Jet (9) =Jowrrev (A) Iu owe (I) $ (2.83) 


The Jones matrix J QWP (3) of a rotated QWP has been evaluated in Eq. (2.80) 
and is unitary, as a QWP represents a loss- and gainless element [68] leading 
to an expression of the type 


Jon(9)=[ a a (2.84) 


Un u 


with the matrix elements 


56 


2.2 Polarization-sensitive optical coherence tomography 


uy = et (cos? +j sin? 9), 
(2.85) 


u = aa ((1- j)sin 9cos 9). 


We choose a right-handed Cartesian coordinate system, with the +z-direction 
being always parallel to the direction of propagation of light. Hence, when 
light is incident on a plane mirror, the reflected light is described by a coordi- 
nate system rotated about the x-axis, and the Jones matrix of a plane mirror at 
perpendicular incidence becomes [68] 


In -f N (2.86) 


The Jones matrix Jey of reverse propagation of any optical reciprocal materi- 
al with Jones matrix J , when taking into account a coordinate system as de- 
scribed above, can be written as [69] 


Iren -( a a with J -(™ >) (2.87) 
-m m, m; mą 
In case of a unitary reciprocal material like the QWP, compare Eq. (2.84), this 
leads to 
u u, 
JawP.rev = r . . (2.88) 
Hy M 


Now, the total Jones matrix of the reference path can be calculated from 
Eq. (2.83), 


ier (3) = Jawp.rev (3)IuJowr (3) 


2, 2 * 
B uy +u; uu, — (uuz ) (2.89) 
p * *2 2 ' 
=u + (uu) u = uy 


7 | cos(29) a 


= -sin(29) cos(29) } 
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This combination of forward-pass of the QWP, the mirror and the reverse pass 
of the QWP acts similarly to a half-wave plate (HWP) with 


Jaw (9) = no a , (2.90) 


where the second row of the Jones matrix is negated as compared to Jef (3). 
With the fast axis of the QWP at 4=22.5° with respect to the horizontal axis, 
the overall reference path matrix is 

v2 


Jer (22.5°) = a (2.91) 


Applied to horizontally polarized input light, this results in 


rol E am 


referring to as linear polarization at —45°. Hence, with equal amplitudes in 
horizontal and vertical direction, the reference field could interfere with both 
horizontally and vertically polarized components of the field backscattered 
from the sample. The polarizing beam splitter in the detection path splits both 
polarizations and guides them to separate detectors. Birefringence of the sam- 
ple would emerge as oscillation of signal strength between the two detection 
polarizations versus measured sample depth. With rotation of the sample or of 
the optic elements in the system, the optic axis of the sample could 
be determined. 


In free-space optical systems, the polarization can be controlled exactly at any 
point of the system but is sensitive with respect to mechanical misalignment. 
In contrast to that, fibre-based systems are more robust, but propagation of 
light in single-mode fibres (SMF) influences the optical phase in both orthog- 
onal polarizations in an undetermined way and thus modifies the polarization 
state. In case of commonly used standard SMF, this impact is highly depend- 
ent on the exact positioning of the fibres, fibre bend stress and temperature. 
Thus, without calibration, the Jones matrix of a fibre is unknown. The use of 
polarization-maintaining fibres could circumvent this calibration step, as the 
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incident polarization state can be maintained. This polarization-maintaining 
behaviour is possible by minimizing cross-talk between signals in the slow 
and the fast axis of the fibre, which is achieved by stress-induced differences 
in the individual propagation constants. 


In OCT practice, polarization-maintaining fibres are difficult to deploy, since 
the strong group delay differences between the slow and the fast axis leads to 
ghost images at a relative amplitude, which corresponds to the polarization 
cross talk, which may, e.g., occur at each fibre-to-fibre connection. From the 
moment of polarization mode coupling, each portion of the respective polari- 
zation travels with a distinct speed and arrives at the detector with a respective 
time delay. This is equivalent to different imaging depths in the sample. With 
these ghost images, OCT signal analysis would become severely disturbed. 
For high-end PM connectors, cross-talk can be kept smaller than -30 dB. Still, 
with dynamic ranges of more than 60 dB, this is insufficient for many OCT 
applications. As a consequence, the use of polarization-maintaining fibres is 
avoided and we use standard single-mode fibres instead. A polarization cali- 
bration process is then mandatory. 


A fibre-based system faces two challenges: First providing a method for char- 
acterization and compensation of the SMF-induced polarization change, and 
second, supplying the sample with incident light, with polarization different to 
the principal state of the sample medium. The latter property is important in 
order to measure material birefringence by observing a difference in propaga- 
tion constants of the material. For a more comprehensive analysis of the mate- 
rial, e.g., when Mueller or Jones matrices are to be measured, the system 
needs to be able to vary the state of polarization incident on the sample. 


In the following we present a system concept, which satisfies the abovemen- 
tioned conditions for a fibre-based Jones matrix PS-OCT, see Figure 2.13. 
Light emitted by the frequency-swept light source (SS) is sent through a first 
polarizer, thus guaranteeing a constant input polarization during one meas- 
urement. This light propagates either towards the sample measurement head or 
to the calibration setup and leaves the fibre section in the sample path (SP) 
with polarization J in. Lhe portion of the input light travelling through the ref- 
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erence path (RP) is adjusted in polarization, such that it later splits in equal 
parts on the two axes of the PBS in the detection path. This is controlled by 
monitoring the received power at each of the single photodiodes. This setting 
ensures that a received sample field can interfere with the reference field, in- 
dependent of its polarization. When a sample is characterized, the sample is 
illuminated with a set of different polarizations, arising from different rota- 
tional positions I of the fast axis of the QWP (QWP 2) with respect to the 
horizontal axis (0° to 165° with a step size of 15) with corresponding Jones 
matrices J owr2(%)- In the sample, the light undergoes a depth-dependent 
polarization change J s(z). This Jones matrix incorporates forward propaga- 
tion to the backscattering plane at depth z, the scattering process itself, and 
backpropagation. Dependent on the material to be investigated, either material 
propagation, e.g., in birefringent material like skin, or the scattering processes, 
e.g., at elongated nanoparticles, dominate the Jones matrix. Backscattered 
light exiting the sample towards the fibre-optic system reversely passes the 
QWP 2. The reverse propagation leads to the Jones matrix J owp2,rev (%) , see 
Eqs. (2.88) and (2.85). The following fibre-optic propagation back through the 
SP and then towards the PBS in the receiver arm is taken into account by the 
matrix J,,,. Finally, the Jones vector as incident on the PBS, and thus meas- 


ured by the detector is 


I = J sys Qwe2,rev (% )I.J ows (% ) Jin £ (2.93) 


The matrices J owr2(%) and J QwP2rev(% } are defined by the rotation an- 
gle 5 of the QWP 2. The matrix of interest J, is unknown, as well as the 


system matrix J,,, and the Jones vector J,, of the input light. Assuming J 


sys sys 


and J;, would be known, J, could be determined from a set of measurements 
of J,, with different QWP-settings. However, both J 
the unknown polarization transformation in the fibre optic system, which 


and J,, depend on 


sys 


makes a preceding calibration step necessary. Once both are determined, the 
calibration holds until fibres are moved or the temperature changes. 
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Figure 2.13: Fibre-based swept-source polarization-sensitive OCT setup. SS: Swept source, 
Pol 0,1: Polarizer, CPL: 3-dB fibre coupler (polarization-insensitive), SP: sample path, 
RP: reference path, FC: fibre collimator, QWP 1,2: quarter-wave plates with Jones matrix 
J QwPI.2 (3): PBS: polarization beam splitter, PD: photodetector, PolC Polarization control- 
ler. Jones vectors and matrices are denoted in blue colour. Jj, : Jones vector of the light exiting 
the fibre section towards the sample, J poi (% ): Jones matrix of the polarizer in the calibration 
setup, Js: Jones matrix of the sample, Jsys: Jones matrix of the fibre optic system connecting the 
FC in the sample path to the PBS, J m: measured Jones vector at the PBS. For polarization 
calibration a calibration setup consisting of a rotatable QWP (QWP 2), a polarizer , and a mir- 
ror (M2) is used in the sample arm. 


In order to perform this calibration, the sample is replaced by a calibration 
setup, comprising a rotatable polarizer and a plane mirror, see Figure 2.13. By 
rotation of the polarizer POL 1 to the angular positions 9, =0° to 150° in 
steps of 30°, and by rotation of the QWP 1 to the angular positions Jg =0° to 
165° in steps of 15°, several polarization states can be generated and the fibre- 
optic system can be probed with these. During this calibration, the Jones vec- 
tor received by the PBS is given by 
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Tii = I ysJ wei .zev (% )Ipottrev (%)IuJ pon ($)Jawrı (% ) Jin > (2.94) 


with the Jones matrices Jpo (9p) of the polarizer, Jpoij rev (Ip) of the polar- 
izer in reverse direction of propagation, and the Jones matrix Jņ of the plane 
mirror. The mirror can be replaced by any strongly reflecting surface, e.g. by 
the surface of a glass cuvette containing the sample. The determination of J sys 
and J,, results from a least-squares fit of the system model in Eq. (2.94) for 
several angular settings of the QWP | and of the polarizer. 
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Nanocomposite materials represent a success story of nanotechnology. How- 
ever, development of nanomaterial fabrication still suffers from the lack of 
adequate analysis tools. In particular, achieving and maintaining well- 
dispersed particle distributions is a key challenge, both in material develop- 
ment and industrial production. Conventional methods like optical or electron 
microscopy need laborious, costly sample preparation and do not permit fast 
extraction of nanoscale structural information from statistically relevant sam- 
ple volumes. Here we show that optical coherence tomography (OCT) repre- 
sents a versatile tool for nanomaterial characterization, both in a laboratory 
and in a production environment. The technique does not require sample prep- 
aration and is applicable to a wide range of solid and liquid material systems. 
Large particle agglomerates can be directly found by OCT imaging, whereas 
dispersed nanoparticles are detected by model-based analysis of depth- 
dependent backscattering. Using a model system of polystyrene nanoparticles, 
we demonstrate nanoparticle sizing with high accuracy. We further prove the 
viability of the approach by characterizing highly relevant material systems 
based on nanoclays or carbon nanotubes. The technique is perfectly suited for 
in-line metrology in a production environment, which is demonstrated using a 
state-of-the-art compounding extruder. These experiments represent the first 
demonstration of multi-scale nanomaterial characterization using OCT. 


3.1 Introduction 


Nanomaterials represent an emerging multi-billion dollar market driven by a 
vast variety of applications that range from mechanical and civil engineering 
to energy storage and life sciences. Examples comprise nanocomposite poly- 
mers with enhanced mechanical or electronic properties [11,12], functional 
coatings [70], flame-retardant materials [71], advanced drug-delivery sys- 
tems [13], and anode materials for Li-ion batteries [72]. These applications 
mostly rely on the unique properties of nanosize particles, namely huge sur- 
face-to-volume ratios, enhanced tensile strengths and superior electrical con- 
ductivities as shown by carbon nanotubes (CNT) or other nanofibers [73], or 
outstanding barrier properties of nanoplatelets [74]. Properties of 


64 


3.1 Introduction 


nanocomposites depend not only on the size and shape of the particles but also 
on their dispersion state, characterized by the degree of agglomeration when 
immersed into a liquid or solid host material. The dispersion state is governed 
by the nanoparticle properties, by the composition and the physical parameters 
of the host material, as well as by the processing route adopted for dispersing 
the nanoparticles in the host material. In order to ensure consistent product 
quality, the dispersion state must be continuously monitored during fabrica- 
tion, which has been identified a key challenge both in nanomaterial develop- 
ment and industrial production [18]. In solid media, dispersion-state character- 
ization mainly relies on microscopic imaging. For small nanoparticles, this 
requires high-resolution techniques such as scanning-electron microscopy 
(SEM) of specially prepared sample surfaces [19], or transmission electron 
microscopy (TEM) of microtome slices [20]. In both cases, sample prepara- 
tion is laborious and costly and not well suited for quality control or process 
development, where sample processing and analysis have to be continuously 
iterated. Moreover, SEM and TEM are limited to small sample volumes, 
which are not necessarily representative of the entire batch. Light microscopy, 
on the other hand, can reduce the experimental effort [21], but is limited to the 
identification of large agglomerates in sufficiently transparent samples and 
areas close to the surface. In contrast to image-based methods, light scattering 
techniques have proven to be viable tools for measuring particle size distribu- 
tions from large sample volumes. Static light scattering (SLS) relies on angle- 
resolved and/or spectrally resolved detection of scattered light [22,23] and is 
based on a rather complex optical setup, especially when large scattering an- 
gles have to be taken into account. Moreover, for solid samples, it is challeng- 
ing to separate scattering inside the sample volume from contributions of the 
rough sample surface. Dynamic light scattering (DLS) methods exploit tem- 
poral fluctuations of interference patterns of scattered light to calculate the 
Brownian motion and the hydrodynamic diameter of particles within the re- 
spective solvent [75]. This technique is limited to low concentrations, and can 
only be applied to liquid media. In addition, both SLS and DLS suffer from 
limitations in dealing with mixtures of particles having vastly different diame- 
ters that may range, e.g., from a few nanometres to hundreds of micrometres, 


65 


3 Multi-scale dispersion-state characterization of nanocomposites using OCT 


as is often the case for nanocomposites with poorly dispersed nanoparticles. 
Partial wave spectroscopic microscopy (PWS) [76] has been used for investi- 
gation of nanoscale refractive index fluctuations, which can be an early indica- 
tor of carcinogenesis, yet without providing a link to the size of the scatterers. 
X-ray diffraction (XRD) [77] or small-angle X-ray scattering (SAXS) [78] 
finally rely on diffraction or scattering of X-rays in the sample and are able to 
reveal the atomic or molecular arrangement inside the material. However, high 
instrumental effort and a small probing region limit the application range of 
these methods to the laboratory environment. Hence, none of the aforemen- 
tioned techniques can meet the stringent requirements associated with indus- 
trial process development and quality control, which comprise robustness of 
the measurement method, fast analysis, the capability to characterize repre- 
sentative sample volumes without laborious sample preparation, good me- 
chanical and thermal robustness of the measurement system, and the possibil- 
ity to integrate the measurement system into the processing line for enabling 
in-line process control. The lack of adequate analysis methods for dispersion- 
state characterization is considered one of the major obstacles towards large- 
scale industrial processing and exploitation of nanomaterials. 


As an alternative, optical coherence tomography (OCT) [1,28,79] was pro- 
posed as a tool for nanocomposite characterization [80,81]. OCT provides 
three-dimensional imaging data from the bulk of the sample and avoids expen- 
sive sample preparation. However, previous demonstrations have been limited 
to image-based analysis of composites containing rather large microparticles 
combined with wavelet-based processing of the image data [81]. Due to the 
limited resolution of OCT imaging, this method cannot provide information 
about the nanoscopic structure of the nanocomposite such as the particle size. 
A further method uses a spectral-domain OCT setup and measures the size of 
particles suspended in a fluid by the temporal variation of the optical phase 
due to particle diffusion [82,83]. Particle sizes range from 15 nm to 625 nm. 
The main disadvantage is that only particles suspended in a fluid can be inves- 
tigated, and that parameters like viscosity and temperature need to be tightly 
controlled. Another technique is angle-resolved OCT [84], where the angular 
dependence of the scattered intensity is evaluated according to Mie’s theory. 
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However, the angular range is limited to about 0.5 rad due to practical reasons, 
and the particle size which can be detected is not smaller than 5 um in diame- 
ter. A further approach is low-coherence spectroscopy, which allows the ex- 
traction of wavelength-dependent scattering coefficients of the investigated 
samples and compares well with Mie scattering calculations [85]. The afore- 
mentioned techniques serve well for the determination of scattering parame- 
ters from nanoparticle samples, but miss the multi-scale approach, including 
imaging of large agglomerates. A technique using a setup similar to spectral- 
domain OCT is super-resolution imaging relying on spectral encoding of spa- 
tial frequency (SESF) [86]. With that approach, sub-micrometre imaging has 
been demonstrated, but exact nanoparticle sizes cannot be extracted due to the 
still limited resolution. In summary, OCT-based quantitative and qualitative 
characterization of composites at the nanoscale still remains to be shown. 


In this paper we demonstrate that OCT represents an attractive tool for fast 
and robust dispersion-state characterization of composite materials over a 
wide range of particle and agglomerate sizes both on the micrometre and on 
the nanometre scale. The technique pursues a multi-scale approach: Using a 
theoretical model of light scattering in the sample, we accurately measure par- 
ticle sizes down to 140 nm. Particle agglomerates with sizes of up to hundreds 
of micrometres can be easily detected by applying dedicated image processing 
techniques to the OCT data. Both methods can be performed in situ, without 
prior sample preparation, in both liquid and solid materials, and are applicable 
to laboratory investigations as well as to in-line process control in industrial 
production. Using a model system of polystyrene nanoparticles dispersed in 
water, we prove the reliability and accuracy of our sizing technique. We fur- 
ther apply the technique to an epoxy resin filled with multi-wall carbon nano- 
tubes (MWCNT). The results of OCT-based scattering analysis show good 
correlation with independently measured material properties, thereby outper- 
forming conventional characterization techniques based on light microscopy. 
Finally, we show that our technique is also perfectly suited for in-line metrol- 
ogy in a production environment. To this end, we integrate our system with a 
state-of-the-art industrial compounding extruder using a dedicated optical 
probe. The OCT system operates reliably during the compounding process and 
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allows to immediately examine the impact of extruding parameters on the dis- 
persion state of the material. We believe that OCT will pave the path towards 
industrially viable nanomaterial characterization and process control. 


3.2 Materials and Methods 


3.2.1 Swept-source optical coherence tomography system 


Optical coherence tomography (OCT) evolved greatly in the past decades. The 
technique provides microscopic resolution in volumetric imaging and highly 
sensitive detection of backscattered optical power. OCT opened a wide field 
of applications reaching from ophthalmology in medical diagnostics [87] to 
particle and defect characterization in material sciences [80,88]. Among vari- 
ous implementations, swept-source OCT (SS-OCT) offers a particularly at- 
tractive combination of highest sensitivity and high imaging speed [28]. 


The SS-OCT setup used in this work is depicted in Figure 3.1 (a). In general, 
an OCT system measures the position and the strength of a multitude of scat- 
terers along a light path in a sample. To this end, the electric field that is 
backreflected from a sample is compared in amplitude and phase to a refer- 
ence field. Both the sample and the reference field are derived from the same 
optical swept-wavelength source (SS). In our experiments, we use an SS laser 
with central wavelength of 1315 nm and a wavelength scanning range of 
1260-1370 nm, a scan rate of 1 kHz, and 10 dBm average output power 
(model s3, Micron Optics Inc., Atlanta, GA, USA). The scan range of the laser 
corresponds to a theoretical depth resolution of 7 um, which compares well to 
the experimentally observed resolution of 11 um. A first fibre-based direc- 
tional coupler (CPL1, splitting ratio 50:50) is used to split the power among 
the sample path (SP) and the reference path (RP). The reference path contains 
a free-space section allowing precise matching of the RP and SP length. 
Backscattered light from the sample and light travelling along the RP is super- 
imposed in a second fibre-based coupler (CPL2, splitting ratio 50:50) and 
coupled to a balanced photodetector (BD, model PDB430C, Thorlabs, Mu- 
nich, Germany). The output current of the BD contains patterns resulting from 
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interference of the backscattered field with the reference field. The BD sup- 
presses intensity noise from the strong RP signal and enables a large dynamic 
range and a high sensitivity, defined by the lowest detectable backscatter from 
the sample, of —110 dB. The electrical signal is digitized by a 16 bit analogue- 
to-digital converter (ADC, model ATS660, Alazar Technologies Inc., Pointe- 
Claire, Canada) and processed in a personal computer (PC). The amplitude 
and the position of the backscatter along the light path can be obtained by 
Fourier analysis of the photocurrent as a function of optical frequency [28]. 
One depth-scan (A-scan) consists of 768 measurement points with 8 um step 
size. We extract the backscatter strength, which denotes the ratio of backscat- 
tered optical power to optical power incident on the sample. 


In the course of this work, the SP of the SS-OCT system is connected to two 
different scan heads: First, to a conventional OCT scan head for off-line char- 
acterization of laboratory samples, Figure 3.1 (b), and second, to a specially 
developed probe head for in-line dispersion characterization during the com- 
pounding process in an industrial extruder, Figure 3.1 (c). This probe is de- 
signed for the harsh environmental conditions at a nanocomposite production 
line and must tolerate vibrations, high temperatures of 250°C, and high pres- 
sures of 200 bar. The probe features a titanium shaft and a sapphire window 
towards the sample, and is designed to allow the adjustment of focal length 
and the axial position of the focus within the sample. 


3.2.2 Model-based dispersion-state analysis and sizing 
of nanoparticles 


Big agglomerates of nanoparticles with dimensions larger than the resolution 
limit of the OCT system can be detected directly by imaging. However, size 
information on nanoscale agglomerates and single particles is relevant as well. 
In this section we show that model-based analysis of OCT backscatter meas- 
urements allows to extract scattering parameters that are correlated with the 
dispersion state of the material so that even small particle sizes can be deter- 
mined. 


69 


3 Multi-scale dispersion-state characterization of nanocomposites using OCT 


a Conventional scan head b 
ss a ent a 
= — Scan lens 55 

- —- 

x, — ) ] 
—— 

FC y —— j 

Celt 6p z Scan mirrors § 


Extruder probe 


Sapphire c 
Adjustable optics Ere 


fn Titanium shaft 
= Adjustable optics Sapphire window 


Figure 3.1: Fibre-based swept-source optical coherence tomography (SS-OCT) setup for labor- 


Sample 


Ed 


atory measurements and in-line process control. (a) Schematic setup of the SS-OCT. SS: 
Swept-wavelength source, scanning range 1260-1370 nm; CPL1,2: fibre-based directional cou- 
pler; SP: sample path; RP: reference path; FC: fibre collimator; PolC: polarization controller; 
BD: balanced photodetector; ADC: analogue-to-digital converter, PC: personal computer. 
(b) Laboratory scan head with galvo-based scanners and a scan lens for 3D-imaging. The red 
line indicates the light beam towards the sample. (c) Extruder probe for in-line measurements. 
The probe consists of a titanium shaft for a high temperature environment and comprises a sap- 
phire window towards the sample. 


The analysis relies on a single-scattering model assuming that incident light is 
scattered at maximum once within the medium, similar to the approach used 
by Kodach et al. [89]. This model is found to be appropriate for weakly scat- 
tering samples and for an analysis of moderate scattering depths [90]. As de- 
picted in Figure 3.2 (a), a light beam with input power Pin enters the sample. 
At each particle (grey), a first portion (blue) of the incident light in the respec- 
tive depth z is scattered back into the numerical aperture of the optical system, 
a second portion (red) is scattered into all other directions or is absorbed, and 
the remaining third portion (black) is propagating deeper into the medium. 
The total scattering o, and absorption cross section o, can be described by the 
extinction cross section o, =0,+0, of a single particle, or by the extinction 
coefficient 4, = No, for an ensemble of particles with volume number density 
N. In analogy, the backscatter is described by the backscatter cross section oj, 
or by the backscatter coefficient 44, = Noy, of the material. 
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The depth-dependent decay dP(z ) / dz = -4P (2) of the forward-travelling 
power P(z) is dictated by the extinction coefficient 44, leading to an expo- 
nential decay of optical power P(z)=P,e “* inside the sample. For a given 
depth z, the power in a depth element öz as scattered back into the numerical 


TMs 


aperture of the optical system amounts to P,(z)= Mp özP(z)e , where 


e “® accounts for the extinction of the backscattered light during backpropa- 
gation. The signal measured in an OCT scan corresponds to the depth- 
dependent backscatter P (z) / Pp. In the presence of a depth-independent 
noise floor R,, the depth-dependent backscatter signal measured by the OCT 


system is 
R,(z)= 1, 8ze +R. 68.1) 


A semi-logarithmic plot Ryg (z)=101g(R,(z)) of this backscatter signal is 
depicted in Figure 3.2 (b). The background noise term R, defines the sensi- 
tivity limit of the OCT system. For real OCT systems, this background noise is 
sometimes dominated by relative intensity noise (RIN) of the swept source, 
which may lead to a depth-dependent noise floor. In contrast to that, the noise 
floor in our system originates from thermal noise of the receiver electronics, 
which exhibits a white power spectrum and is hence constant over the depth 
range of interest. Note that all measurements for particle sizing are taken un- 
der oblique incidence of the OCT beam on the sample surface. This avoids 
occurrence of isolated reflection peaks at the sample surface such that the sig- 
nal model according to Eq.(3.1) can be directly used to fit the 
measurement data. 


In real OCT measurements, the measured backscatter depends on further pa- 
rameters, which need to be determined in a calibration step. This comprises 
the divergence of the measurement beam, the decay of the coherence function 
of the swept laser, and the absorption of the matrix material which surrounds 
the scatterers. These influences are approximated by including two calibration 
factors [80,91] q and Q in the single-scattering model according to Eq. (3.1), 


R (z)=Q 6ze DR, 62) 
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Both calibration factors are determined by comparing measured backscatter of 
a NIST-traceable polystyrene particle standard (246 nm diameter polystyrene 
nanospheres in 0.5 wt.% aqueous dispersion, BS-Partikel GmbH, Wiesbaden, 
Germany) with the model calculations according to Eq. (3.1). 


In the special case of spherical scatterers with a size in the order of the wave- 
length, the scattering cross section o, can be modelled by means of Mie’s the- 
ory [50]. As an example, Figure 3.2 (c) shows direction-dependent scattering 
lobes of water-dispersed polystyrene nanospheres with diameters 143 nm and 
506 nm. The larger sphere (506 nm diameter) shows stronger total scattering 
(red), but less pronounced fractional backscatter as compared to the small 
sphere (143 nm diameter). For a small aperture of the scan head (1° half-angle, 
corresponding to a theoretical and experimental lateral resolution of 28 um 
and 36 um, respectively), the Mie scattering lobes can be parameterized using 
the total scattering cross section o, and the backscattering cross section oj. 
Figure 3.2 (d) shows both cross sections o, and œ as a function of the sphere 
diameter d for polystyrene (PS) spheres (refractive index nps = 1.57) dispersed 
in water (refractive index ny, o =1.33) at a wavelength of 1315 nm. In the lim- 
it of small diameters d, Mie scattering can be approximated by Rayleigh scat- 
tering. In this regime, both scattering cross sections increase in proportion to 
d°. For larger diameters, the spheres show resonances, which lead to dips in 
the backscattering cross section œ, whereas the total scattering cross section 
remains unaffected. The relationship between backscattering cross section oj 
and particle size is unambiguous only if the particle diameter is smaller than 
460 nm, which corresponds to roughly half the material wavelength in the 
polystyrene spheres. 


In real measurements, the total scattering cross section o, and the backscatter- 
ing cross section œ cannot be assessed directly. Instead, only the extinction 
coefficient 44 = No, and the backscatter coefficient 44 = No, are extracted 
from the backscatter signal R,(z). For non-absorbing particles in the Rayleigh 
scattering regime, particle size and concentration cannot be separately evalu- 
ated, since both o, and oj increase with d f Figure 3.2 (d). As an example, for 
particle sizes of less than 150 nm (4/9), we can safely assume Rayleigh 
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scattering, and the backscattering probability p, = © / 0, stays constant within 
10 %. As a consequence, a low concentration of bigger particles cannot be dis- 
tinguished from a high concentration of smaller particles. However, for the 
case of nanomaterial characterization, the total volume of nanoparticles added 
to the sample is usually known. 


90° 


Scattering medium 
— 270° 


Scattering cross section [nm] 


10 100 460 1000 

Sphere diameter d [nm] 
Figure 3.2: Concept of model-based dispersion state analysis and sizing of nanoparticles. 
(a) Single scattering model: Incident light is scattered at maximum once inside the medium. A 
light beam with power Pin enters the sample. At each particle (grey), a first portion (or, blue) of 
the incident light in the respective depth z is scattered back into the numerical aperture of the 
optical system, a second portion (æ, red) is scattered into all other directions or is absorbed, and 
the remaining third portion (black) propagates deeper into the medium. (b) Schematic profile of 
the logarithmic backscattering factor Rag (z) =10 lg(R, (z)). (c) Direction-dependent scatter- 
ing lobes according to Mie’s theory, plotted here for the example of polystyrene (PS) nano- 
spheres (diameters 506 nm and 143 nm, refractive index npg =1.57 at 1315 nm) in aqueous 
dispersion (Ayo =1.33). œ: total scattering (red); œ: backward scattering (blue). (d) Total 
scattering œ and backscattering oj cross sections for non-absorbing polystyrene spheres in 
aqueous dispersion. 
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For an increasing degree of dispersion, the average size d of the agglomerates 
decreases and their volume number density N increases in proportion to d”. 
Together with the d°-dependence of the scattering cross sections in the Ray- 
leigh regime, this leads to an overall decrease of the scattering coefficients in 
proportion to d °, which allows robust separation of particle size d and volume 
number density N. Note that this applies to non-absorbing particles only. If the 
extinction coefficient 44 is dominated by the contribution of absorption rather 
than scattering, an increased degree of dispersion could lead to an increase of 
4. This is due to the fact that the breaking-up of agglomerates exposes more 
particles from the inner region to the incident light. These particles from the 
inner region were formerly shielded from light by the absorbing outer shell, 
and did therefore not contribute to attenuation. If absorption dominates, the 
effect of increasing extinction with increasing number of separated absorbing 
particles can be exploited for the analysis of the dispersion state. We use this 
approach in the section on dispersion analysis of epoxy-CNT composites. 


3.2.3 Image-based dispersion-state analysis 


Although nanocomposites ideally feature a homogeneous distribution of na- 
nosized fillers, in practice, microscale agglomerates cannot be completely 
avoided. The size of the agglomerates could reach several hundreds of micro- 
metres, especially at the beginning of the so-called compounding process, 
which usually exploits shear forces to break up particle agglomerates into their 
nanoscale constituents. Therefore, a dispersion-state analysis suitable for pro- 
cess monitoring should be able to cope with nanometre to micrometre sized 
objects. This section is dedicated to imaging and analysis of agglomerates in 
the micrometre range. 


Optical coherence tomography can be used to visualize agglomerates, if their 
dimensions exceed the spatial resolution dx, dy, 6z of the imaging system in x, 
y and z-direction (about 10 um), and if the backscattering is stronger than the 
background noise. Below these limits, an analysis based on a scattering-model 
has to be applied. To identify agglomerate regions with stronger backscatter 
within the OCT image, we use an image segmentation algorithm based on 
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seeded region growing [92]. For quantitative information related to agglomer- 
ate size and number, two independent parameters are extracted from the seg- 
mented images, namely the area fraction (AF) and the perimeter-to-area ratio 
(PAR) of the agglomerates, see Figure 3.3. The area fraction relates the image 
area covered by all identified agglomerates with individual area A; to the entire 
imaging cross section Ato, 


ar= ch, (3.3) 
A 


tot 


For a given volume fraction of the nanosize filler, and assuming spherical ag- 
glomerate shapes which are accurately detected by an ideal measurement sys- 
tem, the AF would be independent of the agglomerate diameter. This can be 
understood by the following consideration: Assume that each spherical ag- 
glomerate splits and decreases in radius by a factor of v. The total amount of 
material remains constant, therefore the number of (smaller) agglomerates in 
the volume increases by v’, whereas the number of agglomerates in the meas- 
urement plane increases by v’. The average area A; of each individual cross 
section in the measurement plane decreases in proportion to v”. Therefore the 
radius change of the agglomerates has no effect on AF. This would render AF 
as non-indicative for agglomerate analysis. For real measurement systems, 
however, a decreasing average agglomerate size will increase the number of 
agglomerates that are smaller than the detection threshold of the image-based 
analysis technique. In this case, a decrease of AF is observed which is corre- 
lated with a decreasing agglomerate size. In our experiments, the detection 
threshold is set to three times the standard deviation of the background noise 
floor in the image. While the threshold influences the measured percentage of 
the area fraction, our choice of the threshold level suffices to judge the disper- 
sion of a certain sample type. 


As a further parameter, the perimeter-to-area ratio (PAR) relates the sum of all 
agglomerate perimeters s; to the sum of all agglomerate areas A,, 


PAR = DS (3.4) 


DA 
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Assuming agglomerates of equally distributed shapes, the PAR depends only 
on the average size of the agglomerates and is otherwise independent of the 
number density. A change in the average shape of all agglomerates during the 
dispersion process would influence the PAR, but this is not to be expected for 
typical dispersion processes that rely on milling or shear-strain-induced break- 
ing of agglomerates. If the dimensions of each agglomerate decrease by a fac- 
tor of v, the perimeter decreases by v, while the cross sectional area decreases 
by v’. Accordingly, the perimeter-to-area ratio of each agglomerate increases 
by a factor of v, and the same applies to the overall PAR, independently of the 
filler content. Image-based analysis is applied to OCT measurements of nano- 
composites both offline and during production. 


Figure 3.3: Illustration of agglomerate areas and perimeters in a cross-sectional OCT image. 
The total image cross section amounts to Ajo; the quantities A,...A, denote the areas (grey) and 
the quantities s,...s, (red) denote the perimeters of the agglomerates. The area fraction (AF) is a 
measure for total agglomerate content and relates the agglomerate areas A,...A, to the total im- 
age cross section Ajo. The perimeter-to-area ratio (PAR) is a measure for the size of the found 
agglomerates and relates perimeters s)...s, to the areas Ay...A,. 


3.3 Results and discussion 


To prove the viability of OCT as a tool for characterization of nanoparticles 
and nanocomposite materials, we performed a series of experiments compris- 
ing accurate nanoparticle sizing as well as nanocomposite dispersion state 
analysis in liquid and solid materials. Our characterization employs image 
evaluation and a simplified model-based scattering analysis. We show that 
OCT methods are useful for measuring nanoscopic particle sizes as well as for 
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analyzing the dispersion of agglomerates in the micro- and millimetre regime. 
The technique can even cope with highly absorbing CNT-loaded materials and 
is well suited for in-line process control. 


3.3.1 Model-based sizing of nanoscale particles 


Polystyrene (PS) nanospheres (refractive index nps = 1.57) in aqueous disper- 
sion (refractive index ny o =1.33) are used as a model system to prove the vi- 
ability of OCT as a tool for model-based nanoparticle sizing. This model sys- 
tem offers higher refractive index contrast than the sample systems 
investigated in following sections (polyamide/clay: npa = 1.59, nejay = 1.54; 
polypropylene/clay: npp = 1.49, nejay = 1.54). All these samples offer sufficient 
backscattering levels. In this experiment, we characterize dispersed NIST- 
traceable polystyrene particles (BS-Partikel GmbH, Wiesbaden, Germany) 
with diameters of 143 nm and 246 nm. The results are depicted in Figure 3.4. 
The particle size is determined by fitting the calibrated single scattering model 
according to Eq. (3.2) to measured OCT depth scans as described in the previ- 
ous section. We use three different samples S1, S2, and S3 with different di- 
ameters and particle concentrations expressed in weight-% (wt. %) of PS par- 
ticles in the dispersion, see Figure 3.4 (d) and (e) for scanning electron images 
of the dried particles. The nominal size and concentration of the investigated 
dispersed particles amount to 143 nm at 2.0 wt. % for S1, 143 nm at 0.5 wt. % 
for S2, and 246 nm at 0.2 wt. % for S3, respectively. Figure 3.4 (a) depicts 
averaged OCT depth scans (light coloured circles) along with the curves of the 
fitted single-scattering model (bright colours). For the fit, measurement data 
have been taken into account only up to a depth of 1.8 mm (225 measurement 
points) to avoid inaccuracies by multiple scattering [93]. The extinction coef- 
ficient 44 and the backscattering coefficient 44, can be extracted from the fit: 
Sample S1 (143 nm, 0.5 wt.%): La = (85 + 23) m“, 
Loi = (97.0 + 7.4)x10 m; Sample S2 (143 nm, 2.0 wt.%): 
Mo = (352426) m', phy = (53.6 + 5.8)x10° m"; Sample S3 (246 nm, 
0.2 wt.%): 4s = (1257.5) m', 13 = (18.0 +2.0)x10° m". The error 
bounds refer to the standard deviation of the averages. Since the particle con- 
centrations and hence the number densities N are known, we can translate 
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these coefficients directly to the scattering cross sections o, and ©, assuming 
that absorption of the PS particles can be neglected at the measurement wave- 
length around 1315 nm. The following total scattering o, and backscattering 
cross sections o, have been extracted: Sample S1 (143 nm, 0.5 wt.%): 
Os = (25.9+£7.0) nm’, oj) = (29.7 £2.3)x10 ‘nm’; Sample S2 (143 nm, 
2.0 wt.%): 0 = (26.9+2.0)nm’, oy =(41.0+4.4)x10 nm; Sample S3 
(246 nm, 0.2 wt.%): og = (488 +29) nm’, o = (70.1 #7.7)x10 ” nm’. The 
values given include the standard deviation of the averages. The cross sections 
o, and œ can be related to the diameters of the respective particles using 
Mie’s theory [50], see Figure 3.4 (b). The total scattering cross section o; (left 
axis, red) and the backscattering cross section œ (right axis, blue) of polysty- 
rene nanospheres dispersed in water are calculated and plotted as a function of 
the sphere diameter, assuming incident light of 1315 nm wavelength and a 
system aperture of 0.018 (1° half-angle). 
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Figure 3.4: Size determination of sub-wavelength nanoparticles with OCT. (a) Measured OCT 
depth scans of aqueous dispersions of polystyrene (PS) nanospheres (S1: 143 nm diameter, 
0.5 wt.%, grey circles: measurement, black line: fit; S2: 143 nm diameter, 2.0 wt.%, light green 
circles: measurement, green line: fit; S3: 246 nm diameter, 0.2 wt.%, light orange circles: 
measurement, orange line: fit), fitted with the calibrated single-scattering model according to 
Eq. (3.2). (b) Total scattering (red) and backscattering (blue) cross sections, according to Mie’s 
theory for polystyrene spheres in water. Both scattering cross sections are compared with the 
measurements, depicted only for the samples S1 (oj): total scattering cross section, Obı 
backscattering cross section) and S3 (0,3, Op, respectively). Both parameters can be attributed 
to sphere diameters, horizontal axis. (c) Comparison of nominal particle sizes with measured 
sizes. Relative (absolute) deviations from the nominal values are maximal 4% (5 nm). 
(d,e) Scanning electron micrograph (SEM) of the nanospheres having diameters of 143 nm and 
246, respectively. 
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From the measurement fits in Figure 3.4 (a), the scattering parameters 0,1, O41 
for sample S1 and 0,3, o43 for sample S3 are extracted. The sphere diameter 
can then be read from the Mie calculations displayed in Figure 3.4 (b), vertical 
arrows. For clarity, the analysis is only depicted for sample S1, but sample S2 
yields very similar results. 


The table in Figure 3.4 (c) summarizes the results as derived from scattering 
cross section o, and backscattering cross section ©, together with the associ- 
ated standard deviations of the averages, and compares them with the nominal 
particle sizes and their standard deviations. The standard deviation measured 
by OCT can be lower than the standard deviation of the particle size distribu- 
tion, since various particles contribute to one OCT measurement. Measure- 
ment errors may arise from refractive index uncertainties, where an uncertain- 
ty of 0.01 would lead to a size determination error of 2 nm to 4 nm. Relative 
deviations between nominally and measured particle sizes are below 4 %. 


3.3.2 Model-based nanoscale dispersion-state analysis 


Owing to the high sensitivity of the OCT technique, even strongly absorbing 
nanomaterials such as carbon-nanotube (CNT) composites can be analysed. 
This is a key feature, as CNT-loaded composites represent an important mar- 
ket segment [94]. Standard light scattering techniques like DLS or SLS are not 
applicable because the scattered or transmitted optical powers are small. Con- 
ventional CNT dispersion analysis relies either on thin material layers that are 
investigated in a light microscope (LM), or on the determination of macro- 
scopic material parameters like the dielectric permittivity that can also indicate 
the dispersion quality [95]. For our experiments, we use a multi-wall CNT 
(MWCNT) dispersion in an epoxy resin (Araldite LY 556, Huntsman Ad- 
vanced Materials GmbH, Basel, Switzerland). The samples contain 0.12 wt. % 
of MWCNT (NC7000, Nanocyl S.A., Sambreville, Belgium), and were pre- 
pared by 1, 3 and 5 milling cycles in a three-roll mill. We compare the results 
of OCT-based backscattering and extinction parameter analysis with the re- 
sults of standard LM analysis as well as with rheological and dielectric meas- 
urements. For OCT measurements, the CNT-filled resin samples were filled 
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into cuvettes and heated up to a temperature of 50°C for melting the resin, 
thereby avoiding scattering from resin crystals. 


The results of the dispersion analysis are shown in Figure 3.5. Conventional 
dispersion analysis of thin composite layers by light microscopy, see Fig- 
ure 3.5 (a), reveals significant agglomerates after the first milling cycle. After 
three milling cycles, the size of agglomerates reduces. After five milling cy- 
cles, the light microscope images could lead to the conclusion that dispersion 
quality seems to decrease again, but other measurement techniques, like OCT 
analysis and rheological characterization, contradict this finding. Light mi- 
croscopy for dispersion analysis, while being a state-of-the-art technique, suf- 
fers from small imaging areas. This could lead to an erroneous interpretation, 
if the agglomerate concentration is not spatially homogeneous over the whole 
sample. The same samples were investigated using OCT, see Figure 3.5 (b), 
which shows extinction coefficients 44 and the backscattering coefficients 4% 
in a two-dimensional plot. For each of the samples, we take ten measurements, 
each consisting of 5000 depth-scans which were taken while laterally moving 
the sample over a distance of 2 mm. The depth scans are averaged, and the 
extinction coefficient 44 and the backscattering coefficient 4h are extracted by 
fitting Eq. (3.2) to the measurement data. This has been repeated at ten differ- 
ent regions of each sample, each region corresponding to a cross in Fig- 
ure 3.5 (b). Note that the particle sizes cannot be extracted from these data: In 
contrast to the situation for pure nanosphere dispersions used as a reference, 
the CNT agglomerates exhibit a large variety of shapes, and a Mie scattering 
theory based on simple spheres cannot be applied. Yet, the analysis of the ex- 
tinction and backscattering coefficients suffices for a qualitative dispersion 
analysis. The horizontal and vertical error bars in Figure 3.5 (b) represent the 
average (solid dots) and the standard deviation of the measurements obtained 
for each sample. Although it is not possible to exactly discriminate between 
scattering and absorption, we may assume that the high extinction values of 
5 mm up to 15mm! arise mainly from absorption of the CNT and only 
marginally from scattering. This is supported by strong reported absorp- 
tion [96] of 1000 mm ' due to separated CNT in composites with comparable 
concentrations. Consequently, we did not take multi-scattering into account. A 
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clear tendency with respect to the number of milling cycles can be seen: The 
more agglomerates are broken, the more isolated CNT can contribute to ab- 
sorption and hence to the extinction 44, which rises significantly between one 
and three milling cycles, and increases slightly for five cycles. Simultaneous- 
ly, the backscattering coefficient 44 increases with the number of milling cy- 
cles and shows the same tendency as the extinction 44. The backscattering co- 
efficient 44 = œN, however, does not necessarily increase if the dispersion 
quality increases so that the agglomerates become smaller and their number N 
becomes larger. It is impossible to state in general how smaller agglomerates 
contribute to the backscattering inside the light receiving aperture, because op 
strongly depends on the angle dependency of the backscattering for the respec- 
tive agglomerate shape. In our specific case, 44 increases with an increasing 
degree of dispersion, and therefore allows a qualitative judgement of the state 
of the sample dispersion. 


These results of the OCT analysis are confirmed by a rheological characteriza- 
tion of the relative permittivity &. For this purpose, the dispersions were in- 
vestigated using a rheometer (MCR 501, Anton Paar GmbH, Graz, Austria) in 
combination with a programmable LCR bridge for measuring inductance, ca- 
pacitance and resistance (HM 8118, HAMEG Instruments GmbH, Main- 
hausen, Germany). The rheometer consists of two parallel rotatable plates at a 
separation of s = 1 mm. The sample dispersion fills the gap between the plates, 
which are isolated electrically from each other and which are connected to the 
LCR bridge for measuring the capacitance. During a measurement, one of the 
plates performs an oscillating rotation. The amplitude b, of the oscillation at 
the outer radius of the plate, related to the plate separation s results in the 
strain amplitude y =b,/s. Both the measurement setup and the measurement 
procedure are described in detail in Ref. [95]. The relative permittivity is 
measured at a frequency of 1 kHz and increases monotonically from the pure 
epoxy material over | to 3 to 5 milling cycles: This indicates a continuous im- 
provement of dispersion quality, see Figure 3.5 (c), since an increasing num- 
ber of CNT are separated from agglomerates and contribute to the electric po- 
larization of the nanocomposite. Note that the rheological measurements 
themselves introduce shear strain into the samples, which could lead to further 
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exfoliation of CNT from the agglomerates and hence to an increase of the rela- 
tive permittivity during characterization of the samples after three milling cy- 
cles. This additional exfoliation might be not possible for the mixture after one 
milling cycle (insufficient exfoliation) and five milling cycles (largely com- 
plete exfoliation). At higher strain amplitudes (not shown), the separated CNT 
align in parallel to the rheometer plates and thus perpendicularly to the electric 
field. In this case the CNT cease to influence the permittivity. 
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Figure 3.5: OCT characterization of highly absorbing epoxy/carbon nanotube (CNT) compo- 
site materials. Comparison to conventional characterization techniques and macroscopic mate- 
rial properties. The composite was milled over 1, 3 and 5 cycles in a three-roll-mill to improve 
dispersion. (a) Conventional dispersion state analysis by light microscopy of thin composite 
layers. Black areas are attributed to agglomerates. (b) OCT-determined backscattering 4} and 
extinction 44 coefficients, measured in 10 different spatial regions (crosses) of each sample (av- 
erage: solid dots; standard deviations: horizontal and vertical error bars). (c) Rheologic dielec- 
tric characterization of the relative permittivity & as a function of applied shear strain after 1, 3, 
and 5 milling cycles. The relative permittivity increases monotonically from the pure epoxy 
material over 1 to 3 to 5 milling cycles, indicating an increase in the separation of the CNT. 
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3.3 Results and discussion 


In contrast to the microscopy analysis and the permittivity measurement, OCT 
characterizes the heated sample without any further preparation. The results 
correlate well with the permittivity measurement and with the number of dis- 
persion cycles. In contrast to light microscopy, OCT also reveals the increase 
in dispersion quality between the third and the fifth dispersion cycle. This ex- 
periment demonstrates that OCT metrology is useful to characterize even 
highly absorbing nanomaterials, outperforming by far the elaborate and time 
consuming conventional light microscopy method. 


3.3.3 Image-based dispersion-state analysis for 
microscale agglomerates 


For experimentally assessing the ability of the OCT technique to analyse 
nanocomposites with large agglomerates, composites of polyamide (PA; Ba- 
damid B70, Bada AG, Bühl, Germany) and nanoclay particles with 5 wt.% 
concentration were prepared in a compounding extruder [97]. For controlling 
the dispersion quality during fabrication, one batch of samples was prepared 
from a clay with modified surface (1.34, Nanocor Inc., Hoffman Estates, IL, 
USA), the other one was prepared from unmodified clay (PGN, Nanocor Inc.). 
The surface modification is expected to lead to an improved dispersion state 
and to decreased sizes of agglomerates as compared to the unmodified clay. 
After compounding, all samples are granulated without further sample pro- 
cessing. Cross-sectional images (B-scans) are taken from the granules with our 
laboratory OCT system. Representative B-scans from the PA/nanoclay com- 
posite without and with surface modification are shown in Figure 3.6 (a) and 
(b), where the pixel brightness indicates measured backscattering. The sample 
with no modification, Figure 3.6 (a) shows large lengthy bright areas, indicat- 
ing large clay agglomerates that extend over several hundreds of micrometres, 
whereas the sample with modification, Figure 3.6 (b), features small bright 
spots only from which small agglomerates can be inferred. The image seg- 
mentation algorithm detects bright regions automatically; the corresponding 
borders are drawn as red lines. 
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For each sample type, three B-scans are taken, segmented, and the area frac- 
tion (AF) as well as the perimeter-to-area ratio (PAR) are calculated for each 
measurement. Figure 3.6 (c) and (d) display the average AF and the average 
PAR for both sample types along with the corresponding error bars showing 
the standard deviation of the averages. A significant decrease in the AF is ob- 
tained when the surface of the clay filler is modified (AF oa = 2.0 % + 0.4 %) 
as compared to the filler without surface modifications (AFyny. 
mod = 5.3 % + 0.8 %), Figure 3.6 (c). This indicates that a surface modification 
significantly decreases the agglomerate size such that a substantial fraction of 
agglomerates becomes smaller than the size detection threshold of the image 
analysis. This is in accordance with the observation that the cross-sectional 
image of the surface-modified clay composite exhibits a large number of 
bright spots, each of which contains only a few pixels, see Figure 3.6 (b). At 
the same time, the PAR increases significantly with modification of the filler 
surface (PARgamoag = 132 mm | 8.2 mm |, PAR „oa = 178 mm! +4.1 mm ') 
which confirms the decrease in average agglomerate size, Figure 3.6 (d). 


a No surface modification 


a 
S 


Area fraction AF [%] 
Perimeter-to-area ratio (PAR) [mm] 


! No surface mod. Surface mod. No surface mod. Surface mod. 


Figure 3.6: Image-based dispersion-state analysis for microscale agglomerates, showing the 
impact of chemical surface modification on dispersibility. (a) Cross-sectional image (B-Scan) 
of a PA/nanoclay composite with 5 wt-% nanoclay content without surface modification. Bright 
spots (strong backscattering) with red borders indicate agglomerates. Right part: magnified 
section. (b) Cross-sectional image of a PA/nanoclay-composite with surface-modified clay par- 
ticles having the same concentration as the sample in (a). Right part: magnified section. 
(c) Area fraction (AF) covered by detected agglomerates. The circles denote average values of 
the three samples, and the error bars indicate the standard deviation of the average. 
(d) Perimeter-to-area ratio (PAR) of the detected agglomerates. The circles denote average val- 
ues of the three samples, and the error bars indicate the standard deviation of the average. 
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3.3 Results and discussion 


These experiments prove the viability of image-based OCT analysis to charac- 
terize the dispersion state of samples with relatively large agglomerates sized 
from a few micrometres to hundreds of micrometres. The quantitative evalua- 
tion of further dispersion-related parameters like agglomerate shape and num- 
ber would further increase the reliability and robustness of our technique. 


3.3.4 Demonstration of in-line dispersion-state analysis 


Nanocomposite development is hampered by rather long development cycles, 
which are dominated by time consuming off-line characterization. In this sec- 
tion, OCT is demonstrated to be well-suited for continuous dispersion-state 
monitoring in a production environment. 


In the following experiment, a production-scale twin-screw compounding ex- 
truder (Leistritz GmbH, Nuremberg, Germany) is used for production of a 
polypropylene(PP)/nanoclay composite with a mass throughput of 6 kg/h. The 
applied host material is PP (R352-08R, Dow Chemical, Midland, MI, USA) 
and the nanosized filler is a nanoclay (Cloisite 15A, Southern Clay Products, 
Gonzales, TX, USA). The screws introduce shear into the polymer melt filled 
with nanoparticles and thus disperse the particles [97]. With this extruder, 
nanocomposites were dispersed with different revolution speeds. In general, 
increasing speed comes with higher energy input into the material and causes 
a better nanoparticle dispersion [21]. In order to characterize the dispersion 
state of the PP/clay melt during production, the extruder has been equipped 
with the OCT probe for in-line process monitoring, see Figure 3.1 (c). The 
probe is mounted to a sensor port of the extruder located at the end of the ma- 
chine close to the die. The optical window of the probe is in direct contact 
with the main stream of the extruded nanoparticle loaded polymer melt. With 
the OCT probe, A-scans of the medium underneath the probe window are 
measured continuously. The flow of the melt causes a continuous movement 
of the material seen by the OCT measurement beam, thereby replacing a lat- 
eral movement of the measurement beam in conventional measurements. This 
results in a temporal change of the backscatter signal, which can be interpreted 
as 2D image data, where one dimension corresponds to time and the other 
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dimension is the usual imaging depth measured from the probe window, see 
Figure 3.7. 


Figure 3.7 (a), (b) and (c) show the in-line OCT data for the screw revolution 
speeds 200 rpm, 500 rpm and 800 rpm, respectively. Note that, due to non- 
uniform mass flow as a function of depth, it not possible to directly map the 
measurement time (horizontal axis) to spatial coordinates of the sample. 
Bright spots indicate strong backscatter and are attributed to agglomerates. 
The straight horizontal line at z = 0 originates from permanent reflections at 
the probe window. With higher screw speed and accordingly higher shear 
strain inside the material, fewer agglomerates are visible indicating an im- 
provement in dispersion quality. At the highest screw speed of 800 rpm, al- 
most no agglomerates are visible. This impression is also confirmed by quanti- 
tative dispersion-related parameters. Since the measurement is performed as a 
function of time rather than as a function of lateral position, a perimeter-to- 
area ratio in the strict sense as defined by Eq. (3.4) cannot be calculated. 


a 200 rpm b 500 rpm c 800 rpm iB d 
60 


Depth z [um] 
Modified perimeter-to-area ratio PAR [%] 
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Figure 3.7: Demonstration of in-line dispersion-state analysis during operation of a compound- 
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ing extruder. The OCT measurement beam is fixed, the vertical axis of the OCT images indi- 
cates the imaging depth and the horizontal axis corresponds to the measurement time. Bright 
spots: strong backscatter of agglomerates. Straight horizontal line at z ~ 0: permanent reflec- 
tions at the probe window. Material system: polypropylene/nanoclay. Mass flow: 6 kg/h. Clay 
content: 1 %. Bottom: magnification after segmentation. (a) OCT scan at 200 rpm. (b) OCT 
scan at 500 rpm. (c) OCT scan at 800 rpm. (d) Quantitative analysis of the dispersion state by 


area fraction (AF) and modified perimeter-to-area ratio PAR from five OCT scans. The circles 
denote average values of the five scans, and the error bars indicate the standard deviation of the 
averages. 
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3.4 Conclusions 


Instead, we use a modified quantity PAR , which essentially corresponds to 
the PAR except that both, the agglomerate perimeters s; and the agglomerate 
areas A;, are expressed by pixel numbers rather than by physical lengths and 
areas. Note that PAR is a dimensionless quantity in contrast to PAR, which 
has the unit m '. The results are depicted in Figure 3.7 (d), where the circles 
denote average values from five OCT scans, and the error bars denote the re- 
spective standard deviations of the averages. With increasing speed, the area 
fraction AF reduces, and the PAR rises. Both quantities indicate that the ag- 
glomerates inside the material flow become smaller corresponding to a better 
dispersion quality of the composite. These results reveal a clear relationship 
between operation parameters of the machine and OCT-measured dispersion 
parameters. This is the first demonstration of an in-line dispersion-state analy- 
sis in a nanocomposite production line. The results encourage the application 
of our technique to more material systems for controlling multiple 
production parameters. 


3.4 Conclusions 


In this work, we present a novel and a highly attractive approach to nano- 
material analysis applied to nanoparticles and nanocomposites. The method 
uses optical coherence tomography (OCT) and model-based parameter extrac- 
tion. Our approach enables detection and quantification of nanoparticles and 
agglomerates over a wide range of size scales: Image segmentation of OCT 
data sets is well suited for dispersion-state characterization of nanocomposites 
with agglomerates in the micrometre range, whereas model-based scattering 
analysis lends itself to size determination of nanoparticles below the resolution 
limit. We elaborate the measurement technique along with theoretical models 
and demonstrate the viability of the procedure in a series of proof-of-principle 
experiments. A wide variety of material systems is investigated by our exper- 
iments: A first demonstration shows accurate OCT-based nanoparticle sizing. 
An in-situ dispersion-state analysis characterizes strongly absorbing CNT 
composites in liquid media. Finally, we perform in-line monitoring of the 
compounding process in a state-of-the-art production line. Major challenges in 
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industrial applications are the stringent requirements with respect to mechani- 
cal stability and size of the OCT system. These requirements can be met by 
integration of interferometer and detection system on a silicon photonic chip 
[98]. We believe that OCT has the potential to fill a metrology gap in the 
emerging field of nanocomposite technology. We conclude that OCT metrolo- 
gy opens new directions in material analysis, both in laboratory and 
production environments. 


[End of Publication P2] 
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4 Nanoparticle size and shape 
characterization using PS-OCT 


The following chapter is a manuscript intended for publication in a peer- 
reviewed journal. Slight adaptions have been made to match the format and 
the structure of this thesis. 


Nanoparticles of diverse shapes play an increasingly important role in medical 
and technical applications. Their tiny size comes along with unique physical 
properties, which are tightly coupled to their actual size and shape. In conse- 
quence, continuous quality control in nanoparticle fabrication is mandatory. 
Transmission or scanning electron-microscopy is a costly and time-consuming 
method, yielding results from small sample batches only. Here we show that 
polarization-sensitive optical coherence tomography allows to determine the 
size and the shape of nanoparticles in a quick process with small technical ef- 
fort and with no need for laborious sample preparation. The technique relies 
on measuring the depth and polarization-dependent backscatter and on evalu- 
ating the results by a semi-analytical model. We found the backscatter meas- 
urement to be sufficiently sensitive on gold nanorod geometry to enable diam- 
eter determination with 6.2 nm accuracy and maximum deviation in aspect 
ratio measurement of 10 %. These experiments correspond to the first demon- 
stration of quantitative nanoparticle shape determination with optical coher- 
ence tomography. 


4.1 Introduction 


Nanoparticles are attractive for both scientific examinations and technical ap- 
plications. Their size is in or below the range of the optical wavelength and is 
thus only one to four orders of magnitude bigger than an atomic diameter, thus 
leading to huge surface-to-volume ratios. As a consequence, in contrast to 
bulk materials, size and shape of the particles govern their physical properties 
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and become the essential parameter for targeted material design. As an exam- 
ple, the shape-dependent surface plasmon resonance frequency of metallic na- 
norods allows to use them as highly-sensitive refractive index sen- 
sors [99,100]. Moreover, gold and platinum nanoparticles can be used as 
chemical catalysts for which the reaction rate depends significantly on size 
and shape [14,101], and also the cellular uptake of gold nanorods into biologi- 
cal tissue, when used as imaging markers, turns out to be shape- and size- 
dependent [9]. The same is true for soft-matter nanorods as nanocarriers in 
tumour therapy [102]. Nanoclays form layered structures and thus feature a 
large lateral extension and small height. This special shape accounts for their 
outstanding barrier properties [74], when incorporated into a polymer host ma- 
terial. The aspect ratio of carbon nanotubes (CNT) determines their ability to 
build up electrically conductive paths inside isolating matrix materials. Trans- 
parent polymers can thus become conductive and might offer an alternative to 
expensive indium tin oxide (ITO) [12]. All these applications have in common 
that size and shape of the nanoparticles are crucial to the macroscopic material 
properties and thus to the final performance of the product. In this context, 
characterization of the nanostructure is of highest relevance, either indirectly 
by characterization of the final macroscopic properties or directly using imag- 
ing, scattering or filtration techniques. 


State-of-the-art size and shape characterization of nanoparticles mainly relies 
on microscopic imaging. This requires high-resolution techniques such as 
scanning-electron microscopy (SEM) [16], transmission electron microscopy 
(TEM) [17] or atomic force microscopy (AFM) [103,104]. In all cases, sample 
preparation is laborious and costly and might create artefacts, e.g. agglomera- 
tion effects during the drying process in sample preparation. Moreover, SEM, 
TEM and AFM are limited to small sample volumes, which are not necessarily 
representative of the entire batch. Light scattering techniques, in contrast to 
imaging methods, are viable tools for determining particle size distributions 
from large sample volumes. Static light scattering (SLS) uses angle-resolved 
and/or spectrally resolved detection of scattered light [22,23] and needs a ra- 
ther complex optical setup, when large scattering angles are to be taken into 
account. Dynamic light scattering (DLS) methods make use of temporal 
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fluctuations of interference patterns of scattered light to calculate the hydro- 
dynamic diameter of particles within the respective solvent [75]. Relying on 
the observation of the Brownian motion of single particles, the technique is 
limited to rather low concentrations. Both SLS and DLS preferentially provide 
information on particle size, while particle shape cannot be determined [24]. 
X-ray diffraction (XRD) [77] or X-ray scattering (SAXS) [78] can be used to 
reveal the atomic or molecular structure inside the material. However, expen- 
sive instrumentation and a small probing region limit the application range of 
these methods to small batch sizes. Tunable resistive pulse shaping 
(TRPS) [25] measures the particle size by observing the ionic current through 
a small pore. Once a particle passes this pore the ionic flow is affected for a 
certain time. After calibration, this technique allows accurate determination of 
particle size distributions. Information on particle shape however cannot be 
extracted. Particle size distributions can also be measured using centrifugal 
sedimentation techniques, where particles with different dimension are sepa- 
rated by differences in their sedimentation speed in presence of centrifugal 
forces [25]. While this technique is applicable to a broad size range, high ex- 
perimental effort is necessary and no information on particle shape is ob- 
tained. 


Optical coherence tomography (OCT) represents a combination of light scat- 
tering techniques with optical imaging is [1,28,79]. OCT provides three- 
dimensional imaging data from the bulk of a sample and avoids expensive 
sample preparation. Previous demonstrations of nanoparticle analysis using 
OCT have shown size determination with nanometre accuracy [10] and nano- 
composite dispersion-state characterization [10]. Microparticles have been in- 
vestigated using image-based analysis [81]. OCT has been applied to particle 
diffusion measurements and showed particle size-dependent temporal varia- 
tions [82,83]. Angle-resolved OCT has proven to be able to measure particle 
sizes down to 5 um [84], by evaluating the angular-dependent scattering lobes 
with Mie’s theory. Gold nanorods have been differentiated from collagen 
cells, exploiting their strong polarization-cross talk in OCT measure- 
ments [105]. While micro- and nanometre particle size determination has been 
addressed by some OCT-based methods, the exact determination of 
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anisotropic particle shapes, e.g., particle sizes in long and short dimension, 
still needs to be shown. 


In this section we demonstrate that polarization-sensitive OCT (PS-OCT) [59] 
allows accurate determination of both nanoparticle size and shape. By evalua- 
tion of backscatter signals in different polarizations and by comparison to 
electromagnetic field simulations of differently shaped particles, we demon- 
strate precise determination of length and diameter of gold nanorods in aque- 
ous dispersion. Without the need for sample preparation, this technique allows 
to measure nanometre particle dimensions in statistically relevant volumes of 
several cubic millimetres. We determine size and shape of gold nanorods in a 
unimodal aqueous dispersion by comparing measured scattering parameters 
with simulation results. Our samples cover particles which range from 15 to 
50 nm in particle diameter and from 1 to 16 in aspect ratio. In practice, this 
technique can be extended to the analysis of many nanoparticles with prolate 
or oblate shape as long as the refractive index contrast to the surrounding me- 
dia is sufficient and absorption is moderate. 


4.2 | Measurement principle: Size- and shape- 
dependent backscattering 


Light scattering methods allow size and shape measurements beyond the reso- 
lution of optical imaging systems. The analysis presented here links the geom- 
etry of the particles to their scattering properties and relates those to polariza- 
tion-sensitive backscattering measurements with a PS-OCT system. First, we 
introduce the applied scattering model on the basis of nanospheres and their 
size-dependent scattering. Second, we expand the model to non-spherical scat- 
terers using a polarization-dependent scattering model. 


For the analysis, a single-scattering model [89] is used, under the assumption 
that incident light is scattered at maximum once inside the medium, which is 
appropriate for weakly scattering samples and for the analysis of moderate 
depths [90]. Figure 4.1 depicts the scattering model in two exemplary media: 
nanospheres, Figure 4.1 (a), and non-spherical particles Figure 4.1 (c), with 
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their respective polarization-dependent backscattering curves Figure 4.1 (b) 
and Figure 4.1 (d). As shown in the case of nanospheres Figure 4.1 (a), a light 
beam with incident power Pin and polarization Ein enters the sample. At each 
particle (grey), one part of the incident light at depth z is scattered back into 
the numerical aperture of the system, another part is scattered into other direc- 
tions or is absorbed (black), and the remainder propagates deeper into the me- 
dium. Mie’s theory [50] is applicable to the calculation of angle-dependent 
scattering amplitudes of spherical particles with sizes in the range of or small- 
er than the optical wavelength of the incident light. The scattering amplitudes 
are dependent on size of the spheres, refractive index of sphere and surround- 
ing material, as well as on the optical wavelength. Thus, observation of the 
angle-dependent scattering can be used for size analysis of nanospheres. In the 
case of spherical scatterers, the polarization of the output field Eoutm is the 
same as the incident polarization Ein [50,106,107]. This is an indication of 
spherecity of the scatteres, if their size is far below the wavelength and multi- 
ple scattering can be neglected. The total scattering cross section o, and the 
absorption cross section o, of a particle sum up to the extinction cross section 
a. For an ensemble of particles with number density N, the extinction coeffi- 
cient can be expressed as 44 = No. Analogously, the backscattering cross sec- 
tion into the —180° direction is given by ©, and the backscattering coefficient 
is £4 = Noy. The depth-dependent decay dP(z)/dz=-wP(z) of the forward 


travelling power P(z) follows an exponential function P (z) =P e“. Fora 


in 
given depth z, the power is scattered back from inside a depth element &, and 
is captured by the numerical aperture of the system, P,(z)=4,6zP(z)e “, 


where e% 


accounts for the losses in the back propagation path. In an OCT 
measurement the depth-dependent backscattering factor Ro(z) = Pin/Pp(z) + Ra 
is measured, including a system dependent noise floor R,, which defines the 
sensitivity of the system expressed as minimum measurable sample reflectivi- 
ty, see Section 2.1.4. In summary, the scalar single-scattering model, which 


omits polarization effects is given by 


Ro (z) = Moz eM +R, . (4.1) 
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In practice, the measured backscattering signal depends on further effects, 
which need to be determined in a calibration step. This comprises the decay of 
the coherence function of the frequency-swept laser, the divergence of the 
measurement beam, and the absorption of the material, which surrounds the 
scatterers. These influences are taken into account by including two calibra- 
tion factors [10,91] q and Q in the single-scattering model according 
to Eq. (4.1), 


R(z)=0m6z e (4 *9 +R, (4.2) 


Both calibration factors are determined by comparison of measured backscat- 
ter of a NIST-traceable polystyrene particle standard with 246 nm sphere di- 
ameter in 0.5 wt.% aqueous dispersion (BS-Partikel GmbH, Wiesbaden, Ger- 
many) to model calculations according to Eq. (4.1). 


A sketch of a semi-logarithmic plot Rap(z) = 10 1g(R,(z)) of both the polariza- 
tion-maintaining backscattering factor R,(z) (blue) and of the cross- 
polarization backscattering factor R.(z) (red) is shown in Figure 4.1 (b). With 
no polarization conversion, the cross-polarization backscatter is zero, and the 
measured signal is given by the associated noise level. 


Light scattering of particles with non-spherical shape can lead to backscatter 
with altered polarization, see Figure 4.1 (c). In general, light scattering is 
caused by the excitation of dipole oscillations inside the particles and the as- 
sociated radiation of light. The orientation of the dipoles and thus the polariza- 
tion of the emitted light depend on both incident polarization and particle ge- 
ometry [50]. For elongated particles, dipoles are preferentially excited in the 
direction of the long axis of the particle, thus showing a polarization change of 
backscattered with respect to the incoming light. With an ensemble of ran- 
domly oriented particles, all random backscatter polarizations sum up, and an 
average backscatter polarization is obtained, which differs the more from the 
incident polarization, the more the particles’ shape deviate from that of 
a sphere. 
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Figure 4.1: Polarization-sensitive scattering-model-based size and shape analysis of nanoparti- 
cles. (a) Scattering of a linearly polarized laser beam with power Pin inside a dispersion of 
spherical nanoparticles, with aspect ratio AR = 1. Concentration and refractive index contrast 
are low enough to assume single scattering inside the medium. The input polarization Ein is 
considered vertical here (blue colour). Backscattering of spherical particles in the -180° direc- 
tion maintains the polarization, so the backscattering cross section oj, = Ob m equals the polariza- 
tion-maintaining backscattering cross section. The total extinction cross section is given by a. 
Light that is neither scattered nor absorbed propagates deeper into the medium. Poum denotes 
the output power of light in same polarization as of the input light, Ein || Eout,m . (b) A single- 
scattering model leads to a polarization maintaining backscattering factor Rm, describing the 
depth-dependent backscattering factor into the polarization of incidence. In case of spherical 
particles the cross-polarization backscattering factor R. is zero. Together with a system- 
dependent noise level R,, both factors are plotted in dB-scale along the depth axis z. (c) Scatter- 
ing in a medium with non-spherical particles of aspect ratio AR > 1. The input polarization 
(vertical) is only partially maintained, and the backscattering cross section Ob = Ob,m + Ob. iS 
given by the sum of both the polarization-maintaining component op,m and the cross-polarized 
component ©, (red). Accordingly, output light is polarized partially in parallel to the input 
Ein || Eoutm and orthogonal Ein L Eoutc. Poutm and Pou denote the corresponding backscat- 
tered powers. (d) In contrast to the case in (b), the cross-polarized backscattering factor R. is 
non-zero for non-spherical particles. The ratio R./R, can be taken as a measure for non- 
sphericity of the investigated particles. 


Analogously to Figure 4.1 (a), a beam with power Pin and an electric field ori- 
entation Ein is incident on the sample. Interaction of the non-spherical particles 
with the incident light is expressed by the total extinction cross section o;, the 
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polarization-maintaining backscattering cross sections Ob m describing polari- 
zation maintaining backscatter, and the cross-polarization backscattering cross 
section Oj. The backscattered light is hence composed of light with main- 
tained polarization with power Poum and electric-field Eoutm, which is polar- 
ized parallel to the incident field, and of light with power Pou, and electric- 
field Eoutc, which is polarized orthogonally to the incident field. 


The semi-logarithmic depth-dependent backscatter signal of a sample with 
non-spherical particles is depicted in Figure 4.1 (d). In comparison to the sca- 
lar case, the backscattering coefficient Z} = Nō, is vectorial, 


7 Ho,m 
Ly -( v; ) $ (4.3) 
My 
and the backscatter signal is given by 
= R 2 
R, (z) = in (2) = Hom Qôz ea +a)z + R, . (4.4) 
i R: (z) Hoc 


The difference of the two backscattering coefficients, 4m and 4h,, leads to 
offsets between the logarithmic plots R„(z) and R.(z) in Figure 4.1 (d). This 
difference is denoted as polarization cross talk 7= 10 Ig(i%../ Am), which we 
define in dB. The decay of both curves is the same, since the particles in the 
sample are randomly oriented and the extinction is thus dependent on input or 
output polarization. While the backscatter curve in the scalar model contains 
information on particle size and number density, the vectorial polarization- 
dependent model additionally supplies information on particle shape. In the 
following sections, dispersions of particles with different diameters and aspect 
ratios are characterized by measuring their polarization-dependent 
backscattering signals. 


4.3 Measurement results 


To determine the viability of our concept, we analyse aqueous dispersions of 
polystyrene nanospheres and aqueous dispersions of gold nanorods with 
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different aspect ratios. The analysis is performed using a polarization-sensitive 
OCT system, which is described in detail in Section 4.6. Measurements of the 
nanospheres show negligible polarization conversion, which is in accordance 
with theoretical considerations [50,106,107]. For nanorods, we find shape- 
dependent polarization conversion. By comparing the measurement data to 
electromagnetic simulations, we find that diameter and aspect ratio (AR) of 
the rods can be determined with nanometre accuracy. 


In our experiments, aqueous dispersions of gold nanorods with dg = 25 nm 
diameter and average lengths of JG; = 102 nm (sample G1) and loz = 256 nm 
(sample G2; both types of nanorods from Nanopartz Inc., Loveland, Colorado, 
USA) have been characterized with the PS-OCT system. These rods feature 
respective aspect ratios of ARg; = 4 and ARg = 10 and number densities of 
Noi = 5.7x10" m° and Noo = 1.1x10" m’, corresponding to weight concen- 
trations of 0.54 wt. % and 0.026 wt. %, respectively. For comparison, addi- 
tionally a sample (S4) with spherical particles (ARs, = 1), containing an aque- 
ous dispersion of polystyrene nanospheres of diameter of dss = 246 nm, 
number density of Ns4 = 6.4x10'’ m” and weight concentration of 0.5 wt. % is 
investigated. For preparation, all samples were filled into quartz glass cuvettes 
with broad optical transmittance. In Figure 4.2 (a-c), OCT measurement re- 
sults are depicted along with the associated scanning electron microscope 
(SEM) images in Figure 4.2 (d-f) of the dried particles. The measurement re- 
sults consist of two curves, polarization-maintaining backscatter R„(z) and 
cross-polarized backscatter R.(z), both as a function of depth inside the sam- 
ple. Figure 4.2 (a) depicts the results of the microsphere sample (S4), with 
R,(Z) showing an exponential decay inside the sample (z > 0). The data for 
z < 0 is defined by the noise level. By tilted incidence on the sample surface, 
facet reflections could be avoided. The solid curves represent a fit of the scat- 
tering model in Eq. (4.4) to the measured data (light coloured). Applying 
Mie’s theory to R,,(z), the spheres’ sizes can be extracted with nanometre ac- 
curacy, as we have demonstrated previously [10]. The cross-polarized 
backscatter R.(z) is more than 20 dB below R,,(z). Since the polarization ex- 
tinction ratio of the system is around 20 dB, this data can be considered to be 
below the noise level and no scattering model fit is reasonable. The complete 
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polarization maintenance of the backscatter leads to the conclusion that the 
sample particles are completely spherical. The SEM picture, Figure 4.2 (d), of 
the dried particles of sample S4 confirms this finding. The gold nanorod sam- 
ple G1 with AR of 4, see Figure 4.2 (b), shows an exponential decay both for 
the polarization maintaining backscatter R,(z) and the cross-polarized 
backscatter R.(z). 


Sample S4 Sample G1 Sample G2 
Gold nanorods aspect ratio 10 


(a) Polystyrene nanospheres, 246 nm (b), Gold nanorods aspect ratio 4 (c) 


Backscattering signal [dB] 


Backscattering signal [dB] 


Backscattering signal [dB] 


0.4 0.4 


0.2 
Depth [mm] 


Figure 4.2: PS-OCT measurements and scanning-electron microscope (SEM) images of nano- 
particles with different shapes. (a) Sample S4: PS-OCT depth scan of polystyrene 
nanospheres with ds4 = 246 nm diameter in 0.5 wt.-% aqueous dispersion. Rm denotes polariza- 
tion-maintaining backscatter (light red: measurement, bright red: single-scattering model fit), 
and R. denotes cross-polarized backscatter. The polarization maintaining curve is dominating 
the cross-polarized trace by more than 20 dB, hinting to an aspect ratio of the particles close to 
unity. Fluctuations of the cross-polarized part are due to the measurement noise around 
-100 dB. (b) Sample Gl: PS-OCT measurement of gold nanorods with dg; = 25 nm diameter 
and aspect ratio of ARg; = 4 in 0.54 wt. % aqueous dispersion. The signal decay with depth 
corresponds to the total scattering and absorption cross sections of the gold nanorods. The 
cross-polarization curve (grey: measurement, black: single scattering fit) is below the polariza- 
tion-maintaining backscatter, and the polarization conversion amounts to 
nG1= —6.9 dB. The noticeable cross term is in accordance with the expected polarization con- 
version by scattering from non-spherical nanoscale objects. (c) Sample G2: PS-OCT measure- 
ment of gold nanorods with dg = 25 nm diameter and aspect ratio of ARg2 = 10 in 0.026 wt. % 
aqueous dispersion. In comparison to (b), the polarization conversion is stronger and amounts 
to 762 = —3.7 dB. (d-f) Scanning electron microscope images of the respective dried nanoparti- 
cles: (d) Sample S4: polystyrene nanospheres of ds4 = 246 nm diameter, (e) Sample Gl: gold 
nanorods of dg; = 25 nm diameter and ARa; = 4, (£) Sample G2: gold nanorods of doz = 25 nm 
diameter and ARg = 10. As can be seen in (e) and (f), the length of the rods is not completely 
uniform and the nominal values can be regarded as average values. 
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The decay coefficient of both curves is same, as expected in Eq. (4.4), and the 
polarization cross talk is measured to be 7g; = —6.9 dB. This noticeable cross 
talk is attributed to a clear deviation of the scatterers from a spherical shape, 
which is confirmed by the SEM image in Figure 4.2 (e). Compared to G1, the 
sample G2 shows even higher polarization cross talk 762 =—3.7 dB, see Fig- 
ure 4.2 (c). This behaviour is expected as the AR is increased to 10. The SEM 
image in Figure 4.2 (f) confirms the elongated particle shape. Further, a varia- 
tion of rod length around the nominal value of 256 nm can be noticed. 


The following table summarizes the determined scattering coefficients and 
cross sections as well as the polarization cross talk. 


Table 4.1: Scattering coefficients, cross sections and polarization cross talk, determined 
from the measurements in Figure 4.2 using the polarization-sensitive single scattering mod- 
el in Eq. (4.4). Sample S4 is used for calibration of 44 and 44m, for this reason the values of 
fland o are nominal values. 


Sample ih Ho,m a Obm n 

S4 307m” 0.0424m" 4.79x 10° nm” 6.62x10°nm* <-20 dB 
G1 1117m" 0.0331m! 1.96x10°nm” 5.80x 10 ’nm” -6.9 dB 
G2 1059m" 0010m! 963x10’nm 1.54x10'nm -3.7 dB 


4.4 Shape determination by comparison of 
simulation and measurements 


In the following, the measured scattering parameters are compared to simula- 
tion. For the determination of the particles’ size and shape, a set of simulations 
with varying geometries is performed and used for comparison. 


The selected electromagnetic simulation technique relies on a commercially 
available fully vectorial integral equation solver that is part of CST Micro- 
wave Studio [108]. The solver, which follows a Method of Moments (MoM) 
method, a rigorous full-wave technique for solving open-boundary electro- 
magnetic problems, is used for calculating the surface integral of the electric 
and magnetic field at the scattering particle. This method is especially 
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recommended for the calculation of scattering cross sections and thus fits per- 
fectly to our needs. An alternative numerical method, which allows scattering 
calculations of non-spherical particles is basing on the Discrete-Dipole Ap- 
proximation (DDA) [109]. In this method, the specific particle geometry is 
adapted by discrete dipoles of which the mutual excitation is taken into regard 
along with the excitation by the incident field. The technique is used for nano- 
rod investigations, but shows weaknesses in calculation of gold nanorods at 
infrared wavelengths [110], which is the regime of our analysis. A further cal- 
culation approach for non-spherical particles or small agglomerates thereof is 
to make use of the Rayleigh-Debye-Gans theory (RDG) [50]. In this approach, 
the arbitrarily shaped particle is subdivided into small Rayleigh scatterers, 
which are assumed to be unperturbed by each other. In a comparison [111] to 
other calculation techniques, RDG was found appropriate for limited size pa- 
rameters kl < 0.3 and refractive index differences An < 2, with the angular op- 
tical wavenumber k =27/ à and the particles longest dimension /. Although 
the RDG approach is appropriate for some types of non-spherical particles, the 
aforementioned conditions do not hold for our sample types. A promising 
method for calculation of scattering and extinction of particles as used in our 
experiment is the T-Matrix Method [112,113], which relies on the expansion 
of scattered and initial field in suitable vector spherical wave functions, and on 
linking the respective expansion coefficients in a so-called T-matrix. This ap- 
proach is applicable to many particle types and could be a further option in 
future calculations. 


For the CST calculations, the selection of the proper material properties of 
gold turned out to be crucial to the simulation results. Muskens [114] identi- 
fied the complex refractive indices measured by Johnson [115] to be fitting 
best to gold particles of nanometre size. These values were interpolated to a 
wavelength of 1.3 um and used for a Drude material model [116] within the 
simulations. In the simulation environment, the nanorods were built up from a 
cylindrical body and two hemispherical tails. We investigate a variety of dif- 
ferently shaped nanorods with diameters of 15 to 50 nm and aspect ratios of 1 
to 16. For these particles, the absorption cross section and the spatial and po- 
larization-dependent distribution of the scattering power were simulated and 
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averaged over random rod orientations to represent the scattering properties of 
nanorods dispersed in water. Averaging of the scattered powers from the indi- 
vidual particles is possible, since the particle concentration is small and light 
can thus be assumed to be scattered at most once. From integrating these an- 
gular scattering cross sections over the complete solid angle of 47, the scatter- 
ing cross section o, is obtained. The extinction cross section o;= 0, + Oy re- 
sults from taking the absorption cross section o, into account as well. The 
polarization-maintaining backscattering cross section Opm results from integra- 
tion of the polarization-maintaining angular scattering cross section in the 
backward direction over 1° half-angle, which corresponds to the numerical 
aperture of the system. The polarization-resolving simulation allows to detect 
the cross-polarized backscatter separately from the polarization-maintaining 
backscatter, which is used for calculation of the polarization cross talk 7. 


Figure 4.3 depicts the simulation results and compares them to the nanorod 
measurements. The subplots in Figure 4.3 (a-c) show simulated extinction and 
backscattering cross sections as well as polarization cross-talk as a function of 
diameter as well as AR. The term diameter is used here for the short axis of 
the particle, whereas ARxd gives the extent along the long axis, see Fig- 
ure 4.3 (f). For particle diameters between 15 and 40 nm, simulations were 
performed with a step size of 5 nm, and for diameters between 40 and 50 nm, 
the step size was increased to 10 nm. The step size in AR was 1 for 
AR= 1... 14, and 2 for AR = 14... 16, leading to 105 data points in total. For 
display and comparison to measurements, all results have been interpolated 
using a 2D spline interpolation. Figure 4.3 (a) shows the extinction cross sec- 
tion œ, Figure 4.3 (b) shows the polarization-maintaining backscattering cross 
section Op m and Figure 4.3 (c) shows the polarization cross talk 7. Within the 
simulation domain, o; generally rises with both particle diameter and extinc- 
tion ratio. In the vicinity of the plasmonic resonances around AR of 6 and 14 
the extinction is highest. For Om the trend is analogous, although not identical 
to the results in Figure 4.3 (a). The polarization cross-talk 7 is almost inde- 
pendent from diameter but shows strong dependence on AR, especially if AR 
is small. The black lines in the three subfigures indicate the values measured 
for sample G1 and G2. Table 4.2 summarizes the simulation results for the 
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geometries given by samples G1 and G2. The measured values, see also Ta- 
ble 4.1, are denoted for comparison. Measurement and simulation of o, and 
Ob,m Show agreement with respect to the order of magnitude. The polarization 
cross-talk 7 differs at most by 2.71 dB between simulation and measurement. 


Table 4.2: Comparison of simulated and measured extinction and backscattering coeffi- 


cients, Oj, Ob,m, respectively, as well as the polarization cross talk 77. 


Sample fon Obm 7 

G1 simulated 1.26 x 107 nm? 3.18 x 10 ° nm? —4.69 dB 
G1 measured 1.96 x 10° nm? 5.80 x 10” nm? -6.90 dB 
G2 simulated 6.22 x 10°’ nm? 2.31 x 10 nm? -3.67 dB 
G2 measured 9.63 x 10°’ nm? 1.54 x 10 nm” -3.70 dB 


In principle, the determination of pairs of d and AR for a set of two measured 
parameters, e.g. o; and Op m, is possible. However, in practice, each pair out of 
the three measured parameters (9, Obm, 77) would lead to a different solution 
for d and AR, since measurement and simulation provide a limited accuracy. 
For this reason, all three parameters are taken into account by calculating the 
(d, AR)-pair which provides in sum the lowest relative difference of simula- 
tion and measurement. First, the single ratios of each measured parameter to 
the simulated parameters are calculated. Since the extinction and backscatter- 
ing cross sections span a wide range of several magnitudes, comparison on a 
logarithmic scale seems to be favourable. The polarization cross talk 77 is ex- 
pressed in logarithmic dB-units already. Thus, the absolute logarithmic values 
of the ratios read 


fees) 
= t 
Ag, = 101g] ram 
t 
os) 
b, 
Ko, = 101g] |, (4.5) 
opm 
An= (meas) _ (sim) 
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where the superscript denotes either measurement (meas) or simulated (sim) 
data. The single ratios, Eq. (4.5), are divided by their individual maxima with- 
in the simulation domain, and after this normalization the total ratio 
is calculated 


A 


(4.6) 


1 Ao, Abm ‚An 
œ 3| max(Ao,) | max (Aoi, m ) max (A) | 


The factor of 1/3 keeps the result in the [0,1] range but does not affect the cal- 
culation of the particle shape. Figure 4.3 (d-e) shows Aw versus d and AR in 
colour-coded plots for both measurements of sample G1, Figure 4.3 (d), and of 
sample G2, Figure 4.3 (e). The white cross indicates the position of lowest to- 
tal deviation. In the case of sample G1 with nominal diameter of 25 nm and 
nominal AR of 4, see Figure 4.3 (d), Awt is low (blue coloured) in an area 
reaching from AR of 3 and diameter of 50 nm to AR of 6 and diameter of 
15 nm. The rod geometry with lowest deviation from simulation is found at a 
diameter of 26.75 nm and at AR of 4.1. This deviates from the nominal values 
by only 1.75 nm in diameter and 2.5% in AR. In the case of sample G2 with 
nominal diameter of 25 nm and AR of 10, see Figure 4.3 (e), Awt is low for 
AR > 6. In this case, the sensitivity in diameter is less pronounced. The small- 
est total ratio is found at a diameter of 18.8 nm and AR of 9. This amounts to a 
deviation of 6.2 nm in diameter and 10% in AR. 


Table 4.3 summarizes the OCT measurement results, which are in good 
agreement with the particle parameters obtained from SEM imaging. Note that 
the OCT measurement captures backscatter from a multitude of particles at 
once, thus averaging out deviations of individual particles from the nominal 
shape values. 
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(a) Extinction cross section a [nm?] (b) Backscattering cross section opm [nm?] (C) Polarization cross talk 7 [dB] 
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Figure 4.3: Particle shape determination by comparison of OCT-measured scattering parame- 
ters and electro-magnetic scattering simulations. Subfigures (a-c) show averaged scattering 
parameter simulations for an aqueous dispersion of randomly oriented gold nanorods. The hori- 
zontal and vertical axes denote aspect ratio AR and diameter d, respectively. (a) Total extinc- 
tion cross section o;, (b) polarization maintaining backscattering cross section oO} m, and (c) po- 
larization cross talk 7. Extinction cross section o; and backscattering cross section Ob m increase 
generally with particle diameter d and aspect ratio AR and show high values at the plasmonic 
resonances, e.g., around AR = 6. The crosstalk 7 increases with AR. The black lines indicate 
the values extracted from the OCT measurements of sample G1 and sample G2, which contain 
nanorods with nominal diameter dyominal,ı 2 = 25 nm and nominal aspect ratios ARyominaı,ı = 4 and 
ARyominal2 = 10. (d,e) Total ratio Aw of measurement and simulation comprising all three pa- 
rameters (Oj, Obm, 7). (d) Total ratio Aw of the simulations and the measured parameters ex- 
tracted from sample Gl, see Eq. (4.6). (e) Total ratio ^w of the simulations and the measured 
parameters extracted from sample G2, see Eq. (4.6). The white crosses show minimum total 
ratios and thus indicate the estimated particle diameter dmeasi2 and aspect ratio ARmeası,2. Ta- 
ble 4.3 lists the results for both samples. (f) Schematic drawing of both nanorod geometries G1 


and G2, along with the nominal diameters dpominai,1.2 and the respective aspect ratios AR pominal,! 
and AR pominal,2- 
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Table 4.3: Estimated particle diameters dmeas and aspect ratios ARmeas compared to the nom- 
inal specifications dpominal, ARnominal. The absolute and relative deviations between meas- 
urement and nominal values are given by ga and sar. 


Sample d & AR EAR 

G1 nominal 25 nm 

G1 measured 26.8 nm an) 4.1 ULOS W) 
G2 nominal 25 nm 10 


6.2 nm (25 %) 1.0 (10 %) 


G2 measured 18.8 nm 


4.5 Discussion 


The ramp-up of nanotechnology-driven industry demands fast and reliable 
techniques for quality control. In this work, we demonstrated that polarization- 
sensitive OCT offers precise size and shape measurement of nanoparticles, 
along with no sample preparation. The large probing volume of several cubic 
millimetres and the prospect of integrability into production lines, as well as 
the instantaneous output of characterization results are unique among the na- 
noparticle shape characterization tools. 


This work for the first time demonstrates quantitative characterization of size 
and shape of nanoparticles using PS-OCT. Nevertheless, the technology has 
some practical limitations, which we will discuss in the following. The tech- 
nique is constrained to the analysis of light scattering media, which are neither 
fully transparent, nor strongly absorbing, however the range of possible mate- 
rials is broad, as has been shown in the analysis of strongly absorbing CNT- 
loaded materials [10]. A further restriction is given by the aspect ratio, which 
should not exceed AR = 10, when exact quantitative shape determination is 
aimed for. Aspect ratios below this limit cover most practically relevant nano- 
particles from nanospheres to nanorods, but excludes, e.g., the differentiation 
of high-AR nanowhiskers from long nanorods. For larger aspect ratios, see 
Figure 4.3, the extinction cross section, the backscattering cross section and 
the polarization cross talk showed to be less sensitive to further aspect ratio 
changes. In our experiments, the particle size distributions have been narrow. 
Particles with broader size and shape variations could possibly be 
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distinguished by trend, but a quantitative shape determination down to nano- 
metre accuracy, as has been shown here, would be hampered. 


In conclusion, we have shown that PS-OCT allows to accurately determine 
size and shape of monodisperse gold nanorods having 25 nm diameter and 
lengths of 102 nm and 256 nm. We used a polarization-sensitive OCT setup, 
which we calibrated with respect to both power level and polarization. We 
measured three parameters: extinction cross section, backscattering cross sec- 
tion and polarization cross talk. Shape and size information were extracted by 
comparison to electro-magnetic scattering simulations of randomly dispersed 
particles, and the results are in fair agreement with the particle dimensions ob- 
tained from SEM-based imaging. We believe PS-OCT offers an attractive op- 
tion for simple and rapid nanoparticle characterization with reasonable instru- 
mentational effort as needed, e.g., in in-line or on-line quality control. 


4.6 Methods 


The swept-source PS-OCT system used in this work is illustrated in Fig- 
ure 4.4. In general, an OCT system measures the strength and the position of a 
multitude of scatterers along a light path in a sample. To this end, the electric 
field that is back reflected from the sample into the numerical aperture of the 
system is superimposed with a reference field originating from the same swept 
laser source. The superposition leads to optical frequency dependent interfer- 
ence patterns on a photodetector. The variation of signal amplitude with laser 
frequency is connected to the depth-dependent back reflection profile of the 
sample by a Fourier transform [28], see also Chapter 2.1.1. In our experiment 
we use a swept-source laser with central wavelength of 1315 nm and a wave- 
length scanning range of 1260-1370 nm, a scan rate of 1 kHz, and 10 dBm 
(10 mW) average output power (model s3, Micron Optics Inc., Atlanta, GA, 
USA). A first fibre-based directional coupler (CPL1) is used to split the power 
among the sample path (SP) and the reference path (RP). The reference path 
contains a free-space section allowing precise matching of the SP and RP 
lengths. Backscattered light from the sample and light travelling along the RP 
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is superimposed in a second fibre-based coupler (CPL2) and coupled to a po- 
larization-sensitive balanced photodetector (PS-BD, model INT-POL-1300, 
Thorlabs Inc., Newton, NJ, USA) offering 15 MHz electrical bandwidth. Po- 
larization controllers (PolC) are used for matching the polarization in both de- 
tector arms, and for setting the reference arm polarization to 45°, thus provid- 
ing the same reference arm power in both orthogonal polarizations. The PS- 
BD consists of two polarizing beam splitters (PBS) and a balanced detector 
(BD) for each of the polarizations. The balanced detection principle suppress- 
es relative intensity noise (RIN) from the strong RP light and enables a large 
dynamic range and a high sensitivity defined by the smallest measurable sam- 
ple reflectivity of -110 dB. The electrical signal is digitized by a 16 bit ana- 
logue-to-digital converter (ADC, model ATS660, Alazar Technologies Inc., 
Pointe-Claire, Canada) and processed by a personal computer (PC). The sam- 
ple path comprises fibre collimator (FC), rotatable quarter-wave plate (QWP) 
for polarization variation, a pair of galvo scan mirrors for lateral 2D beam 
scanning, a scan lens and the sample under test. The SP exists in two configu- 
rations: First the polarization calibration setup (top) with additional rotatable 
polarizer (Pol) and plane mirror instead of a sample, and second the measure- 
ment setup (bottom). 


For the describing the polarization of light within the PS-OCT system, we use 
the Jones formalism [117], where the electric field E (a,t) of polarized light is 
represented by the Jones vector J . The Jones vector consists of two complex 
elements Jy, J, with the normalized field amplitudes cos(a), sin(a) , of each 
of the orthogonal field components in x and y direction, as well as their com- 
mon phase ô and individual phases -ø and +9. Thus, the electric field is 


E(o,t)= Ale) F = A eie)? (e? cos(a)é, +e*!? sin(a)é,). 


The transmission through a linear optical system is described by a Jones ma- 
trix J, transforming input Jones vectors to output Jones vectors. If multiple 
optical systems are concatenated, their individual Jones matrices are multi- 
plied, and reverse transmission [69] of a unitary, optical reciprocal element 
with forward propagation matrix J results in 
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m =m; . my: my 
Jiev = , with J = : (4.7) 
-m m, m m, 
The aim of a PS-OCT measurement is the determination of the backscattering 


of the sample as a function of depth and polarization. This can be expressed by 
a Jones matrix 


(4.8) 


with the backscattering coefficient R;;(z) denoting backscattering of incident 
polarization i into backscattered polarization j, where i, j can each be x (hori- 
zontal) or y (vertical orientation). The division by the total backscattering 
Rs (z) is necessary to normalize the Jones matrix. According to Eq. (4.4), the 
sample backscatter of randomly oriented single-scattering particles only dis- 
tinguishes polarization maintaining R,,(z) from cross-polarized backscatter 
R.(z), thus in this case the sample Jones matrix can be simplified to 


_ 1 | Rm(z) Re (z) 
Lp nl mol a 


The measurement setup was built up with standard single-mode fibres (SMF) 
instead of polarization-maintaining fibres (PMF). In PMF, both principal po- 
larization axes feature different group refractive indices, which imply different 
group velocities, which in turn cause spurious signal peaks in the OCT meas- 
urement. Using an SMF instead avoids spurious peaks but does not maintain 
the incident polarization. Hence, a polarization calibration step is necessary 
before performing a sample measurement. This is done using the calibration 
setup depicted in Figure 4.4 (top). In this step, the Jones matrix Jsys describing 
the propagation of backscattered light through the fibre-optic system towards 
the PS-BD is determined as well as the input polarization Je of the light, 
which is emitted from the SP fibre towards the sample. 
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Figure 4.4: Setup schematic of the fibre-based polarization-sensitive swept-source optical co- 
herence tomography setup used in this work. SS: swept-wavelength source, scanning range 
1260-1370 nm; CPL1,2: fibre-based directional coupler; SP: sample path; RP: reference path; 
FC: fibre collimator; PolC: polarization controller; PS-BD: polarization-sensitive balanced de- 
tector, consisting of two polarization beam splitters (PBS) and two balanced detectors (BD); 
ADC: analogue-to-digital converter; PC: personal computer. For sample measurements, the 
sample arm comprises a free-space propagation section with a rotatable quarter-wave plate 
(QWP), scan mirrors, scan lens and sample. For calibration purposes, the sample is replaced by 
a plane mirror, and a rotatable polarizer (Pol) is added to the free-space section in the sample 
arm. Jones vectors and matrices are denoted in blue. J in: Jones vector of the input polarization 
leaving the FC towards sample or towards calibration setup; Jawrı,: Jones matrix of the 
QWP1,2; Jpoi: Jones matrix of the polarizer; Js: Jones matrix of the sample; Jsys: Jones matrix of 
the fibre-optic system from FC at the sample or calibration measurement head to the PBS at the 
PS-BD; J m Jones vector, measured by the PS-BD. 


109 


4 Nanoparticle size and shape characterization using PS-OCT 


The polarization measured at the PS-BD results from the concatenation of all 
transmitted optical elements 


Taci = Jys JowPi.sev (%) Josie (%) Ju Jpoı (%) Jowri (%) Jin > 
(4.10) 


which describes propagation of J in through the quarter-wave plate with 
Jowrı (39 )at an angle A of the fast axis with respect to the horizontal direc- 
tion, the polarizer with Jp, (Ip) at an angle % with respect to the horizontal 
direction, the back-reflection at the plane mirror described by Jm, the propaga- 
tion back through polarizer Jpoj rev (%) and QWPI JowPi.ev (%) and finally 
through the fibre optic system with matrix J,,, towards the detector. 


To determine both J,,, and J. multiple measurements of J 


in» | are neces- 


m,cal 
sary and performed at angular positions of the QWP1 of 99 = 0° to 165° in 
steps of 15° combined to angular positions of the polarizer of % = 0° to 150° 
in steps of 30°. With 72 measurements in total, the system is over-defined, and 
a least-squares based fit allows robust determination of both J,,, and Jin. The 
calibration is valid until fibres are moved or until a temperature drift causes a 
change of the fiber birefringence. In practice, the system is calibrated right 
before a measurement session. Here, all sample measurements use the 
same calibration. 


The configuration for the actual sample measurement is shown in Fig- 
ure 4.4 (bottom). In contrast to the calibration setup, the polarizer is removed 
and the mirror is replaced by the sample. The measured Jones vector is 
given by 


Jy = Joys J QWP2,rev (%) Js Jowr2 (%) Fine (4.11) 


Now, only J, is unknown, and for a robust determination thereof, J m iS meas- 
ured at various positions of the QWP2 from Jp = 0° to 165°, with a step size 
of 15°. A least-squares fitting procedure reveals J,, which we verified by using 
a linear polarizer as sample and determined its angular position with an accu- 
racy of 3°. 
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This chapter corresponds to publication [P1] with slight adaptations to match 
the format and the structure of this thesis. With the preceding chapters OCT 
has proven to be a powerful tool for the characterization of nanoparticles and 
nanocomposite materials. However, current systems are not suitable for many 
industrial environments, where often robust, small-sized and cost-effective 
systems are required. These challenges can be addressed by integrated-optical 
systems, with silicon as the material platform of choice, see Chapter 1.3. 


[Begin of Publication P1] 
Adapted with permission, © 2016 Optical Society of America. 
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Miniaturized integrated optical coherence tomography (OCT) systems have 
the potential to unlock a wide range of both medical and industrial applica- 
tions. This applies in particular to multi-channel OCT schemes, where scala- 
bility and low cost per channel are important, to endoscopic implementations 
with stringent size demands, and to mechanically robust units for industrial 
applications. We demonstrate that fully integrated OCT systems can be real- 
ized using the state-of-the-art silicon photonic device portfolio. We present 
two different implementations integrated on a silicon-on-insulator (SOI) pho- 
tonic chip, one with an integrated reference path (OCTint) for imaging objects 
in distances of 1 mm to 6 mm from the chip edge, and another one with an ex- 
ternal reference path (OCText) for use with conventional scan heads. Both 
OCT systems use integrated photodiodes and an external swept-frequency 
source. In our proof-of-concept experiments, we achieve a sensitivity of 
-64 dB (-53 dB for OCTint) and a dynamic range of 60 dB (53 dB for 
OCTint). The viability of the concept is demonstrated by imaging of biologi- 
cal and technical objects. 


5.1 Introduction 


Optical coherence tomography (OCT) [1,28] looks back on great advances in 
the past decades. Micrometer-range resolution in volumetric imaging opened a 
wide field of applications reaching from ophthalmology in medical diagnostics 
[87] to particle and defect characterization in materials science [80,88]. While 
current OCT systems are still assembled from discrete components, optical 
integration offers the prospect of compact low-cost mass-producible imple- 
mentations. Especially in endoscopic applications as well as in industrial sur- 
face and distance metrology, integrated systems can comply with the demands 
in size, mechanical robustness, temperature stability and cost. An arrangement 
of multiple parallel OCT systems becomes feasible, allowing the simultaneous 
recording of depth scans at different positions (B-scan). Silicon photonics is a 
particularly attractive platform for integration of OCT systems, lending itself 
to large-scale photonic-electronic integration on the basis of mature high-yield 
CMOS processes that are offered by widely available foundry services [118]. 
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Early in 2014, we introduced the first integrated silicon-photonic OCT system 
with surface imaging capabilities, comprising the interferometer and on-chip 
photodetectors [119]. In the same year, Yurtsever et al. demonstrated a sili- 
con-photonic OCT system and showed a sensitivity of -62 dB, yet with off- 
chip photodetectors [120]. Recently, a silicon photonic optical combiner with 
integrated photodiodes was used as a polarization-diverse receiver for a fiber- 
based swept-source OCT system [121]. However, for completing the OCT in- 
terferometer, external couplers and circulators are used. Other realizations of 
integrated OCT systems rely on silicon nitride and silicon oxynitride wave- 
guides, where both swept-source [122] and spectral-domain implementations 
have been shown. These demonstrations comprise, e.g., spectral-domain OCT 
chips based on arrayed waveguide gratings (AWG) [123], which have been 
used for characterization of biological samples [124]. Yet, silicon nitride and 
silicon oxynitride integration platforms are limited to passive components only 
such as waveguides and power splitters, and hence must be always comple- 
mented by external photodetectors. This limits the potential for large-scale 
integration. 


In this paper we demonstrate the first fully integrated silicon photonic OCT 
system, where receiver and interferometer components are monolithically in- 
tegrated on a single chip, and where only the swept-laser source (SS) and op- 
tional fiber-based delay lines are external to the chip. The experiments build 
upon the state-of-the-art silicon photonic device portfolio and demonstrate the 
advantages in size, cost and robustness inherent in monolithic integration. 
Based on our recent results [119] we investigate two system concepts, one 
with a fixed internal integrated reference path (OCTinr), and one with external 
reference path with variable length (OCText). With the compact OCTint sys- 
tem, a lens at the chip edge focuses the swept-source light to the sample and 
collects the backscattered radiation. The OCText system allows the use of a 
conventional OCT scan head positioned remote from the chip. We achieve 
sensitivities down to —64 dB (-53 dB for OCTint) and demonstrate imaging of 
both technical and biological samples. The sensitivities of our current OCT 
implementations are limited by relative intensity noise (RIN) and on-chip 
backscatter. Further improvement is possible by using optimized integrated 
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components and by reducing reflections at the chip edges. Each of the OCT 
circuits occupies an on-chip area of less than 0.4 mm’. 


This chapter is organized as follows: In Sect. 5.2 the concept and the imple- 
mentation of the OCTint and OCText systems is introduced. In Sect. 5.3 the 
sensitivity and the dynamic range of the systems are discussed. Sect. 5.4 is 
dedicated to the evaluation of the performance of both systems by demonstrat- 
ing two-dimensional and three-dimensional imaging of technical and biologi- 
cal samples. A summary and an outlook are given in Sect. 5.5. 


5.2 Silicon photonic OCT systems and 
experimental setup 


Both the OCTint and OCText systems use the principle of swept-source OCT 
[28]. In general, an OCT system measures the position and the strength of a 
multitude of scatterers along a light path in a sample at various distances z - z, 
from a reference plane at z,. To this end, the electric field E, that is backscat- 
tered from a sample is compared in amplitude and phase to a reference field 
E,. Both the sample and the reference field are derived from the same external 
optical swept-wavelength source. In an interferometer arrangement, sample 
and reference fields are superimposed on a (balanced) photodetector, the out- 
put current of which contains patterns resulting from interference of the 
backscattered field with the reference field. From these patterns, the amplitude 
and the position z — z; of the backscatter along the light path can be obtained 
by Fourier analysis of the photocurrent [28]. The backscatter strength denotes 
the ratio of backscattered optical power to optical power incident on the sam- 
ple. For an optimum sensitivity, the optical length difference between the re- 
turn path from the sample and the reference path should be small compared to 
the coherence length l, of the source [47]. 


The integrated OCT configurations presented in this work are depicted in Fig- 
ure 5.1 and Figure 5.2 and are described in more detail in the subsequent sec- 
tions. In each of the implementations, the photonic integrated circuits (PIC) 
comprise passive waveguides and 3-dB couplers (CPL) based on multimode 
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interference (MMI), forming an interferometer structure which is connected to 
on-chip germanium photodiodes (PD). Silicon photonics allows particularly 
compact implementations: Both OCT systems were co-integrated on the same 
chip along with a large number of further optical circuits for various applica- 
tions — the on-chip area occupied by each of the OCT circuits is less than 
0.4 mm”, and is dominated by the contact pads of the photodetectors which 
could be further reduced in size. As an external light source, we use a com- 
mercially available swept-source laser (Santec HSL-2100-ST) featuring a cen- 
ter wavelength of 1315 nm and a scanning range of 110 nm. The source is 
coupled to the PIC by either a lensed fiber (LF) or a polymer microlens (PL). 
All waveguides leading to chip facets are equipped with tapers narrowing to- 
wards the facet, designed for the full width at half maximum (FWHM) spot 
size of 1.9 um of the lensed fibers. An external polarization controller (PolC) 
is used to precisely excite the quasi-TE mode of the on-chip waveguides. The 
interference signal is recorded by two balanced on-chip photodiodes, which 
are specified to have a bandwidth of 20 GHz. The electrical signal from the 
photodiodes is tapped by an RF probe, combined for balanced detection, sub- 
sequently amplified, recorded by an analog-to-digital converter, and further 
processed on a computer. While the OCTint system as described in Sect. 5.2.1 
has the reference path integrated on-chip and thus allows measuring samples 
directly in front of the chip edge, the OCText system as described in 
Sect. 5.2.2 has an external reference path and allows imaging with a standard 
OCT scan head, meters away from the chip. Both PIC were fabricated on a 
silicon-on-insulator (SOI) wafer with a 2 um thick buried oxide and a 220 nm 
thick device layer using 248 nm lithography in the framework of a standard 
multi-project-wafer (MPW) run. For fabrication, the OpSIS service was made 
use of through JME A*STAR in Singapore, using design elements from the 
OpSIS OI25A PDK [38]. 


5.2.1 OCT chip with internal integrated reference path 
(OCTint system) 


The experimental setup and the silicon photonic integrated circuit (PIC) of the 
OCTint system are shown in Figure 5.1 (a). The reference path (RP) of the 
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OCT interferometer is integrated on the chip. The length of the reference path 
is designed such that the reference position z, of the depth scan coincides with 
the output facet of the PIC, i.e., sample and reference light interfering at the 
power combiner CPL2 have traveled the same distance through integrated on- 
chip transport waveguides. 


(a) [ss | oO SMF PU sp | 
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Figure 5.1: OCTint system configuration with integrated reference path: Experimental setup 
and photonic integrated circuit (PIC). (a) Schematic of setup. SS: swept-source laser, PolC: 
polarization controller, LF: lensed fiber, CPL1,2: 3 dB couplers with ports designated by 1...4, 
BL: ball lens, SP: sample path, RP: reference path, PD1,2: photodiodes with anodes (A1,2) and 
cathodes (C1,2), RF amp: RF amplifier, ADC: analog-to-digital converter, PC: personal com- 
puter. The photodiodes are contacted with RF probes and their photocurrents are subtracted for 
balanced detection. (b) PIC microscope image with optical input (LF), optical port with free- 
space path (via BL) to and from the sample, along with electrical connections (via RF probes). 
The OCTint system was co-integrated with a large number of additional optical circuits used 
for other purposes — the occupied on-chip area is indicated by a red frame and amounts to less 
than 0.4 mm’. 
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All transport waveguides are uniform and the associated modes feature identi- 
cal effective refractive indices which are invariant along the propagation di- 
rection. This leads to inherent compensation of dispersion in the on-chip sec- 
tions of the reference and the sample path. The swept-source laser has 
+11 dBm (12.6 mW) of output power and is edge-coupled to the PIC by a 
lensed fiber having a 1.9 um FWHM-spot size diameter, matched to the spot 
size of the waveguide tapers. This leads to an attenuation of 5 dB at the input 
facet and to an on-chip power of approximately +6 dBm (4 mW), averaged 
over a 100 nm wavelength span from 1270 nm 
to 1370 nm. 


A polarization controller (PolC) aligns the polarization according to the quasi- 
TE mode of the chip waveguides. In the sample path, a silica ball lens (BL) 
with a diameter of 1 mm is used to focus the radiation emerging from the chip 
facet and to collect the backscattered light from the sample. Figure 5.1 (b) 
shows a micrograph of the PIC along with the lensed fiber, the ball lens and 
the RF probe. 


5.2.2 OCT chip with external reference path 
(OCText system) 


The experimental setup and the silicon photonic integrated circuit (PIC) of the 
OCText configuration with external reference path are depicted in Fig- 
ure 5.2 (a). In contrast to the OCTint structure in Figure 5.1, the reference path 
(RP) is routed out of and fed back into the chip. The external reference path 
with adjustable length allows a long sample path (SP) and the use of a conven- 
tional scan head consisting of a fiber-collimator (FC) with angled physical 
contact connector (APC), a pair of galvo scan mirrors and an OCT scan lens. 
The reference path comprises a free-space section for precise length adjust- 
ments. Polarization controllers (PolC) align the polarization according to the 
quasi-TE mode of the chip waveguides. Efficient fiber-chip coupling is ena- 
bled by a novel kind of microlenses that can be fabricated in-situ by using di- 
rect-write three-dimensional laser lithography [125]. Figure 5.2 (b) shows a 
microscope image of the chip facet with attached polymer lenses. Fiber-chip 
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coupling loss was estimated by a combination of transmission and photocur- 
rent measurements. The lensed ports on the right-hand side couple to standard 
single-mode fibers (SMF) with typical losses of 4 dB (6 dB at maximum) over 
a 100 nm wavelength span from 1270 nm to 1370 nm. The lensed port at the 
left-hand side has a loss between 7 dB and 10 dB in the same wavelength 
range, due to an offset error between lens center and chip waveguide. With 
this coupling scheme and an available laser source power of +11 dBm 
(12.6 mW), we estimate a wavelength-averaged on-chip power of +3 dBm 
(2 mW). The average power incident on the sample was measured to be 
—5 dBm (0.3 mW), and the average reference path power received by each of 
the photodiodes is estimated to be -17 dBm (20 uW). 
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Figure 5.2: OCText configuration with long external reference and sample paths: Experimental 
setup and photonic integrated circuit (PIC). (a) Schematic of setup. SS: swept-source laser, 
PolC: polarization controller, SMF: single-mode fiber, PL: polymer lens, CPL1,2: 3 dB cou- 
plers with ports designated by 1...4, SP: sample path, RP: reference path, PD1,2: photodiodes 
with anodes (A1,2) and cathodes (C1,2), APC-FC: fiber collimator (FC) with angled physical 
contact connector (APC), RF amp: RF amplifier, ADC: analog-to-digital converter, PC: per- 
sonal computer. The photodiodes are contacted with RF probes and their photocurrents are sub- 
tracted for balanced detection. (b) PIC microscope image of the chip’s right edge. The OCText 
system is co-integrated with a large number of additional optical circuits used for other purpos- 
es — the occupied on-chip area amounts to less than 0.4 mm”. SP: silicon waveguide for sample 
path, C1,2 and Al,2: contact pads for photodiode readout with RF probes, PL: Three polymer 
lenses between chip waveguides and a standard single-mode fiber array. 
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5.3 Sensitivity and dynamic range 


Noise and parasitic backscatter limit both the sensitivity and the dynamic 
range (DR) of the system. The sensitivity is defined as the minimum measure- 
able sample reflectivity, which is determined by the noise level of the system. 
The dynamic range, in turn, denotes the maximum achievable signal-to-noise 
power ratio (SNR) and is smaller than or equal to the sensitivity. In our im- 
plementations, noise is dominated by either the relative intensity noise (RIN) 
of the swept source or by parasitic backscatter — thermal noise of the photodi- 
ode RF amplifier is irrelevant in all cases. 


Figure 5.3 illustrates the different signal and noise contributions as applicable 
to the two OCT systems. In general, the relative intensity noise of a laser 
source relates the variance of the optical power fluctuations to the mean power 
squared [46]. The relative intensity noise is generally described by its one- 
sided power spectrum RIN( f- fe) centered at the carrier frequency f., and 
can be obtained from a spectral analysis of the photocurrent in the baseband. 
The spectrum RIN( f- fa) is specified in units of Hz", or, in logarithmic 
form, as dB Hz". The total RIN results from an integration over all frequen- 
cies above the carrier frequency, 


RIN ot = | RIN f=Te df, 
! ( ) 6.1) 


RIN (f - fe) =101g(RIN(f - f,)x1Hz) (in dBHz™'). 


In swept-source OCT, balanced detection is normally used. The photocurrent 
then only exhibits RIN contributions that result from the interference of the 
sample field and the reference field, whereas the individual RIN fluctuations 
of the sample and the reference signal cancel if ideal balancing is assumed 
[126]. As a consequence, a reflection peak at z =z, is surrounded by an ac- 
companying noise profile RIN, (z = Ze). Assuming an idealized SS laser for 
which the optical frequency increases linearly with time, the depth offset 
Z — Z, from the reflection peak is proportional to a baseband frequency offset, 
which in these experiments is f- f, = (z-z.)x3.55 MHz/mm. This value 
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was determined by inserting a mirror into the sample path and by measuring 
the frequency shift of the photocurrent spectrum in MHz as a function of mir- 
ror position in mm. Using the relation RIN, (z-z,)=RIN(f)xdf/dz, the 
spatial relative intensity noise profile is then given by 


RIN, (z2-2,)=RIN((z- 2,)x3.55 MHz/mm)x 3.55 MHz/mm. (5.2) 


When expressed in logarithmic form, it is convenient to consider the RIN con- 
tribution within a spatial element that corresponds to the depth resolution of 
öz=8 um of the OCT system in vacuum, 


RIN, ag(Z—Z_)=10lg(RIN,(z-z,)x8um)  (indB($um)'). (5.3) 


For a linear sweep of the optical frequency, the depth resolution of öz=8 um 
is equivalent to a z x 3.55 MHz/mm=28.4kHz frequency resolution of the 
photocurrent spectrum. The position-dependent RIN, ( Zs Z) of the SS laser 
can be directly measured using a conventional fiber-based OCT setup, see 
Figure 5.3 (a). To this end, a fully reflecting (full refl.) mirror is introduced at 
position z,—Z,=Ilmm, leading to a spatial relative intensity noise profile 
RIN, ag (2-2,) centered at z.-2,=1mm (thin blue curve). The thick blue 
curve provides a schematic of these data for a spatial resolution of öz=8 um. 
In the immediate neighborhood of this reflection peak at z,—Z,=0Omm, we 
measure RIN, ag (0)=—60 dB (8 um) | , which is by 60 dB smaller than the 
reflection peak. Further away from the reflecting mirror, e.g., at 
Zy—z,=8mm, we find RIN, gp(z—z,)=—80 dB(8 um), which is 
smaller than the reflection peak by 80 dB. If the mirror is replaced by a 
—50 dB reflector, the noise is about 100 dB down (gray curve), dominated by 
residual RIN due to non-ideal balancing of the photodetectors. The spurious 
peaks in the gray curve in Figure 5.3 (a) stem from weak reflections in the set- 
up. As a comparison, we measure the one-sided relative intensity noise spec- 
trum RIN( f ) of the Santec swept-laser source with a photodiode and an elec- 
trical spectrum analyzer. We find  RINgg(f,)=—105 dB Hz"! 
=-60 dB (28.4 kHz) ' for f,=3.55 MHz, which corresponds to 
RIN, as (Za — z1) =-60 dB (8 um) at a distance of Z, — zı =1 mm from the 
peak of a fully reflecting mirror, in perfect agreement with the results obtained 
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from the OCT measurement. Similarly, at a larger baseband frequency 
fo = 30 MHz, we measure RINgp (fp) =—125 dB Hz™ =-80 dB 
(28.4 kHz)", corresponding to RIN, u (2-2) =-80 dB (8 um)" at a 
distance of Z} — Z,=8.5 mm from the peak. 


Besides RIN, distributed on-chip backscatter can be a reason for sensitivity 
limitations of integrated OCT systems. On-chip backscatter is mainly caused 
by random sidewall roughness or other irregularities of high index-contrast 
waveguides, which couple forward and backward-propagating fundamental 
modes to each other [127]. Interference of the backscattered light with the ref- 
erence light or with light from spurious reflections on the PIC leads to a con- 
tinuous measurement background. This unwanted interference signal may be 
suppressed by the common-mode rejection of a balanced detector if both 
backscattered light and reference light are coupled to the same input port of 
the combiner CPL2 in front of the balanced detector. Coherent backscatter is 
hence only visible in case of imperfect common-mode rejection or if backscat- 
tered light and reference light are coupled to different input ports of the com- 
biner. 


OCTint In case of the OCTint system, Figure 5.1 and Figure 5.3 (b), the 
chip edge reflection (Zedge,int — Zr =9 mm) amounts to around 3 % (-15 dB) 
and is by 24 dB larger than the light returned via free-space propagating from 
a fully reflecting (full refl., return loss 0 dB) plane mirror. This corresponds to 
single-pass coupling losses between the chip and the free-space beam of ap- 
proximately 20 dB. Since the depth information in SS-OCT results from a 
Fourier analysis of the real-valued time-domain data, negative depths cannot 
be discriminated from the true positive depths (depth degeneracy). Therefore 
the on-chip backscatter appears also at depths z — z; > 0. In the OCTint system, 
distributed on-chip backscatter (black) dominates over the RIN associated 
with the strong reflection from the chip edge for small depths z — z; < 5 mm, 
corresponding to the optical path lengths of the on-chip waveguides. For big- 
ger depths z — z, > 5 mm, the influence of on-chip backscatter can be neglect- 
ed, and the sensitivity is limited by the RIN background (blue) caused by the 
chip edge peak. This RIN background is RIN, gp (Zg — 20)=-80dB (8 um) | 
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below the maximum of the chip edge peak, leading to a sensitivity of —56 dB. 
The total background consists of the sum of on-chip backscatter and RIN (am- 
ber dotted line). Since relative intensity noise originating from the chip edge 
dominates all noise contributions that could possibly come from other objects 
(red: full refl.; green: low refl.), the DR is 56 dB. 


For the OCTint system, sensitivity can be drastically improved by reducing 
the chip edge reflections using, e.g., angled chip facets [128], or by improving 
coupling efficiency to and from the chip, using, e.g., polymer microlenses as 
for the OCText system. Assuming a reduction of facet reflectivity to be- 
low -30 dB and a reduction of the single-pass coupling losses to less than 
5 dB, as in the OCText system, we believe that sensitivity improvements of 
approximately 45 dB can be achieved. 


OCText The situation is different in the case of the OCText system, Fig- 
ure 5.2 and Figure 5.3 (c), where the chip edge peak at Zedgeext is far outside the 
source coherence window (long reference path, Zedge,ext — Zr = —8 m). The asso- 
ciated RIN, (zı = Caer) in the measurement window near zı — Zr = 1 mm 
corresponds to RIN(28.4 GHz), which is many orders of magnitude smaller 
than the RIN levels measured at MHz frequencies as described above. For the 
OCText system, RIN from the chip edge reflection can hence be safely ne- 
glected as a relevant source of signal impairment. Moreover, on-chip 
backscatter (black) is less prominent for the OCText system than for the 
OCTint system. Due to the long fiber-based RP and the limited coherence 
length of the source, on-chip backscatter does not coherently interfere with the 
reference light. For a perfectly balanced pair of photodiodes, on-chip 
backscatter should hence not be visible at all. This is different for the more 
realistic case of imperfect balancing, where the photocurrent may also contain 
signal components that result from interference of optical signals that are cou- 
pled to the same port of CPL2 of the OCText system. In this case, spurious 
reflections at the chip edges or at CPL 1 may act as a parasitic reference signal 
and generate an interference signal with distributed on-chip backscatter, thus 
limiting the sensitivity. In our implementation, the parasitic reference signal is 
dominated by light reflected at the chip edges of both SP and RP. 
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To obtain an estimate of the signal background associated with on-chip 
backscatter, we perform an independent measurement to determine the overall 
on-chip backscatter of the SP and RP waveguides to approximately 
Pwc.ap = 50 dB within the depth resolution of 8 um. To translate this into a 
signal background in the OCT measurement, we need to estimate the strength 
of the parasitic reference signal. To his end, the chip edge reflection of each of 
the SP and RP waveguide ports is assumed to amount to -15 dB (3 %) as in 
the case of the OCTint system, and the on-chip input power is estimated to 
+3 dBm (2 mW). Taking into account 1 dB of excess loss for each of CPL1 
and CPL2, the parasitic reference light coupled to each of the photodiodes thus 
amounts to —21 dBm (8 uW), which is öyg =4 dB less than the reference light 
from the off-chip reference path. Moreover, the parasitic reference light and 
the on-chip backscatter are coupled through the same port | to CPL2, and their 
interference signal is therefore attenuated by the common-mode rejection ratio 
(CMRR) of the balanced receiver, which amounts to Zug =15 dB. Taking into 
account the insertion loss of the lensed ports towards the sample, agg =4 dB, 
and referring the backscatter to an ideal mirror in the off-chip sample path, the 
backscatter floor in the OCT measurement can be estimated to be 
PBc.aB = Pwa.aB + 244p — gn — Gag = —61 dB. This backscatter floor limits the 
sensitivity to — 61 dB and is the relevant limitation for low sample reflectivi- 
ties (green: low refl.). For strong reflectivities, the total background consists of 
the sum of on-chip backscatter and RIN (amber dotted line) associated with 
the reflection peak. For a fully reflecting mirror (red: full refl.) the associated 
RIN, ae (21-2) =-60dB (8 um) "determines the DR of 60 dB. 


For the OCText system, the sensitivity limitations can be mitigated in various 
ways: On-chip backscatter can be reduced by optimizing the fabrication pro- 
cesses of the waveguides — our current experiments were based on a standard 
silicon photonic MPW process for optical communications, where backscatter 
has not been considered a crucial aspect that requires special attention. More- 
over, the chip edge reflections can be reduced to values well below —30 dB, by 
using, e.g., angled waveguide facets [128], which would lead to lower parasit- 
ic reference signals. Additionally, improved balancing of the coupler (CPL2) 
in front of the balanced detectors leads to an increased CMRR and thus further 
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improves the sensitivity of the OCText system. Using more advanced power 
combiners with tunable splitting ratio, e.g., based on Mach-Zehnder interfer- 
ometers and phase shifters, we believe that CMRR values of significantly bet- 
ter than 30 dB can be achieved, leading to overall sensitivity improvements of 
more than 30 dB. 
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Figure 5.3: Sensitivity and dynamic range (DR) derived from a measured OCT scan in a con- 
ventional setup, and from A-scan schematics for configurations OCTint and OCText. Horizon- 
tal axes: Measurement depth z — z, referred to a reference position zr. Vertical axes: Backscatter 
signal and noise power relative to the power reflected from an ideal mirror at position zy. All 
relative intensity noise data are related to the OCT resolution bandwidth of 28.4 kHz, or, equiv- 
alently, to the OCT depth resolution of 8 um. (a) Measurement of RIN background of the swept 
source used for the OCTint and OCText system: The profile of the RIN, ıs background is de- 
rived from the measured backscatter signal of a fully reflecting mirror (full refl.) positioned at 
Zı— Zr = 1 mm in a conventional fiber-based OCT setup (thin blue curve: measurement, thick 
blue curve: schematic approximation). The spatial RIN, ap(z— Ze) profile for a resolution of 
8 um is about 80 dB down and results from the interference of reflected source RIN with the 
reference field. If the fully reflecting mirror is replaced by a partially reflecting mirror (power 
reflection factor —50 dB), the noise level corresponds to a reflectivity of —100 dB (gray curve, 
spurious peaks originate from weak multiple reflections in the setup), because the balanced 
receiver suppresses the source RIN, which is transmitted now mainly along the reference path. 
(b) Noise and backscatter background for the OCTint setup. At small distances z — zr < 5 mm, 
on-chip backscatter (on-chip backsc., black) is dominant. At larger distances z — zr > 5 mm, RIN 
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from the strong reflection at the chip edge (blue, Zedge, int = Zr) dominates the noise level for any 
sample reflectivity. The total background consists of the sum of on-chip backscatter and RIN 
(amber dotted line). Both a fully reflecting mirror (full refl., red) and a partially reflecting mir- 
ror (low refl., green) lead to similarly shaped, but shifted RIN contributions which are smaller 
than the RIN related to the chip edge reflection peak. The minimum detectable reflectivity (sen- 
sitivity) is -56 dB and represents also the DR of the system. (c) Noise and backscatter back- 
ground for the OCText setup with schematic backscatter curves. The chip edge peak at Zedge, ext — 
Zz is far outside the source coherence window (long reference path, 
Zedge, ext — Zr = 8 m). However, on-chip backscatter (on-chip backsc., black) falls within the 
coherence window and limits the sensitivity to -61 dB. RIN from a fully reflecting mirror (full 
refl., red) is responsible for the noise floor 60 dB below the peak. The DR is therefore 60 dB. 
The total background consists of the sum of on-chip backscatter and RIN (amber dotted line). A 
weaker mirror reflectivity (low refl., green) can be measured as long as it is larger than —61 dB. 


5.4 Performance evaluation and 
application demonstration 


The performance of the integrated OCT systems was assessed by measuring 
sensitivity, dynamic range and depth scanning range with a plane mirror as a 
test object placed at various depth positions z — z,. In addition, two- and three- 
dimensional OCT images of technical and biological samples are taken. 


5.4.1 OCT chip with internal integrated reference path 
(OCTint system) 


The OCTint system is investigated by using a movable plane mirror as a sam- 
ple. For calibrating the backscatter measurements, a fully reflecting sample 
mirror is placed at a position z-z, =8mm away from the chip edge 
Zedge,int — Zr = 0. Figure 5.4 (a) depicts measured depth scans. For various mirror 
positions between 5 mm and 12 mm (brightly colored scans) the scans are su- 
perimposed, each curve resulting from averaging 100 scans. Because the depth 
information results from a Fourier analysis of the scan traces in the frequency 
domain, negative depths z- zr; < 0 located inside the chip cannot be discrimi- 
nated from the true positive depths outside the chip (depth degeneracy). 
Therefore the on-chip backscatter (black part of the scan) appears also in 
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space-inverted form at depths z — z; > 0 (gray part of scan), and the backscatter 
from the mirrors could be also seen at z — z, < 0 (lightly colored scans). Note 
that the depth axis was derived assuming light propagation in vacuum, such 
that internal backscatter appears over a 5 mm long depth range even though 
the scattering on-chip waveguides are approximately 1.4 mm long. The meas- 
urement range is hence limited to 2-2, >5mm. Because the edge reflection 
cannot be removed, the associated fluctuations due to relative intensity noise 
determine the noise floor in the measurement window and therefore the sensi- 
tivity, as has been explained in Sect. 5.3 and Figure 5.3 (b). For mirror posi- 
tions between 5 mm and 10 mm, the noise level is measured to be between 
50 dB and 53 dB below the mirror reflection leading to a dynamic range of up 
to 53 dB. The minimum measureable sample reflectivity would be —53 dB and 
defines the measurement sensitivity. This is in good agreement with the 56 dB 
DR and the -56 dB sensitivity estimated for the OCTint system in Section 5.3. 
The spatial variation of the measured peak reflectivities is attributed to the fi- 
nite coherence length of the swept source and the finite Rayleigh range of the 
focused Gaussian beam after the ball lens. The spatial autocorrelation function 
of the swept-source laser in vacuum is shown as a black broken line in Fig- 
ure 5.4 (a). A 10 dB coherence length of l; =2x6 mm can be inferred. The 
drop in measured reflectivities for positions outside the focus is described by a 
defocusing function, black dotted line in Figure 5.4 (a). If the ball lens focus is 
positioned properly, here at 10 mm from the chip facet, the resulting defocus- 
ing function partially compensates the signal drop caused by the limited co- 
herence length. 


The axial 3 dB-resolution was measured to be 11 um, which is slightly worse 
than the 7 um estimated from the wavelength scanning range of the source. 
Since dispersion is inherently compensated by equal lengths of on-chip wave- 
guide sections, the reduced resolution is attributed to an imperfect knowledge 
of the frequency-time relation of the swept-frequency source. In general, the 
frequency of the tunable laser does not increase linearly with time during the 
rising slope of a frequency sweep. This can be compensated by digital signal 
processing if the nonlinearity of the frequency sweep is known exactly. To this 
end, the time-dependent increase of frequency must be extracted from a 
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calibration measurement using a mirror as sample reflector [129]. Ideally, for 
a system with no dispersion imbalance between the sample and the reference 
path, this time-frequency mapping should lead to an axial resolution that is 
only limited by the tuning range of the swept source. In practice, however, re- 
sidual errors of the measured time-dependent increase of frequency lead to a 
depth-dependent broadening of the reflection peaks in z-space. 
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Figure 5.4: Backscatter measurements on the OCTint chip with integrated reference path. 
(a) Reflections from a plane mirror at 8 (positive) depth positions z — zr. A depth z — zr = 0 cor- 
responds to the chip edge. The black broken line is the spatial autocorrelation function of the 
light source with a 10 dB coherence length of l = 2x6 mm in vacuum. The beam is focused by 
a ball lens (BL) at z — zr = 10 mm. The defocusing function (black dotted line) describes the 
depth-dependent variations of the power reflected from a plane mirror and coupled back into 
the on-chip waveguide. The drop in spatial coherence (black broken line) partially compensates 
the defocusing function. The 8 axial A-scans of the mirror backreflections are superimposed 
(brightly colored scans, normalized to the maximum of the 8 backreflections). The finite isola- 
tion of ports 1-4 and 2-3 of coupler CPL1, multiple reflections, and irregularities inside the chip 
lead to backscatter at negative depths (z — zr < 0, black part of scan). Because the depth infor- 
mation results from a Fourier analysis of the scan traces in the k-domain, negative depths can- 
not be discriminated from the true positive depths (depth degeneracy). Therefore the on-chip 
backscatter appears also at depths z—z, > 0 (gray part of scan), and the backscatter from the 
mirrors could be also seen at z — zr < 0 (lightly colored scans). The backscatter in combination 
with the depth degeneracy limits the measurement range to z — zr > 5 mm. The measured mirror 
reflections are uniform with variations of less than 3 dB for mirror positions between 5 mm and 
10 mm. The noise is measured to be between 50 dB and 53 dB below the mirror reflection lead- 
ing to a dynamic range of up to 53 dB. A minimum measureable sample reflectivity of -53 dB 
defines the measurement sensitivity. (b) Two-dimensional B-scan of a piece of pumice. The 
gray line emphasizes the sample surface. The porous surface structure is clearly visible on both 
the B-scan and the photograph of the cross-section (inset). 


127 


5 Optical coherence tomography system on a silicon photonic chip 


The imaging capability of the OCTint system is tested by two-dimensional 
measurements of a piece of pumice, Figure 5.4 (b). Strong scattering prevents 
measurable reflections from deeper ranges inside the material, but the porous 
surface structure is clearly visible (see the photograph in the inset). 


5.4.2 OCT chip with external reference path 
(OCText system) 


The OCText system is also investigated by using a movable plane mirror as a 
sample. For calibrating the backscatter measurements a sample mirror with a 
1 % reflection factor (20 dB return loss) is placed at a position z — zr = 1 mm 
within the measurement window 1 mm < z — z, < 10 mm. The chip edge is lo- 
cated at Zedge, ext — Zr =—8 m and hence does not influence the measurements. 
Figure 5.5 (a) shows the superposition of single depth-scans (A-scans) for a 
single —20 dB sample reflector placed at 1 mm increments within the meas- 
urement window. The sample reflector comprises a plane mirror and a neutral 
density filter with 10 dB single-pass attenuation. Each curve results from av- 
eraging 100 scans. The noise level, which defines the sensitivity of the system 
with respect to a 0 dB reference return loss, is measured to be —64 dB. This is 
in fair agreement with the sensitivity of — 61 dB estimated in Sect. 5.3 based 
on an analysis of effects of spurious on-chip backscatter. 


The imaging depth of the OCText system is mainly limited by three effects: 
The limited coherence length of the source (black broken line), the defocusing 
of the measurement beam which leads to a drop of power coupled back into 
the fiber for positions outside the focus (defocusing function, black dotted 
line), and by residual errors of the time-frequency mapping of the swept 
source. These residual errors lead to a depth-dependent broadening and a de- 
crease of the maxima of the reflection peaks in z-space. The total decay caused 
by limited coherence, defocusing and imperfect time-frequency mapping 
amounts to approximately 3.3 dB / mm, see Figure 5.5 (a). 


The depth resolution of the system is 30 um, which is significantly worse than 
the 7 um expected from the tuning range of the swept source, also when tak- 
ing into account the effects of imperfect time-frequency mapping. The rather 
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bad axial resolution is caused by unbalanced dispersion in the sample and ref- 
erence path. This originates from the fact that the sample and the reference 
path are composed of different fiber lengths and different path lengths of the 
on-chip SOI waveguides. Dispersion-induced resolution limitations can be 
overcome by carefully designing the dispersion maps of the two paths. 


Using a perfect reflector, the measured dynamic range is 60dB at 
Z—Z,= 1mm (not shown) and is limited by the laser RIN. This DR agrees 
very well with the DR of 60 dB estimated in Sect. 5.3, see also Figure 5.5 (c). 
The dynamic range is not limited by the RIN originating from the reflection at 
the chip edge, which is located at Zedge, ext — Z = —8 m far outside the measure- 
ment range. The corresponding RIN(28.4GHz) is not measureable. 


-20 dB reflector at 
different positions 


Backscatter [dB] 
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Figure 5.5: Backscatter measurements obtained with the OCText chip using an external refer- 
ence path. A depth z—z,=0 denotes an approximate geometrical position 8 m off the chip 
edge. (a) Axial scans of a —20 dB reflector as test sample placed at various depth positions z — 
Z, Showing a signal decay of 3.3 dB / mm. The black broken line is the spatial autocorrelation 
function of the light source with a 10 dB coherence length of l. = 2x6 mm in vacuum. The posi- 
tion 1 mm coincides with the focus of the scan lens. The black dotted line represents the influ- 
ence of defocusing on the reflected power coupled back into the on-chip waveguides. The sam- 
ple is placed at distances from 1 mm to 10 mm within the scanning range, and the resulting 
scans are superimposed. For each curve, 100 subsequent scans are averaged. The measurement 
sensitivity is —64 dB for objects placed in a region 0.5 mm < z — zr < 2 mm. (b,c,d) Three- 
dimensional C-scans of biological and non-biological objects. The images are averages of 100 
scans. (b) shows a piece of pumice, (c) a part of a decayed leaf of cornus sanguinea, and (d) a 
cross section of a reel of tape. 
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The applicability of the system to sample imaging is demonstrated with 3D 
volume scans of biological and non-biological samples. Figure 5.5 shows 
three-dimensional scans (C-Scans) of a piece of pumice (b), a part of a de- 
cayed leaf of cornus sanguinea (c), and a reel of adhesive tape (d). 


5.5 Summary and outlook 


We demonstrate silicon photonic swept-source OCT systems with integrated 
photodetectors. We design, implement and characterize two different OCT 
configurations with on-chip internal reference path (OCTint) and with external 
reference path (OCText). We demonstrate the viability of the systems by im- 
aging both technical and biological samples: Configuration OCTint enables 
imaging directly in front of the chip facet. For this system, we measure sensi- 
tivities between -50 dB and -53 dB over 5 mm axial scanning range. As a 
demonstration, we show a cross-sectional B-scan of a scattering pumice sam- 
ple. Configuration OCText supports the use of a remote standard OCT scan 
head. For this system, a sensitivity of -64 dB was achieved — the best value 
reported so far for a fully integrated optical OCT processing unit. The OCText 
system achieves a better sensitivity than the OCTint system, although a larger 
number of fiber-chip interfaces add to the coupling loss. However, our directly 
written [125] optimized polymer microlenses minimize this loss. For the 
OCTint setup, aberrations caused by the use of a standard ball lens lead to 
suboptimal coupling. The limiting factors for both systems are coupling losses 
and on-chip backscatter. The OCTint exhibits a measured dynamic range of 
53 dB, limited by the RIN contribution from the strong chip edge reflection. 
The dynamic range of the OCText is 60 dB, limited by the RIN contribution 
from a fully reflecting object. 


We believe that our miniaturized OCT systems are suitable for a wide range of 
both medical and technical applications, including also advanced schemes 
such as Doppler OCT [130], optical coherence elastography (OCE) [131], or 
photo-thermal OCT [132]. By combining silicon photonic OCT engines with 
chip-scale swept-source lasers [133], fully integrated OCT systems come 
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within reach. We report on a first proof-of-concept demonstration, with vast 
potential for improving the performance: Polymer microlenses, as used in 
OCText, could be used for the OCTint system as well, enhancing power cou- 
pling and reducing chip-edge reflections. Thus, a higher dynamic range and an 
improved sensitivity could be achieved with this system. For OCText, the ef- 
fects of on-chip backscatter can be significantly improved by optimized fabri- 
cation processes and device design. In both cases, sensitivities better than 
-90 dB can be achieved based on technological and design improvements. 
Moreover, we believe that compact and highly integrated OCT systems with 
tens or even hundreds of parallel spatial sensors will open completely new ap- 
plication spaces and create attractive market opportunities. 


[End of Publication P1] 
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In the course of this work, OCT systems have been developed, both fibre- 
based implementations and photonic integrated circuits on a silicon chip. The 
fibre-based system has been designed for either single-polarization or polari- 
zation-sensitive operation, as well as for in-line dispersion-state monitoring of 
nanocomposites during compounding. In addition to system setup, a light scat- 
tering model for nanomaterial analysis has been developed. Combined with 
Mie’s scattering theory accurate particle size determination has been demon- 
strated. An extension of the model towards different polarizations enabled the 
precise characterization of gold nanorods with respect to their diameter and 
aspect ratio. The micro and nanostructure of composite materials has been 
characterized in a probing volume of several cubic millimeters, using image 
processing techniques on OCT measured data for larger features, and by ap- 
plying scattering model analysis, which allowed to obtain nanoscale infor- 
mation on the sample structure. Besides material analysis, photonic integration 
of OCT systems during this work has led to silicon photonic OCT circuits, 
comprising interferometer and balanced detector on chip. These systems al- 
lowed imaging of biological and technical samples. In the following, the key 
achievements of this work are summarized and an outlook to possible future 
activities is given. 


For nanomaterial analysis with OCT, a polarization-resolved light scattering 
model has been developed. The model is adequate for weakly scattering sam- 
ples and links the scattering properties of the single particles to the backscat- 
tering signal obtained in the OCT measurement. This allows the analysis of 
scatterers beyond the resolution limit of OCT [P2, P3, P5]. At higher sample 
depths or for strongly scattering samples, multiple-scattered light adds to the 
measurement signal. Future work might address this effect by models that ac- 
count for multiple scattering. 


A polarization-sensitive optical coherence tomography system has been 
developed using standard single-mode fibres, free-space polarization optics, 
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and a polarization-sensitive receiver. For experiments during nanocomposite 
production, a dedicated OCT probe for measurements at a compounding ex- 
truder with rugged titanium shaft, adjustable optics and temperature and 
scratch-resistant sapphire window has been developed and applied to produc- 
tion process monitoring [P2]. 


Dispersion-state analysis of nanocomposites has been one main application 
of the light scattering model. The dispersion state, which is crucial to the func- 
tionality of nanocomposites, has been characterized based on both, image pro- 
cessing of OCT depth-scans and model-based evaluation of the backscatter 
signals. Image analysis has been applied to particle agglomerates which are 
bigger than the resolution of the system. For a quantitative evaluation, the im- 
age data has been segmented and two independent parameters have been ex- 
tracted: Area fraction and perimeter-to-area ratio, both of which show strong 
correlation with the dispersion effort in sample production [P2]. The scattering 
model analysis has been applied successfully to CNT-epoxy samples prepared 
in three-roll mills and even outperformed conventional light microscopic anal- 
ysis methods [P2, P6]. In addition, the extracted parameters were strongly cor- 
related with trends observed for electrical permittivity of the sample. Further 
macroscopic material parameters, like chemical barrier properties, tensile 
strength or electrical conductivity are possible design goals of nanocomposite 
manufacturers. The direct correlation of the OCT analysis to these parameters 
could significantly advance material development. 


In addition to sample characterization in a laboratory environment, a dedicated 
extruder probe was developed and demonstrated for OCT-based in-line 
characterization in a fabrication environment. The OCT sample head with 
extruder probe was mounted on a compounding extruder at a pilot plant at the 
Fraunhofer Institute for Chemical Technology (ICT). During nanocomposite 
processing, the dispersion state has been monitored in real time, revealing 
good correlation with the processing parameters [P2]. This technique is not 
limited to compounding extruders but can be applied to a range of composite 
material production methods. 
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Apart from nanocomposites, the analysis of single dispersed nanoparticles 
of different sizes and shapes has been conducted in this work. Making use of 
the size-dependent scattering of nano and micro particles, dispersed polysty- 
rene nanospheres of different diameter and concentration have been character- 
ized using OCT [P2, P3, P5]. By fitting the calibrated single scattering model 
to the measurements based on Mie’s scattering theory, the particle diameters 
have been determined with an accuracy of below 4 % [P2]. In addition, aspect 
ratio and diameter of dispersed gold nanorods have been estimated. For this 
purpose, polarization-sensitive OCT measurements have been fitted to numer- 
ical scattering calculations. The rod diameters have been determined with 
25 % deviation and the aspect ratio matches the nominal value by at maximum 
10 % deviation. The application of this technique to monitor fabrication pro- 
cesses of particles with critical demands on shape or size could lead to greatly 
reduced development cycles of nanoparticles for technical or 
medical applications. 


In addition to fibre-optic OCT systems and material analysis with OCT, the 
integration of OCT on a silicon photonic chip has been demonstrated in this 
work. In the framework of a multi-project-wafer run [38], silicon photonic in- 
tegrated OCT systems have been fabricated at A*STAR IME in Singapore. 
Our demonstrations represent the first fully integrated swept-source OCT on 
silicon, including an interferometer and a balanced photodetector. Two differ- 
ent configurations have been fabricated, one with external reference arm for 
long-distance imaging, and one with internally integrated reference arm allow- 
ing sample measurements at the chip edge. The systems offer up to —64 dB 
sensitivity and enable imaging of biological and technical samples [P1, P4, 
P7-10]. The sensitivity is still below those of fibre-based systems, however, 
improvements in coupling losses and reflections offer the perspective to en- 
hance the sensitivity to at least -90 dB. The on-chip integration brings sub- 
stantial advantages in robustness, cost per unit, and footprint and is hence an 
essential step towards low-cost sensors for industrial applications. Highly in- 
tegrated OCT systems with many parallel channels could even obviate lateral 
mechanical beam scanning, and thus represent a further step towards powerful 
OCT characterization of large sample volumes. 
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A. Mathematical and physical 
definitions 


A.1 Fourier transform 


This section serves as reference for the Fourier transforms used within this 
work and derives important correspondences. All definitions refer to Kiencke 
and Jäkel [134]. The Fourier transform F{-} of a signal y(t) in the time do- 
main to a signal y( f ) in frequency domain has been favoured over the often 
seen transformation to angular frequency @, as it gets along without addition- 
al asymmetric pre-factors. 


Note that temporal and spatial Fourier transforms and their inverse transforms 
assume a time dependency of the form ren). as it is common in electrical 
engineering. In physics, the symbol i=-1 is often used for the imaginary 
unit, and a time dependence of the form exp(-ior) is assumed. Both nota- 
tions describe the same physical situation. The formulations are complex con- 
jugate to each other. 


The Fourier transform and the inverse Fourier transform are defined as 


y(t) = F {y(t)} = fro e?rfigr, 
K (A.1) 
I= FHI = [re ar. 


The concatenation of inverse and forward Fourier transforms leads to the orig- 
inal signal, shown here for the spectrum 


139 


Appendix 


+00 +00 +00 


y(f)= Í y(t ye Prfigr= Í Po Fe ‘af’ e?rfigr 
= Í iF) f € -j2a(f-f')t dt df’ = a jair- f Jay" a 
o(f-f’) 
=y(f) : 


A multiplication in time domain corresponds to a convolution in 
frequency domain 


FOREN [WU nt -P ya". aay 


oo 


Time shift is linked to linearly increasing phase of the spectrum 
F {y(t —to)} =e 7o y(f), (A.4) 


and a linear phase increase in the time domain corresponds to a spectral shift 
in the frequency domain, 


Fl y(t)} =5(F - fo). (AS) 
The cosine-function corresponds to a pair of delta-functions 
1 
Feos(Zrfr))=Z(ölt + fo) + 5(F -= fo) - (A.6) 
A rectangular window function 
1, ees 
rect ,,( ) 2 (A.7) 
At 
0, || > 
2 
corresponds to a sinc function 
z$ n( fAt 
F {recta (t)} =A ZNS zu ) = At sinc( fAr). (A.8) 
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A Hann window that weights a rectangular function with a raised cosine func- 
tion of one period 


Ymar (t) = rectu ol + cos )| (A.9) 


corresponds to three overlapping -shifted sinc-spectra, which partially sup- 
press the side-lobes of each other [135,136] 
At| . 1. 1, 
Ff Yma (r)} = S| sine( 4, + „sine(zfAr + m) + „sine (zfAr _ a)| 
(A.10) 


A Gaussian pulse of 1/e -width of 1/ Va in the time domain transforms into a 
Gaussian spectrum of 1/e -width of Ja i; m 


Fee Te a (A.11) 


A.2 Discrete Fourier transform 


In the discrete Fourier transform (DFT), signals in both domains, here referred 
to as time and frequency domain, are discretized and limited in number of el- 
ements. The discretization in time domain leads to a periodic spectrum, and 
the limited number of observation points in time domain leads to discrete 
sampling points in the frequency domain [134]. The total observation time 


is split into M equal sampling time steps ôty. The sampling frequency 
1 
=—=2A (A.13) 
Ja Bt, Ya 
is the reciprocal of the sampling time 6f, and is twice the observable band- 
width Af;. The frequency resolution is obtained by dividing the sampling 
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frequency fy into M intervals, resulting in the reciprocal of the total 
observation time, 


pp E (A.14) 
M tM Aty 


With counting variables m and i for points in time and frequency domain, 
respectively, the DFT and the inverse DFT are defined as 


M-1 
y= > Ve 2zim/M 
mS (A.15) 


y; el 2rzim/M ‘ 


In this relation, y,, are the elements of the data vector in time domain and Y, 
the elements in frequency domain. 


A.3 Convolution and correlation 


The convolution of two signals is given by 


+00 


yi(t)*y2(t)= f(z) y2(t-z)dr . (A.16) 


00 


In this thesis random processes are assumed to be ergodic, see A.6. The tem- 
poral correlation [44] of two signals is then defined as 


+00 


By (1)= f xi (y(t +r). (A.17) 


00 


A.4 Delta distribution 


The delta distribution [137], or Dirac delta function 6 (x), is defined by its 
measure with a test function ® (x) 
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+00 
(8,0)= | 6(x)®(x)dx=0(0), (A.18) 
yielding the value of the test function at the position x =0. The test function 
®(x) needs to be continuous and indefinitely often continuously differentia- 
ble. The function ®(x) either needs to be zero outside a limited region, or the 
function and all its derivatives need to approach zero for |x| — œ, faster than 
any power |x|” approaches zero. The second case includes the Gauss- 


x 


functione ~ . The delta distribution can be described as follows, 


ö(x)=0, in—œ<x<0 and 0< x< +0, 


a (A.19) 
Í ô (x)dx =I. 
The integral of a delayed delta function with a test function results in the test 
function at the value of the delay, 
+00 +00 
| Elx- xa )®(x)dx= | E(x,- x)D(x)dr=P(x,);  (A.20) 


oo oo 


A.5 Gaussian beams 


A Gaussian beam is a solution of the paraxial Helmholtz equation and com- 
bines wave- and ray-optical elements [138]. It is characterized by a beam 
waist with radius Wo at position Zp and a Gaussian-shaped decay of the com- 
plex field amplitude Ag (7,z) with increasing radial distance. Here, we assume 
linear polarization and Ag (7,z) denotes the complex amplitude of the domi- 
nant transversal field component. Its radial and axial dependence is described 
by 


lt i jkr? 
Ac (r z) = ea Ea e 2Ralz) tes (Ka). (A.21) 


Here, Ag o is the maximum field amplitude, 
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(A.22) 


is the confocal parameter, 


w(z)= nfe) (A.23) 


is the beam radius over axial position z, 


Rg (z)=(z-z) ur 2 a) (A.24) 


Z—Zf 


is the radius of curvature of the phase fronts, and 


(2) zarean 2E) (A25) 


denotes the Gouy phase, which describes a continuous phase shift of the mode 
from - 7/2 to +7/2, when passing the focal point. 


A.6 Coherence of optical fields 


Coherence of optical fields describes the correlation of their field vector com- 
ponents either at different positions in space (spatial coherence) or at different 
times (temporal coherence), or for different vector components qı, q2, with 
q= {x y, z} [139]. The coherence function of the first order 


Gaa (te iito B) = E, (4,7) Ey, (1,5) (A.26) 


is given by the correlation of two complex electric fields E, (t,,%), 


Ep (b, 2) in directions q1, q2, at the points in time t,, f, and at the spatial posi- 


tions 7,7). The normalized complex coherence function 
Eno (ti. f; h» ry) 


as (A.27) 
= Ge (4, h; ty, »)/ 


Ga, nN (ti, Ait, yale (ty, h;t, A) 


A Mathematical and physical definitions 


states, if two fields are coherent to the first order, then | 8a. (t ñ; bs, 2) =], 


or if two fields are incoherent to first order, |8,,,4, tts At, 5) =0. Interme- 


diate values define partially coherent fields. 


Up to now, ensemble averages are considered. If averaging over a sufficiently 
large ensemble of realizations of a random process leads to the same result 
than averaging one realization of the process over time, the process is called 
ergodic. Stationarity is a prerequisite for ergodicity, but is not sufficient [134]. 


In most practical cases, ergodicity can be assumed. Then, not the absolute 
points in time 7, t»), but only the time difference T =t, —t; is relevant and 
Eq. (A.26) results in 


IL, „(er th) = Ey, (6A)Ep (t-15), (A.28) 


which is normalized to 


Yara Eh) EA) En (OFF) Ey, gn (0472-4) 


(A.29) 


Since the coherence function is represented by a correlation of £} (t,,7) and 
Ey, (t, 2); the Wiener-Khintchine theorem [26] can be applied to calculate 
their mutual spectral density T M al f.%.%) by a Fourier transform of 


Eaa (T, 7,7%). This results in 


+00 


Delan anne ir. (A.30) 
Temporal coherence describes the correlation in time at one position in space, 
thus the normalized coherence function is given by y(t, N» ii). The coherence 
time Tm, defines the maximum time difference le] < T, for which | y(r, AN ) z1 
holds. In accordance with this, the coherence length /,=cgt, defines the 
propagation length in which the field varies sinusoidally with constant phase 
offset and amplitude and can thus be assumed to be one longitudinal 
mode [44]. The reciprocal of the coherence time is the spectral width 


of =Vt,. 
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When spatial coherence is investigated, the time difference is set to r=0, and 
the spatial coherence function results, (0, N, h). The coherence radius 
T, S li = ñ| gives the maximum spatial distance of 4,7 for which 
lv (0, N, 7 )| x1. From this, a coherence area F, =ar can be defined. 
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C. Glossary 


C.1 List of acronyms 


A-Scan One-dimensional (1D) OCT depth scan 
Al,2 Anode of a diode with contacts 1,2 
ACF Autocorrelation function 

ADC Analogoue-to-digital converter 

AF Area fraction 

AFM Atomic force microscopy 

APC Angled physical contact connector 
AR Aspect ratio 

AWG Arrayed waveguide grating 

B-Scan 2-D OCT cross sectional scan 

BD Balanced detector 

BL Ball lens 

BS Beam splitter 

C-Scan 3-D OCT volumetric scan 

C1,2 Cathode of a diode with contacts 1,2 
CCD Charge-coupled device 

CMOS Complementary metal-oxide-semiconductor technology 
CMRR Common-mode rejection ratio 

CNT Carbon nanotube 

CPL Optical 3-dB coupler 

DDA Discrete-Dipole Approximation 
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DFT 
DLS 

DR 

FC 
FD-OCT 
FT 
FWHM 
G1, G2 
HWP 
IDFT 
ITO 

LF 

LM 
LP-modes 


MMI 
MoM 
MPW 
MWCNT 
NEP 
NIST 
OCT 
OCText 
OCTint 


PA 
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Discrete Fourier transform 

Dynamic light scattering 

Dynamic range 

Fibre collimator 

Fourier-domain optical coherence tomography 
Fourier transform 

Full width at half the maximum 

Gold nanorod samples in aqueous dispersion 
Half-wave plate 

Inverse discrete Fourier transform 

Indium tin oxide 

Lensed fibre 

Light microscope 


Linear polarized modes in fibres or cylindrical 
waveguides 


Multimode interference 

Method of Moments numeric field calculation approach 
Multi-project-wafer 

Multi-wall carbon nanotube 


Noise-equivalent power 


National Institute of Standards and Technology 
Optical coherence tomography 
Integrated OCT system with external reference arm 


Integrated OCT system with internal integrated 
reference arm 


Polyamide 
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PAR 
PAR 
PBS 
PC 

PD 
PIC 
PL 
PMF 
Pol 
PolC 
PP 

PS 
PS-BD 
PWC 
QWP 
RDG 
RIN 
RMS 
RP 

S1, $2, S3, S4 


SAXS 
SD-OCT 
SEM 
SESF 


Perimeter-to-area ratio 

Modified pixel-wise PAR 

Polarization beam splitter 

Personal computer 

Photodiode 

Photonic integrated circuit 

Polymer microlens 
Polarization-maintaining fibre 
Polarizer 

Polarization controller 

Polypropylene 

Polystyrene 

Polarization-sensitive balanced detector 
Partial wave spectroscopic microscopy 
Quarter-wave plate 
Rayleigh-Debye-Gans Theory 
Relative intensity noise 

Root mean square 

Reference path 


Nanoparticle samples with different concentration and 
nanosphere diameter 


Small-angle X-ray scattering 
Spectral-domain optical coherence tomography 
Scanning electron microscopy 


Spectral encoding of spatial frequency 
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SLD 
SLS 
SMF 
SNR 
sol 

SP 

SS 
SS-OCT 
TD-OCT 
TE mode 
TEM 
TIA 
TRPS 
XRD 


Superluminescent diode 

Static light scattering 

Single-mode fibre 

Signal-to-noise ratio 

Silicon-on-insulator 

Sample path 

Swept-laser source 

Swept-source optical coherence tomography 
Time-domain optical coherence tomography 
Transverse electric mode 

Transmission electron microscopy 
Transimpedance amplifier 

Tunable resistive pulse shaping 


X-ray diffraction 


C.2 List of symbols 


Greek symbols 


Symbol Unit 

a Orientation angle of the electric field rad 

4 Strain amplitude 1 

a Normalized coherence function of ergodic | 1 
fields in spatial directions q1, 9 

T, (7 ) Complex coherence function of the Vim 
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source field 
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ô fa 
ôP, (t) 


or, Ôx, Ôy 
Öl 


Coherence function of ergodic fields in 
spatial directions q1, q2 


Mutual spectral density of two fields in 
spatial directions q1, q2 


Dirac delta distribution 
Common phase in Jones vector 


Difference in parasitic reference power to 
nominal reference power (OCTex?) 


Instantaneous source bandwidth 
Frequency step size after sampling 


Mean-free optical power fluctuations 


öR(t)=R(t)-A(t) 
Lateral resolution of an OCT system 


Sampling-related resolution in sample re- 
turn time 


Measurement resolution of the sample 
response time 


Quantization voltage step size 


Axial (depth) resolution of an OCT sys- 
tem 


Total ratio of measurement and simulation 
in particle shape determination 


Optical measurement bandwidth in 
FD-OCT 


Electrical measurement bandwidth 
Total frequency range after sampling 


Unambiguity range in sample return time 
due to sampling 


V’ m? 


V?’s/m? 


rad 


dB 
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At 


Ag 


E0 


EAR 


Ea 


166 


Temporal width of an exemplary 
rectangular function 


Coherence limited measurable time range 
of sample responses 


Coherence related imaging depth 
Focal measurement range 
Wavelength span of an OCT light source 


Ratio of simulated and measured scatter- 
ing cross section, polarization maintaining 
backscattering cross section and polariza- 
tion cross talk, respectively 


Optical phase difference of light returning 
from sample and reference arm 


Vacuum permittivity 
= 2 -12 
£o = (wos) »8.854x10°"? As/(Vm) 


Relative accuracy in particle aspect 
ratio determination 


Accuracy in particle 
diameter determination 


Relative permittivity 


Common-mode power rejection ratio of a 
balanced detector (linear) 


Common-mode power rejection ratio of 
on-chip receiver (dB) 

Focal power response function of an 
OCT system 


Polarization cross talk in 
PS-OCT measurements 


5 


dB 


dB 
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5g; I Rotation angle of the fast axis of the QWP | rad 
Ch) and of the transmission axis of the 
polarizer (9p) with respect to the 
horizontal axis 

Pyy (7) Temporal correlation of the signals y,, y>- 

K(z) Real-valued normalized baseband coher- |1 
ence function of the source field 
Optical wavelength in vacuum, Å = co / f |m 

= Central wavelength of an OCT m 

light source 

Ho Magnetic permittivity, N/A? 
Ly = 42x10! N/A? 

My, m> Hb, c Polarization maintaining (m) and cross- m! 
polarized (c) backscattering coefficient 

L, Relative magnetic permittivity 1 

Leb Extinction 44 and backscattering m! 
coefficient 4w 

v Factor describing fractions of 1 
a particle diameter 

p Responsivity of the photodiode A/W 

PBG.dB Backscatter floor of OCText dB 

PWG.dB Distributed on-chip waveguide backscat- |dB 
ter within resolution z 

Oa,s,t,b Absorption- O, total scattering- o,, m? 
extinction- o;, and backscattering 
cross section Op 

Opnm> Ob, Polarization maintaining (m) and cross- m? 
polarized (c) backscatter cross section 
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T Time difference s 
(A Reference mirror return time s 
To Coherence time s 
o Relative phase in Jones vector rad 
ON Constant phase offset in an rad 
SS-OCT signal 
PG (z) Gouy phase in a Gaussian beam rad 
(0) (x) Arbitrary test function on x, indefinitely 
often differentiable and outside a limited 
region identically zero. 
o Optical angular frequency, © =2z7f sg! 
O) Optical angular carrier frequency of the s” 
light source 
Latin symbols 
Symbol Unit 
a Parameter of a Gauss function s” 
dap Lensed port power loss (OCText) dB 
A Electric field amplitude V/m 
A, Ay ay Electric field amplitude vector and x,y,z | V/m 
components thereof 
Ay Slowly-varying random amplitude ofthe | V/m 
source field 
A, (r,z) Complex field amplitude of a Gaussian | V/m 
beam with maximum Ago 
A; Individual agglomerate area m’ 
Atot Imaging cross section m 
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w 8,8 & 
© 


l 
© 


Er Gee 


Zout? Zout,m? =out,c 


Confocal parameter of a Gaussian beam 


Mechanical oscillation amplitude of 
a rheometer 


Magnetic induction field vector 


Vacuum speed of light, 
Co =299 792 458 m/s 


Diameter of a scatterer 


Nanorod diameters of samples G1 
and G2 


Electric displacement field vector 

Unit vectors 

Electric field vector 

Electric field emitted by the light source 


Electric field measured at 
the photodetector 


Electric field vector of polarized incident 
electric field 


Electric field vector of polarized output 
electric field E,,,, polarization main- 
taining (m) and cross-polarized (c) 


Electric field reflected by the reference 
mirror, measured after second pass of the 
beam splitter 


Electric field backscattered from the 
sample, measured after second pass of 
the beam splitter 


Frequency 


Auxiliary frequency variable, replacing f 


V/m 


V/m 


V/m 


V/m 
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to Optical carrier frequency of the Hz 
light source 

Fa Baseband frequency in Hz 
RIN measurement 

fa Sampling frequency Hz 

fe Electrical carrier frequency Hz 

Ta Centre frequency of a FD-OCT Hz 
frequency span 

Fig OCT signal frequency Hz 

Fo Electric field-covered area of the source | m? 


emitted beam 
F Coherence area m 


Sia Normalized coherence function of fields | 1 
in spatial directions q1, q2 


Coherence function of the first order of V?/m? 
fields in spatial directions q1, q2 


h(t), h(t) Field impulse responses of sample and |s 
reference path 


h(f).h,(f) Complex field transfer function of sam- | Hz | 
ple and reference path 


H Magnetic field vector A/m 
i Counting variable, e. g. in a DFT 1 
iy OCT photocurrent related to sample- A 


reference interference 


I, Effective OCT photocurrent related to V’/m 
sample-reference interference, 
i, = 2Rel E,E,| ~ i, 
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Lot vers 


1135, Tz, Iız 


Í hors 


Is, 


Fourier transformed effective 
OCT photocurrent 


Photocurrent in general or from PD 1 or 
PD 2, respectively 


Maximum photocurrent 
Total noise current 
Quantization noise current 
RIN-current 


OCT signal current after 
balanced detection 


Equivalent photocurrent according to 
ADC or PD limitation 


Shot noise current 
Detector noise current 


Light intensities depending on polariza- 
tion: total Zot horizontal Jor, vertical Zer, 
at 45° I4s, at 135° 1135, right-circular Irz, 
and left-circular /ıLz 


Intensity of a Gaussian beam with max- 
imum J, , at the focus 

Imaginary unit, j=v-1 

Jones vector and x,y components thereof 


Jones vector exiting the fibre towards the 
sample (Jin) and entering the 

PS-BD (J, ) 

Jones vectors towards sample (s) and to- 


wards reference (r) arm in the free-space 
PS-OCT setup in Figure 2.12 


V?’ m? 


> 


> > >>> 


W/m? 


W/m? 


jean 


= 
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Srey 


Jowe: Jnwe- Spot 


Js Jys Jm 


Po 
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Jones matrix; transformation of 
Jones vectors 


Jones matrix of the reference arm in the 
free-space PS-OCT setup in Figure 2.12 


Jones matrix of reverse propagation 
through any optical reciprocal material 


Jones matrices of a QWP, a HWP, a Po- 
larizer (Pol), a sample (s), a fibre optic 
system (sys) and a plane mirror (M) 


Angular optical wavenumber, k = 277/A 


Optical wave vector and x,y,z 
components thereof 


Real-valued normalized coherence func- 
tion of the source field 


Real-valued coherence function of the 
source field 


Coherence length 
Nanorod lengths of samples G1 and G2 
Counting variable, e. g. in a DFT 


Jones matrix elements of an optical re- 
ciprocal material 


Number of sampling points in a DFT 


Refractive index. Refractive index of 
sample n,, and reference path n,. 


Quantization error of a voltage signal 
Volume number density 


Backscattering probability 
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Pin 
P P 


out? ~ out,m? ~ out,c 


Sr 


Probability density of quantization noise 
inside quantization interval 


Optical power 
Optical power emitted by the source 


Optical power scattered back from 
the sample 


Optical power at photodetector 
Optical power incident on the sample 


Optical output power (Pout), with main- 
tained polarization (m) and 
cross-polarized (c) 


Optical power incident on one detector 
photodiode, which is returned from the 
sample path or the reference 

path, respectively. 


Material polarization vector 


Decay calibration factor in an 
OCT measurement 


Elementary charge, gy = 1.602 x 10° C 


Power calibration factor in an 
OCT measurement 


Radius with respect to the optical axis 


Spatial vector, e.g. propagation direction 
of an electric field 


Coherence radius 


Amplitude reflectivity depth-profile of 
the reference path 


= 


As/m? 
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Ry(z) 


Rap(Z) 
R,(z) 


Rs 
Ro Ry» Ru» R 


xx? “yx? “xy? “yy 


RIN(f - f.) 


RINag (f = fe) 
RIN 


tot 
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Amplitude reflectivity of the 
reference mirror 


Amplitude reflectivity depth-profile of 
the sample 


Rotation matrix, with rotation angle 9 


Theoretical backscattering factor from a 
scattering sample, excluding influences 
of the measurement system 


Semi-logarithmic OCT backscatter signal 


Radius of curvature of the phase fronts of 
a Gaussian beam 


Backscatter factor with maintained polar- 
ization (m) and cross-polarized (c) 


Noise floor in an OCT measurement, ex- 
pressed in minimum measurable power 
reflectivity of a sample 


Power reflectivity of the reference mirror 


Depth-dependent power reflectivity pro- 
file of the sample 


Power reflectivity of a sample reflector 


Backscattering coefficient R,; denoting 
scattering of i- into j-polarization 
Spectral RIN density around carrier 
frequency f, 

RIN in 1 Hz bandwidth around carrier f, 


Total RIN integrated over 
entire spectrum 
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A 


T, 


sw 


Ug, Ug 


Uy, U, 


Umax,PD ’ U max, ADC 


Usig BD 


Spatial RIN profile around reflection 
at Ze 


Spatial RIN within OCT depth 
resolution dz 


Rheometer plate distance 
Perimeter of a single agglomerate 
Tuning slope of a swept-source laser 


Sensitivity of the OCT system, expressed 
in sample power reflectivity 


Poynting vector 


Stokes vector describing the polarization 
of light 


Time 

Frequency ramp start time 
Reference path return time 
Sample path return time 
Absolute temperature 
Swept-source ramp duration 


Input signal u, and output signal u, of a 
quantizer. Here: voltage signals 


Matrix elements of a unitary 
Jones matrix 


Voltage saturation limits of photodetec- 
tor and ADC, respectively. 


OCT signal voltage after 
balanced detection 
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Ze 
Zedge,int 
Z 


edge,ext 


Zf 
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Speed of light in sample and reference 
arm, respectively 


Waist radius of a Gaussian beam 


Radius of a Gaussian beam over 
axial position 


Spectral windowing function 
Temporal windowing function 


Lateral spatial coordinate in an 
OCT measurement 


Arbitrary variable in 
exemplary functions 


Lateral spatial coordinate in an 
OCT measurement 


Arbitrary time-dependent function 


Correspondent function to y(t) , after 
Fourier transform, with dependence on 
frequency f. 


Time-data vector element in a DFT 
Frequency-data vector element in a DFT 
Depth axis of an OCT measurement 
Depth position of a sample mirror 

RIN probing depth 

Depth position of a reflection peak 
Facet position of OCTint chip 

Facet position of OCText chip 


Beam waist position of a Gaussian beam 


B 5B BB B E B 
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Depth axis of an OCT measurement 
(geometrical depth) 


Depth position of the beam splitter 
Depth position of the reference mirror 
Depth position of the sample 
Impedance of free space 


Transimpedance gain of the 
photodiode amplifier 


5 


Oo 8 BaB 


V/A 
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